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ABSTRACT 

Wu,  I-pai.   Ph.D.,  Purdue  University,  June  1963.   Hydrology 
of  Small  Watersheds  in  Indiana  and  Hydrodynamics  of  Overland  Flow. 
Major  Professor:  J.  W.  Delleur. 

A  study  is  made  of  the  hydrology  of  watersheds  less  than  200 
square  miles  in  area  located  in  the  State  of  Indiana  for  which  flows 
are  recorded  by  the  U.S.G.S.  A  statistical  frequency  analysis  of  the 
peals  flows  was  made  by  means  of  Gumbel's  extreme  value  theory.  A 
geomorphological  study,  involving  the  quantitative  determination  of 
five  watershed  characteristics,  combined  with  the  use  of  multiple 
correlation  techniques, serves  to  establish  a  formula  for  estimating 
the  peak  discharge  of  ungaged  watersheds  in  Indiana. 

The  shape  and  the  peak  discharge  of  the  synthetic  hydrograph  were 
determined  for  those  areas  where  no  stream  gaging  station  is  available. 
A  mathematical  expression  for  the  hydrograph  containing  two  parameters 
which  can  be  correlated  with  readily  obtainable  watershed  characteris- 
tics gives  the  theoretical  basis  for  establishing  the  synthetic  hydro- 
graph.  A  complete  procedure  for  the  design  of  the  storm  hydrograph 
for  small  ungaged  watersheds  is  presented.  The  basic  mathematical 
expression  for  the  hydrograph  is  considered  appropriate  for  general 
application. 
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The  research  concludes  with  a  study  of  the  hydrodynamics  of 
overland  flow  based  on  the  momentum  and  the  continuity  equations.  A 
set  of  two  partial  differential  equations  of  hyperbolic  type  was  derived 
for  the  overland  flow  produced  by  the  precipitation.  The  equations  may 
be  solved  by  the  method  of  characteristics.  Approximate  solutions 
also  were  given  for  both  the  laminar  and  turbulent  sheet  flow.  They 
are  simple  to  calculate  and  the  results  compare  well  with  Izzard's 
experimental  measurements. 


PART  I 
DESIGN  PEAK  DISCHARGE 
FOR 
SMALL  WATERSHEDS  IN  INDIANA 


INTRODUCTION 

The  determination  of  the  required  waterway  area  of  a  bridge  or 
the  selection  of  the  size  of  a  culvert  are  problems  which  require  an 
accurate  estimate  of  the  peak  flood  discharge  that  will  pass  through 
the  structure.  The  determination  of  this  discharge  is  more  difficult 
for  small  watersheds  because  the  majority  of  them  are  ungaged.  There 
are  very  few  gaged  small  watersheds  on  which  to  base  an  estimate.  Par- 
ticularly in  the  state  of  Indiana,  there  is  very  little  information  on 
watersheds  less  than  200  square  miles.  There  are  only  twelve  water- 
sheds of  less  than  100  square  miles,  and  seventeen  watersheds  with  an 
area  between  100  and  200  square  miles  for  which  the  U.S.G.S.  is  cur- 
rently reporting  flows.  A  research  program  was  initiated  at  Purdue 
University  to  obtain  reliable  methods  for  estimating  the  peak  dis- 
charge for  a  safe  and  economic  design  of  highway  drainage  structures, 
serving  watersheds  in  Indiana  of  less  than  200  square  miles  but  larger 
than  20  square  miles. 

The  existing  methods  are  empirical  and  fail  to  take  into  account 
the  factors  upon  which  the  runoff  depends.  Ki  unison  in  I9U6  and  Chow* 
in  1962  have  given  a  complete  list  of  empirical  formulas  which  include 
the  watershed  characteristics.  The  most  frequently  used  formulas  are 

those  of  Talbot*5'  published  in  1887  and  Meyer *  '  in  1879  and  the 

( 5) 
Rational  formula  originally  derived  by  Mulvanyw'  in  1857-  Talbot's 

formula  was  originally  intended  for  locations  in  Illinois.  It  estimates 


the  waterway  area  from  the  watershed  area.  The  formula  is: 

A  =  CM5/  (1) 

where  A  is  the  required  waterway  area  in  square  feet,  M  is  the  watershed 
area  in  acres,  and  C  is  a  coefficient  varying  between  1/5  and  1  depend- 
ing on  the  slope  and  character  of  the  watershed.  The  selection  of  the 
coefficient  depends,  among  other  things,  on  the  experience  of  the 
designer.   Due  to  the  various  factors  that  affect  the  runoff  other  than 
the  watershed  area,  the  value  of  the  coefficient  C  cannot  be  accurately 
determined  to  represent  all  the  watershed  characteristics.  Talbot's 
formula  is  unsatisfactory  for  a  safe  design  of  a  hydraulic  structure. 

Yulev   in  1950  developed  a  similar  formula  for  locations  in 
Indiana.   The  required  waterway  area  is  expressed  as  a  function  of  the 
two- thirds  power  of  the  drainage  area.   It  is 

A  =  Ca2/5  (2) 

where  A  is  the  waterway  area  in  square  feet,  a  is  drainage  area  in 
square  miles,  and  C  is  a  general  coefficient  which  varies  with  the 
watershed  topography:   C  is  0.3  to  0.7  for  flat  land,  0.7  to  1.3  for 

rolling  land,  and  1.3  to  2.2  for  hilly  land. 

(7) 
Benson    in  1959  found  that  the  mainstream  slope  is  next  in 

importance  to  drainage  area  among  the -.factors  which  affect  the  runoff. 

The  slope  considered  was  for  that  part  of  the  mainstream  located  between 

85  to  10  percent  of  the  total  distance  above  the  gaging  point.  The  fol 

lowing  empirical  formula  wa6  found  for  the  New  England  region; 

Q=a.AbSC  (3) 

where  Q  is  peak  discharge  in  cfs,  A  is  drainage  area  in  square  miles, 


S  is  the  85  to  10  percent  slope  of  the  mainstream,  and  a,  b,  c,  are  the 
regression  coefficients  varying  with  the  recurrence  interval  of  flood. 

The  factors  affecting  the  runoff  may  be  grouped  into  three 
categories:   the  storm  characteristics,  the  geomorphological  and  the 
geological  characteristics  of  the  watersheds.  The  dependence  of  the 
runoff  on  the  geomorphology  of  the  watershed  is  analyzed  in  this  study. 
The  effect  of  the  storm  characteristics  and  of  the  soil  types  are  con- 
sidered in  the  following  study  entitled  "Design  Hydrographs  for  Small 
Watersheds  in  Indiana . "    The  validity  of  the  unit  hydrograph  theory 
will  ultimately  be  analyzed  by  a  theoretical  study  of  overland  flow 
hydrodynamics . 


STATISTICAL  ANALYSIS 

The  following  two  methods  give  a  linear  relationship  between  the 
observed  variable,  i.e.,  the  discharge,  and  the  recurrence  interval  or 

a  function  of  the  recurrence  interval. 

(8 ) 
Geyer,  in  19*K),  derived  a  mathematical  expression  relating  the 

flood  magnitude  to  the  exceedance  interval 

Y  »  LR  X  (M 

where  Y  is  the  flood  magnitude  in  sec-ft.  that  has  an  exceedance  in- 
terval of  t  years,  and  L,  R,  K  are  constants  for  a  particular  stream. 
Taking  logarithms  twice  on  both  sides  of  Eq.  k,   one  obtains 

log  (log  L  -  log  Y)  »  K  log  t  +  log  log  R  (5) 

where  L  is  then  determined  by  trial  and  error  so  that  Eq.  5  represents 
a  straight  line  when  log  (log  L  -  log  Y)  is  plotted  against  log  t. 

Gumbel,   '  10'  ^'   in  19^1,  developed  the  extreme  value  theory 
which  is  used  to  analyze  the  observed  extremes  and  to  forecast  further 
extremes.  The  theory  states  that  the  probability  *(x)  for  the  dis- 
charge x  to  be  the  largest  among  n  independent  observations  is  given  by 

*(x)  -  e_e  X  -  exp  (-e"y)  (6) 

as  n  ■*  00 
where  e  is  the  base  of  Napierian  logarithms  and  y,  termed  the  reduced 
largest  value,  is  given  by 


y  -  an(x  -  un)  (7) 

where  aQ,  uq  are  two  extreme  parameters,  u  is  a  certain  expected 

largest  value  having  the  return  period  n,  and  its  probahility 

F(uQ)  is  1  -  -  ;  aQ  is  defined  as  nf(u  )  where  f(u  )  is  the  initial 

distribution  given  by  f(u  )  =  F'(u  ).   A  probability  paper  designed 

(12) 
for  extreme  value  was  proposed  by  Powell;  '        The  observed  variate  x 

is  traced  on  the  ordinate,  the  largest  value  y  is  traced  on  the  abscissa, 

both  in  linear  scales.  The  value  of  the  probabilities  «>(x)  is  given  in 

Eq.  (6)  or  by 

*<*>  •  irh.    •  <8> 

*<!•  (8)  gives  the  "plotting  position"  where  m  is  the  rank  of  the 
yearly  maximum  in  increasing  order,  and  n  is  the  number  of  years 
of  observation.  The  return  period  is  given  by 

T-  rrWr  •  (9) 

Both  the  probability  *(x)  and  the  return  period  T  are  laid  off  on 
auxiliary  horizontal  scales  on  the  probability  paper.  Eq.  (7)  relating 
the  observed  variate  x  to  the  reduced  largest  value  y  may  be  written 
as 

x  ■  —  y  +  u  (10) 

a       n 

from  which  it  appears  that  there  is  a  linear  relationship  between  x  and  y. 
Thus,  theoretically,  the  plot  of  x  vs.  y  should  be  a  straight  line  on 
probability  paper. 


(13) 

Benson/  "  in  1950,  made  use  of  historical  data  in  flood- 
frequency  analysis.  The  study  was  made  on  the  basis  of  annual  peak 
discharges,  and  data  were  plotted  on  probability  paper  based  on  Gumbel's 
extreme  value  theory. 

The  present  study  includes: 

1.  A  statistical  analysis  of  existing  peak  discharge  data  and  the 
determination  of  the  25-year  flood  for  32  gaged  small  watersheds 
by  Gumbel's  extreme  value  method. 

2.  A  geomorphological  analysis  of  16  of  these  small  watersheds. 

3.  A  multiple  correlation  of  the  25- year  flood  and  the  geomorphological 
characteristics.  This  is  based  on  the  assumption  that  the  State  of 
Indiana  is  an  area  sufficiently  homogeneous,  so  that  variations  of 
storm  characteristics  and  of  physical  properties  of  soils  are  not 
significant  variables  compared  to  the  variables  describing  the 
geomorphological  characteristics.  The  watersheds  considered  in 
this  study  have  an  area  between  20  and  200  square  miles  approxi- 
mately. These  watersheds  are  called  small  because  they  refer  to 

the  smaller  group  for  which  flow  records  are  reported  by  the  U.S.G.S. 
They  are  large  enough,  however,  so  that  the  land  use  or  the  type  of 
vegetation  is  not  an  important  variable,  and  may  be  disregarded. 
This  method  has  an  advantage  that  the  estimation  of  design  dis- 
charges is  based  on  the  discharge  records  themselves,  thus  eliminating 
the  need  of  relating  rainfall  and  runoff  by  means  of  a  number  of 
variables  difficult  to  evaluate. 


Flood  Frequency  Analysis  for  Small  Watersheds  in  Indiana 
Thirty- two  watersheds  distributed  throughout  the  whole  state 
were  selected  to  study  for  flood  frequency  analysis.  Fig.  1  is  a  map 
of  Indiana  showing  the  thirty- two  watersheds,  and  Table  1  lists  the 
names  of  the  watersheds,  their  assigned  number  and  their  areas.  For 
convenience,  these  numbers  will  be  used  in  this  report  instead  of  the 
name  of  the  watersheds.  There  are  only  three  watersheds,  the  area  of 
which  is  less  than  50  square  miles,  12  watersheds  under  100  square  miles, 
and  29  watersheds  under  200  square  miles  for  which  flow  records  are 
available  from  the  U.S.G.S.   All  of  these  were  included  in  this  fre- 
quency analysis.   In  order  to  have  a  good  distribution  of  the  water- 
sheds over  the  whole  state,  three  additional  larger  watersheds  were 
included.  A  bar  diagram  was  plotted  in  Figure  2  to  show  the  time 
period  of  records  for  each  watershed. 

Instantaneous  annual  peak  flow  was  used  for  the  extreme  value 
analysis  instead  of  the  maximum  mean  daily  flow,  since  it  is  the  peak 
discharge  that  is  desired  for  hydraulic  structures  design.  Data  for 
this  instantaneous  annual  peak  discharge  were  obtained  from  the  U.S.O.S. 
office  in  Indianapolis,  Indiana.  The  data  are  analyzed  on  a  water  year 
basis — that  is,  from  October  to  September.  The  data  were  examined  for 
man-made  changes,  such  as  new  reservoirs,  soil  and  water  conservation 
works,  and  changes  of  gaging  sites.   It  was  found  that  no  correction 
was  needed. 

The  historical  floods  were  Included  in  the  analysis  as  suggested 

(!*>)  (Ik) 

by  Benson,    '   Dalrymple,    '  in  order  to  extend  the  short  time  records 

available.  The  recurrence  interval  of  historical  floods  may  be  computed 

for  the  following  two  cases: 


Number 


Table  1 
List  of  Watersheds,  Their  Area  and  Assigned  Number 

Watershed  ~   ZZ  Watershed 

Gaging  Station  Area 
(sq.mi. ) 

S-  1    Bean  Blossom  Creek  at  Do Ian,  Ind.  100 
S-  2    Clifty  Creek  at  Hartsville,  Ind.  88.8 

8-  3    North  Fork  Vernon  Fork  near  Butlersville,  Ind.     87.3 
S-  k        Hart  Ditch  at  Munster,  Ind.  69.2 

S-  5    Salt  Creek  near  McCool,  Ind.  78.7 

S-  6    Little  Calumet  River  at  Porter,  Ind.  62.9 

S-  7    Cedar  Creek  at  Auburn,  Ind.  93. 0 

S-  8    West  Creek  near  Schneider,  Ind.  U6.3 

8-  9    Iroquois  River  at  Rosebud,  Ind.  30.3 

S-10    Bice  Ditch  near  8outh  Marion,  Ind.  22.6 

S-ll    Big  Slough  Creek  near  Collegevllle,  Ind.  8^.1 

8-12    Carpenter  Creek  at  Egypt,  Indiana  H8.1 

1  Tippecanoe  River  at  Osvego,  Ind.  115 

2  Mississineva  River  near  Ridgeville,  Ind.  130 

3  Wildcat  Creek  at  Greentovn,  Ind.  162 
k          Cicero  Creek  near  Arcadia,  Ind.  131 

5  Fall  Creek  near  Fortville,  Ind.  172 

6  Eagle  Creek  at  Indianapolis,  Ind.  179 

7  Young  Creek  near  Edinburg,  Ind.  109 

8  Blue  River  at  Carthage,  Ind.  187 

9  Sand  Creek  near  Brewer sville,  Ind.  1% 

10  North  Fork  Salt  Creek  near  Belmont,  Ind.  120 

11  Patoka  River  at  Jasper,  Ind.  257 

12  Busseron  Creek  near  Carlisle,  Ind.  228 

13  East  Fork  White  Water  River  at  Richmond,  Ind.  123 
Ik          Silver  Creek  near  Seller eburg,  Ind.  188 

15  Big  Indian  Creek  near  Corydon,  Ind.  129 

16  Kankakee  River  near  North  Liberty,  Ind.  152 

17  Singleton  Ditch  at  Schneider,  Ind.  122 

18  Deep  River  at  Lake  George  Outlet  at  Robert,  Ind.  125 

19  Pigeon  Creek  at  Hogback  Lake  Outlet  near 

Angola,  Ind.  102 

20  Laughery  Creek  near  Farmers  Retreat,  Ind.  2U8 
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FIG.     I-         LOCATION      OF      GAGING       STATIONS        FOR        STUDIED 
WATERSHEDS 


Fig.  2      PERIOD    OF  RECORD    OF    INSTANTANEOUS 
ANNUAL    PEAKS     AT    GAGING    STATIONS. 

Watershed 
NO. 

PERIOD    OF   RECORD 

* 

n 

1900 
1910 

1920 

1930 

1940 

1945 

1950 

1955 

o 

ID 
01 

S-  1 

14 

S-2 

12 

S-3 

18 

S-4 

17 

S-5 

15 

S-6 

15 

S-7 

17 

S-8 

9 

S-9 

II 

S-IO 

II 

S-ll 

10 

S-  12 

10 

1 

_ 

II 

2 

13 

3 

16 

4 

5 

5 

18 

6 

1 

21 
17 

7 

8 

— 

10 

9 

12 

10 

| 

14 

II 

1 

1 

12 

12 

16 

13 

, 

10 

14 

5 
17 

15 

16 

9 

17 

II 

18 

13 

19 

14 

20 

19 

*    n=    Length    of   Record    in    Years 

Table  1* 
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Example  of  Computation  for  Frequency  Study 
Little  Calumet  River  at  Porter,  Indiana 


Year 

Instantaneous  Annual 
Peak  Discharge 
Q  cf  s 

Rank 

Q 
cfs 

* 

y 

19*5 

2,1*1*0 

1 

1*90 

0.0625 

-1.01979 

19l*6 

715 

2 

521 

0.1250 

-0.75210 

I9H7 

2,lUo 

5 

690 

0.1875 

-0.51520 

19*8 

1,960 

k 

715 

0.2500 

-0.52663 

19*9 

690 

5 

81*8 

0. 5125 

-0.15111* 

1950 

1,720 

6 

1,060 

0.5750 

0.01956 

1951 

1,560 

7 

1,170 

0.1*575 

0.19051* 

1952 

1,060 

8 

1,560 

0.5000 

0.56651 

1955 

521 

9 

1,570 

O.5625 

0.55275 

195* 

1,170 

10 

1,1*20 

0.6250 

0.75501 

1955 

5,110 

11 

1,720 

0.6875 

0.98165 

1956 

1,570 

12 

1,960 

0.7500 

1.21*590 

1957 

8U8 

15 

2,lU0 

0.8125 

1.57196 

1958 

1*90 

H* 

2,1*1*0 

0.8750 

2.0151*2 

1959 

1,1*20 

15 

5,no 

0.9575 

2.71*065 
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GEOMORPHOLOGICAL  STUDY 

(15) 
Langbein  and  others  Jl   in  I9U7  Indicated  that  river  floods  were 

the  results  of  many  causes.  One  of  the  primary  objectives  of  scientific 
hydrology  is  the  segregation  and  evaluation  of  the  causative  factors. 
The  climatic  factor  and  the  soil- vegetation  complex  are  variables  that 
exercise  their  principal  influence  on  the  volume  of  runoff.  The  to- 
pography of  drainage  basins  is  a  reasonably  permanent  characteristic 
which  influences  mainly  the  concentration  of  time  distribution  of 
discharge  from  a  drainage  basin.  The  topographic  characteristics  of 
drainage  basin  cited  were:  area  of  basin,  drainage  density,  area- 
distance  distribution,  length  of  basin,  land  slope,  area-altitude 
distribution,  and  area  of  water  surfaces. 

Strahler,  *   '  in  1952,  expressed  the  area-altitude  relation  by 
a  hypsometric  analysis.  The  hypsometric  curve  relates  horizontal 
cross- section  area  of  a  drainage  basin  to  relative  elevation  above 
basin  mouth.  By  the  use  of  dimenslonless  parameters,  the  curves  can 
be  described  and  compared  irrespective  of  true  scale.  The  area  under 
the  curve  can  be  used  to  find  total  land  mass  and  the  mean  relief  of 
watersheds. 

Strahler,  '  ''  in  1957,  made  a  quantitative  analysis  of  watershed 
geomorphology  and  showed  that  the  linear  scale  measurements  include 
length  of  stream  channel  of  given  order,  drainage  density  and  relief. 
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Surface  and  cross- sectional  area  of  basins  are  length  products.  Dimen- 
sionless  properties  include  stream  order  numbers,  stream  length  and 
bifurcation  ratio,  maximum  valley- side  slopes,  mean  slopes  of  water- 
shed surfaces,  channel  gradients,  relief  ratio  and  hypsometric  curve 

properties. 

(7) 
Benson,     in  1959,  cited  the  following  basin  characteristics: 

drainage  area,  channel  slope,  land  slope,  tributary  channel  slope, 
watershed  shape  factor,  mean  elevation,  percentage  of  lakes,  swamps 
and  reservoirs,  and  drainage  density.  The  shape  factor  of  the  water- 
shed was  represented  by: 

L2 

j—     ,  equivalent  of  basin  length  divided  by  basin  width. 

•r       ,  basin  length  divided  by  drainage  area. 

I  aL,   the  summation  of  small  subdivisions  of  the  drainage  area, 
each  multiplied  by  the  distance  of  travel  to  the 
gaging  point. 

In  this  study,  the  first  geomorphological  factor  considered  is 
the  drainage  area.   It  is  the  projected  area  of  watershed,  also  called 
the  catchment  area.  Obviously,  the  bigger  the  area  of  catchment  the 
larger  the  amount  of  runoff.  However,  the  rate  of  runoff  is  largely 
dependent  on  the  slope  of  the  land,  the  drainage  density,  and  the 
slope  of  the  streams.  A  steep  land  slope  with  high  drainage  density 
and  larger  slope  of  stream  will  give  a  higher  rate  of  runoff  than 
those  which  have  smaller  slope  and  lower  drainage  density.  The  shape 
of  the  watershed  is  also  a  factor  affecting  the  runoff,  because  it 
affects  the  time  of  concentration.  With  a  three-dimensional  concept  in 
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mind,  the  principal  geomorphological  factors  which  can  affect  the 
amount  of  peak  discharge  are  listed  as  following. 

1.  The  drainage  area 

2.  The  drainage  density 

3.  The  land  slope 

U.  The  main  stream  slope 
5.  The  watershed  slope. 

Evaluation  of  Geomorphologi cal  Factors 

1.  Drainage  Area 

Area  of  the  watershed  is  directly  measured  from  topographic 
maps  with  a  planimeter.   It  is  expressed  in  square  miles. 

2.  Drainage  Density 

Drainage  density  is  defined  as  the  total  length  of  streams  in 
the  watershed  divided  by  its  total  area,  that  is  the  length  of 
streams  per  unit  area  of  the  watershed.   It  can  be  expressed  as 
follows : 

Drainage  Density  (see  Figure  J) 

D  -  f*  (U) 

where  ZL  is  the  total  length  of  streams,  and 
A  is  drainage  area. 
The  length  of  stream  is  expressed  in  miles  and  is  obtained  from 
drainage  maps.   Since  the  area  is  expressed  in  square  miles,  the 
drainage  density  can  be  expressed  as  miles  per  square  mile. 
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Land  Slope 

Since  the  land  slope  is  changing  from  place  to  place  in  a 
watershed,  it  is  difficult  to  find  a  quantitative  value  to 
represent  the  land  slope  of  a  whole  watershed.  A  new  parameter 
introduced  here  to  replace  the  land  slope  is  the  mean  relief  of 
land.  The  mean  relief  is  defined  as  the  total  volume  of  land 
mass  above  the  outlet  of  a  watershed  divided  by  its  projected 
area.  This  can  be  evaluated  quantitatively  by  using  the  so-called 
hypsometric  curve  developed  by  Langbein  and  by  others  which  give 
a  dimensionless  relationship  between  the  horizontal  cross- sectional 
drainage  basin  area  and  the  elevation.  Fig.  k   shows  a  watershed 
and  its  horizontal  projection.  The  maximum  elevation  H'  may  be 
obtained  from  topographic  maps,  and  the  cross- sectional  area  A 
of  the  watershed  can  be  measured  by  planiraeter.   Similarly,  the 
projected  area  above  any  height  h  may  be  obtained  from  topographic 
maps.  The  dimensionless  plot  of  the  relative  area  a/A  against  the 
relative  height  h/H'  is  called  a  hypsometric  curve,  and  its  general 
aspect  is  as  shown  in  Fig.  5.  From  the  hypsometric  curve,  the 
total  volume  of  land  mass  and  the  mean  relief  can  be  calculated. 
Since  the  area  under  the  hypsometric  curve  can  be  easily  measured, 
this  can  be  expressed  as 

/  !  d(  -Sr7  )  «  or-  (12) 


A  "v   H 

o 


which  is  rewritten  as 


h 
~   /  adh  =  a'  (15) 


AH 
o 
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Hence  the  total  volume  of  land  mass  is 


h 
V  =  /  adh  =  AH'  q«  (ik) 

o 

and  the  mean  relief  is 

H  -  |  »  q'H'  (15) 

It  is  thus  seen  that  the  mean  relief  is  equal  to  the  product  of  the 
area  under  the  hypsometric  curve  and  the  maximum  height  over  the 
outlet . 

h.     Main  stream  slope 

The  slope  of  the  main  stream  can  be  obtained  from  the  topo- 
graphic map  at  several  points  along  the  length  of  the  stream.  Usually, 
the  upper  stream  reaches  are  steeper,  and  the  downstream  reaches  are 

flatter.  The  mean  slope  is  calculated  by  means  of  the  formula  in- 

( 18 ) 
troduced  by  Taylor  and  Schwarz: 

si 


1    ♦  A  *  4-  +  ...4- 

Sf     8g     8|       8| 


(16) 


whe 


re  n  is  numbers  of  equal  reaches,  S.,  S   S  ,  ...  and  S 
are  the  slopes  of  each  reach.   It  is  based  on  the  assumption 

that  the  quantity   ,  which  appears  in  the  Manning  formula, 

is  the  same  in  all  reaches. 

5.   Watershed  shape  factor 

The  watershed  shape  factor  is  the  ratio  of  the  main  stream 
length  to  the  diameter  of  an  imaginary  circular  watershed  of 
equal  area.   Figure  6  shows  two  watersheds  with  the  same  area  but 
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different  shape,  the  actual  irregularly  shaped  watershed  and  the 
imaginary  circular  one.  AO  is  measured  along  the  main  stream  up  to 
the  watershed  boundary  line,  and  A'O'  is  calculated  from  the  known  area 
of  the  actual  watershed. 


I  Area  of  actual  watershed  .      . 


(18) 


The  shape  factor  is  determined  as 

AO 


A'  0' 


Geomorphological  Factors  of  Small  Watersheds  in  Indiana 
The  following  Table  5  shows  quantitatively  the  geomorphological 
factors  of  16  small  watersheds  in  Indiana,  the  only  ones  for  which 

topographic  maps  were  available.  The  stream  length  and  stream  density 

(19) 
were  obtained  from  the  drainage  maps  of  Indiana. 

Multiple  Correlation  of  Peak  Discharge 
and  Geomorphological  Characteristics 

A  multiple  correlation  was  derived  between  the  25-year  flood 
and  the  geomorphological  variables  considered  above.   Such  a  correla- 
tion presumes  that  the  area  of  application  is  meteorologically  and 
geologically  homogeneous,  otherwise  the  state  should  be  divided  into 
different  zones  in  which  homogeneous  conditions  exist. 

Since  the  State  of  Indiana  is  relatively  flat,  the  orographic 
precipitation  is  not  a  factor  in  the  larger  storms.   Convectional 
precipitation  resulting  from  most  thunderstorms  usually  has  a  duration 
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Table  5 

Watershed  Characteristics  and  25- Year  Annual 

Instantaneous  Peak  Discharge  of 

16  Watersheds  in  Indiana 


25- Year  Annual 

Watershed 

Instantaneous 
Peak  Discharge 

Watershed  Characteristics 

number 

Area 

Mean 

Drainage 

Shape 

Main 

Q  cfs. 

A  8q.rai. 

Relief 

Density 

Factor 

Stream 

H.  ft. 

d  mi/ 
sq.rai. 

f 

Slope  . 
S  x  10" 

S-l 

11,800 

100 

216 

10.66 

2.63 

9.84  x  1 

S-2 

12,900 

88. 

8 

270 

7-58 

3.00 

20.88 

S-5 

2,950 

78. 

7 

101 

6.57 

1.75 

9.05 

S-6 

5,500 

62. 

9 

110 

8.00 

1.12 

21.10 

S-7 

1,650 

93. 

0 

79 

5.10 

1.47 

8.29 

1 

880 

115 

65.^ 

3-55 

l.Ui 

2.64 

6 

17 , 000 

179 

195.2 

7.88 

2.2 

13. 4o 

7 

10,500 

109 

86 

7.02 

1.94 

10.39 

9 

21,500 

156 

250 

9.76 

2.85 

10.63 

10 

19,000 

120 

257 

11.20 

2.18 

9.90 

11 

15,500 

257 

181.5 

13.95 

2.87 

2.95 

12 

8,700 

228 

99.8 

15.20 

1.95 

5.^3 

Ik 

15,500 

188 

195.8 

10.1*7 

1.55 

6.21 

15 

22,500 

129 

251 

8.70 

2.56 

10.16 

18 

1*,800 

125 

8^.7 

4.50 

1.91 

6.05 

19 

760 

102 

66.1 

3.16 

1.93 

7.93 

.-k 
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which  is  less  than  the  time  of  concentration  for  watersheds  of  more 
than  20  square  miles.  Thus  frontal  and  cyclonic  precipitation  are 
the  cause  of  the  large  storms  producing  peak  runoff.   Therefore,  it 
could  be  assumed  that  the  climatological  condition  is  homogeneous  over 
the  State  of  Indiana.  That  is,  the  probability  of  being  subjected  to 
a  storm  of  a  given  frequency  is  almost  equal  for  all  watersheds  in 
Indiana . 

By  studying  the  soil  map^    of  the  State  of  Indiana,  it  was 
found  that  the  permeability  varies  from  soil  to  soil.  Although  this 
would,  in  theory,  disprove  the  assumption  that  the  geological  conditions 

are  homogeneous  throughout  the  state,  one  additional  multiple  correla- 

(21) 

tion  including  the  maximum  intake  rate      of  the  soil  as  a  variable 

indicated  that  it  was  not  significant  compared  to  the  geomorphological 

variables. 

(22) 
Multiple  correlation      is  a  statistical  method  to  find  the 

relationship  between  one  dependent  variable  and  a  number  of  independent 
variables.   If  a  linear  relationship  exists,  the  method  of  fitting  is 
to  make  the  sum  of  the  squares  of  the  deviations  of  actual  observa- 
tions from  the  theoretical  linear  relation  a  minimum.  This  is  called 
the  method  of  least  squares. 

The  general  formula  for  multiple  correlation  is 

y  =  a  +  b^  +  b2x2  ...  bkxk  (19) 

where  y  is  the  dependent  variable,  x  ,  x  ,  x  ,  . ..  x   are  the 
variables,  and  a,  b  ,  b   b   ...  bfc   are  constants  obtained  by 
solving  the  following  simultaneous  equations: 
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b^x2)  +  b2S(X;Lx2)  +  b5S(Xlx3)      ...   +  bkS(xlXk)     =     S(X;Ly    ) 
b1S(x2x1)   +  b2S(x|)   +  b5S(x2x5)      ...   +  bkS(x2xk)  =     S(x2y    ) 


b1S(x5x1)   +  b2S(x5x2)  +  b5S(x^)      ...   +  bkS(X5xk)  =     S^y    ) 
blS(Vl>  +  b2S(xkX2>  +  *3*l\*J      ••'    +  V(xk}  "     S(V    } 


and 

a  =     y  -   b1x1  -   b2x2  -   b^   ...   -   b^  (20) 

where  S(xy)  is  a  symbol  which  means  Zxy  -  - 

so  that 


S(x2) 

=     Zx, 

Ux^2 

n 

S(x1x2) 

-        ZX,X2       - 

Zx,   Zx2 
n 

Sfv_*.   ) 

as        T.V     "K 

Zx,   Zx, 

1       k 

The  general  formula  for  estimate  of  variance  is 


S(y)2  -  b^^y)  -  b2S(X2y)  ...  -  bkS(xky  ) 


S2y,  xv   x2,  x^  ...  xk  = — "  ;  (k+l) 


where  k  is  the  number  of  independent  variables. 


(21) 


Tnis  serves  as  a  measure  of  the  degree  of  correlation.  Generally 
the  smaller  the  variance,  the  better  the  correlation. 

Tne  following  Table  6  shows  the  results  of  correlation  between  the 
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25-year  instantaneous  peak  discharge  and  the  geomorphological  factors 
of  small  watersheds  in  Indiana.  The  regression  formulas  in  exponential 
type  as  shown  in  Table  6  were  obtained  by  means  of  a  logarithmic  trans- 
formation and  the  formulas  (19)  and  (20).  The  standard  deviations 
were  calculated  by  means  of  Eq.  (21). 
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A  whole  set  of  multiple  correlations  was  listed  in  Ta"ble  6. 
By  studying  the  standard  deviations,  it  is  easy  to  find  the  significance 
level  of  geomorphological  factors  which  affect  the  discharge.   If  A  is 
assumed  to  be  the  first  in  significance  to  influence  the  discharge,  then 
the  H  will  be  the  second  one,  and  S,  D,  f  will  be  the  third,  fourxh, 
and  last,  respectively.  The  last  regression  formula  in  Table  6, 

Q  =  0.05565  A0'*715  H°'78M+  D0-82^  f°-4i6°  S0^01  (22) 

with  the  least  standard  deviation  is  of  course  the  best  expression  ob- 
tained from  multiple  correlation.  If  all  the  geomorphological  factors 
with  their  powers  together  are  combined  as  a  "basin  characteristic",  B  , 

B  =  a0'9715  E°'l3kk   D°*3254  f°-kl6°   sO-5901  (23) 

y 
then,  the  regression  formula  can  be  expressed  as  Q  =  cB,  where  c  is  a 

constant.  This  is  the  equation  of  a  straight  line  on  log- log  paper 

as  shown  in  Figure  7.   A  9%   confidence  interval  was  calculated  and 

plotted.  True  mean  will  fall  within  the  limits  with  a  probability 

of  y5  percent. 

To  test  the  assumption  of  geological  homogeneity  an  additional 

multiple  correlation  was  made  including  a  variable  representing  soil 

(21) 
characteristics.   The  maximum  intake  rate     of  the  soil  expressed 

in  inches  per  hour  was  selected  as  the  new  variable  in  addition  to  the 

five  geomorphological  variables  used  before.  The  six  variable  multiple 

correlations  which  would  take  care  of  the  nonhomogeneous  condition  of 

the  watershed  geology  did  not  show  much  improvement  in  the  degree  of 

correlation  to  the  peak  discharge,  having  a  standard  deviation  somewhat 

larger  than  that  for  Eq.  (22).  This  is  because  the  variation  in  the 
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intake  rate  is  small  for  the  watersheds  studied.  Until  additional 
data  on  runoff  from  small  watersheds  become  available,  Eq.  (22)  is  the 
best  relationship  that  can  be  derived  between  peak  flow  and  physio- 
graphic factors. 

Construction  of  Correlation  Chart 
for  Predicting  Future  Flood 
The  regression  formula,  Eq.  (22)  or  its  graphical  representation 
in  Figure  7  may  be  used  to  calculate  the  flood  discharge  if  the  geo- 
morphological  factors  are  known.  A  nomographic  representation  of  the 
regression  formula,  Eq.  (22),  was  prepared.   From  this  it  is  possible 
to  obtain  directly  the  flood  discharge  from  the  known  geomorpho logical 
factors. 

The  regression  formula  has  the  form 

Q  =  KA8  Bb  CC  Dd  Ee  (2k) 


Letting 


then 


Mx  =  AaBb 

M2  =  AaBbCC  =  M  CC 

Mx  =  ASBbCCDd  =  M0Dd 


<4     =  KAaBbCCDdEe  =  KM  Ee 


The  above  equations  may  be  plotted  with  B,  C,  D,  E  as  parameters 
and  b,  c,  d,  e  are  known  exponents  of  the  regression  formula  as  shown 
in  Figure  8.  The  diagrams  of  Fig.  8  may  be  combined;  the  "A  "  scale 
is'  changed  from  "A"  scale  so  that  the  area  A  can  be  introduced  directly 
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FIG      8-  GRAPHICAL      SOLUTION      OF      EQ       <?4) 
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FIG.      9-        GRAPHICAL       SOLUTION        OF      THE 
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into  the  chart  along  with  the  other  factors  necessary  to  find  the  flood 
discharge,  as  shown  in  Figure  9. 

The  above  procedure  was  used  to  prepare  the  nomograms  of  Figures 
10  and  11.  These  may  be  used  instead  of  Eq.  22  or  of  Figure  7  to  esti- 
mate the  2^-year  peak  discharge  from  the  geomorphological  characteris- 
tics.  Fig.  10  was  prepared  for  watersheds  less  than  100  square  miles 
and  in  Fig.  11,  the  correlation  is  extrapolated  to  cover  watersheds  up 
to  300  square  miles.  Two  examples  are  presented  to  show  the  use  of 
the  charts: 

(a)  Watershed  No.  S-6 

Watershed  area 
Mean  relief 
Drainage  density 
Shape  factor 
Main  stream  slope 

Following  the  dotted  line  in  Fig.  10,  the  flood  discharge  is  read  from 
the  chart  directly  as  3,750  cfs.,  while  the  flood  predicted  by  the 
frequency  study  is  3,300  cfs. 

(b)  Watershed  No.  S-l 

Watershed  area  (A)  100  sq.  mi. 

Mean  relief  (H)  216  ft. 

Drainage  density  (D)  10.66  mi./sq.  mi. 

Shape  factor  (f)  _k                 2.63     _k 

Main  stream  slope  (S)  x  10~        9.8U  x  10~ 

Following  the  dotted  line  in  Fig.  11,  the  flood  discharge  reading  is 
12,1*50  cfs.,  while  the  flood  predicted  by  the  frequency  study  is 
11,800  cfs. 


(A) 

62.9  sq.    mi. 

U) 

110  ft. 

(D) 

8  mi./sq.   mi. 

(S)   x   10   * 

1.12             _k 
21.10  X  10 
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FIG.     II  -         WORKING       CHART       FOR       REGRESSION       FORMULA 
EQUATION       (22) 
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A  Simple  Approximate  Estimation  of  Peak  Discharge 
for  Small  Watersheds  in  Indiana 
As  shown  in  the  previous  paragraphs,  equation  (22)  is  the  best 
one  among  the  ten  formulas  listed  in  Table  6  for  predicting  the  peak 
discharge  for  small  watersheds  in  Indiana  since  it  has  the  least  stan- 
dard deviation.  Although  a  graphical  representation  of  the  correlation 
formula,  Eq.  22, is  given,   its  use  is  time-consuming  because  of  the 
tedious  work  required  for  the  determination  of  the  five  geomorphological 
factors  which  enter  into  the  formula  or  in  the  correlation  charts  as 
independent  variables.   In  particular,  the  determination  of  the  drainage 
density  D  and  of  the  mean  relief  H  is  time-consuming.   For  the  multiple 
correlations  listed  in  Table  6,  it  appears  that  the  formula  containing 
A,  H  and  S  has  a  standard  deviation  of  0.211  which  is  close  to  the 
standard  deviation  of  0.190  for  Eq.  22  and  may  therefore  be  used  as  a 
good  approximation  in  practical  design.   This  formula  eliminates  the 
need  of  calculating  the  drainage  density  D,and  the  watershed  shape 
factor  f  which  is  the  least  significant  variable  is  omitted.  The  formula 

Q  =0.0022  A1'4620  h1-5°55s0.6938  (25) 

still  requires  the  determination  of  the  mean  relief  H  of  the  watershed. 
However,  there  is  an  approximate  way  of  determining  H  by  estimating  a' , 
the  area  under  the  hypsometric  curve.  The  average  hypsometric  curves 
have  been  plotted  for  small  watersheds  in  Indiana  in  Fig.  12,  from 
which  it  appears  that  the  a' -values  vary  from  O.k   to  0.8.  From  the 
pattern  of  the  hypsometric  curves  it  is  possible  to  estimate  a'    by 
calculating  only  one  value  of  a/A  for  a  certain  elevation  ratio, 
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say  h/H1  =  0.5-  Thus  two  area  measurements  are  needed:   the  total  area 
of  the  watershed  A,  and  the  area  at  an  elevation  half-way  between  the 
maximum  elevation  and  the  mouth.  The  following  Table  7  gives  the  a' 
values  for  corresponding  a/A: 

Table  7 
The  q'  Values  Corresponding  to  a/A 
for  b/H'  =0.5 


a/A  a' 

0.3  O.lj-0 

o.u  0.U5 

0.5  0.50 

0.6  0.55 

0.7  O.60 

0.8  0.70 

0.9  0.80 


For  a  rough  estimate  of  a' ,  it  may  be  assumed  that  it  varies  between 
O.k-   to  0.5  for  U-shaped  valleys  and  from  0.6  to  0.7  for  V-shaped  valxeys 
in  Indiana.  With  the  value  of  a'  determined  from  Table  7,  the  mean 
relief  H  can  be  easily  calculated  from  Eq.  15,  and  the  peak  discharge 
can  be  obtained  by  Eq.  25. 

Thus,  the  peak  discharge  can  be  calculated  by  three  methods.  The 
first  is  the  simplest;  a'  is  obtained  from  Table  7,  from  which  by  Eq. 
15  and  Eq.  25  a  first  approximation  of  the  peak  discharge  is  obtained. 
In  the  second  method  H  is  determined  by  plotting  the  hypsometric  curve 
and  Eq.  25  is  used  as  a  second  approximation  of  the  peak  discharge. 
The  third  method  requires  the  evaluation  of  the  five  geomorphological 
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factors,  after  which  the  peak  discharge  is  calculated  by  Eq.  22  or  by 
means  of  Figures  10  and  11. 

An  example  is  given  below  and  solved  by  the  three  methods. 

Watershed  No.  S-l 

First  Method: 

Geomorphological  factors 

A  =  100  square  miles 
s  «=  $.&k 
H'  »  U00  feet 
at  h/H'  =  0.5,  a/A  =  0.575,  from  Table  7,  a'  =  0.55 
H  =  U00  x  0.55  =  212  feet 

by  Eq.  25,  the  25-year  peak  discharge  is 
Q  «  9,280  cfs. 

Second  Method: 

Geomorphological  factors 

A  =  100  square  miles 

S  =  Q.&k 

H  =  216  feet 

by  Eq.  25,  the  25- year  peak  discharge  is 

Q  *  9,U8o  cfs. 

Third  Method: 

The  25-year  peak  discharge  as  found  by  using  the  regres- 
sion formula,  Eq.  22  or  Figures  10  and  11  as  shown  in  the 
previous  example  relating  to  those  charts,  is 

Q  =  12,^50  cfs. 

while  the  25- year  peak  discharge  estimated  from  the  frequency 
study  was  11,800  cfs. 

Comparing  the  results  obtained  by  the  three  methods  to  the  peak 
discharge  obtained  from  the  frequency  study,  it  appears  that  the  percentages 


uo 


of  error  are  respectively  -21. hi,   -l^.fi,   and  +5.5*.   It  should  be 
remembered,  however,  that  these  percentages  will  vary  from  one  water- 
shed to  the  next  as  methods  1  and  2  are  calculated  by  means  of  one 
regression  formula  whereas  method  3  is  obtained  by  a  different  regres- 
sion formula. 

Relationship  between  25- year  Peak  Discharge 
and  the  Peak  Discharge  for  Other  Frequencies 
In  the  preceding  paragraphs,  the  peak  discharges  from  small 
watersheds  were  obtained  for  a  recurrence  interval  of  25-years  which 
was  based  on  the  average  life  of  small  highway  drainage  structures. 
However,  it  may  be  desirable  to  estimate  the  peak  discharge  for  other 
return  periods  so  that  the  design  engineer  may  have  a  greater  freedom 
of  choice.  Hence  the  relationship  between  the  peak  discharge  for  any 
frequency  and  the  25-year  peak  discharge  was  derived  for  small  water- 
sheds in  Indiana.  The  theoretical  relationship  is  based  on  Gumbel's 
extreme  value  theory.  Figure  13  shows  two  theoretical  straight  lines 
for  any  two  watersheds.  The  differences  between  the  25-year  peak 
discharge  and  the  n-year  peak  discharge  for  the  two  watersheds  are  ob- 
tained from  Eq.  (10),  and  are: 

">1  *  -\* 

AQ2  =  -~  Ay  (26) 

where  ,  are  the  slopes  of  the  straight  lines,  and  the 

a±  '     Op  *  ' 

increment  Ay  of  the  reduced  variate  corresponds  to  the  selected 
increment  of  frequency.  For  a  fixed  frequency  n  the  increment  Ay  is 
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a  constant.  A  general  form  thus  can  be  written  for  all  the  watersheds 


as 


AQ   =  ^Ay  (27) 


or 


a  Qn  -  a  Q25  =  Ay 


thus 


a  <=tn  -  a  Q2«-  +  Ay 


«n  -  V1  +  H^~  »       <28) 
Taking  logarithms  on  both  sides 

log  Q^   =  log  Q25  +  log  (1  +  a^     )  (29) 

Eq.  29  is  a  linear  on  log- log  paper,  if  the  last  term  is  con- 
Ay 


stant.  An  examination  of  the  last  term  showB  firstly  that  the 


a  ^25 


is  small  compared  to  "1"  and  secondly  that  the  value  of  a  Q?e.  for 
most  of  the  small  watersheds  in  Indiana  varies  from  h   to  6  and  an 
average  value  of  5  can  be  used  for  a  Q?c.  for  small  watersheds  in 
Indiana.  Then,  Equation  29  becomes 

log  Qn  =  log  Q25  +  log  (1  +  &)  (50) 

Values  of  Ay  for  several  design  frequencies  are  given  in  Table  8. 
Values  of  Ay  for  other  frequencies  may  be  obtained  from  probability 
paper. 
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Table  8 
Values  of  dy  for  Different  Frequencies 


n 

Ay 

10 

-0.93 

50 

0.70 

75 

1.00 

100 

1.1*0 

Fig.  Ik   is  a  plot  of  Eq.  50  on  logarithmic  paper,  giving  the  rela- 
tionship between  the  25-year  peak  discharge  and  the  peak  flow  for  fre- 
quencies of  10,  50,  75  and  100  years.  Such  peak  discharges  may  be  de- 
termined from  Figure  lb   or  from  Eq.  (50)  if  the.  25-year  discharge 
has  previously  been  determined  from  Figures  11  and  12  or  from  Eq.  (22). 
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DISCUSSION  AND  CONCLUSIONS 

1.  The  paper  considers  watersheds  between  20  and  200  square  miles 
in  area  in  the  state  of  Indiana.  All  available  past  observations 
of  annual  peak  discharge  were  plotted  on  a  probability  paper  using 
Gumbel's  extreme  value  theory.  As  there  is  a  linear  relation  be- 
tween the  observations  and  the  reduced  largest  value  y,  the  best 
fit  straight  line  was  then  obtained.  Future  floods  with  different 
frequencies  were  obtained  by  extending  the  straight  line.  Table  3 
gives  the  predicted  flood  of  25,  50,  75  and  100  years  of  frequency 
for  52  gaged  watersheds  in  Indiana. 

2.  The  quantitative  study  of  geomorpho logical  factors  of  small 
watersheds  in  Indiana  combined  with  the  use  of  multiple  correla- 
tion techniques  gives  an  indirect  determination  of  peak  discharge. 
This  is  based  on  the  assumption  that  the  climatological  and  geo- 
logical conditions  are  homogeneous  throughout  the  state.  Hence, 
the  geomorphological  characteristics  are  the  dominant  factors 
which  affect  the  peak  discharge  from  small  watersheds.  The  geo- 
morphological factors  considered  significant  are:  the  watershed 
area,  the  drainage  density,  the  mean  relief  of  watershed,  the 
main  stream  slope,  and  the  shape  factor  of  the  watershed. 
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3.      Correlation  charts  (Figs.  10  and  U)  were  prepared  to  obtain 
the  25-year  peak  discharge  directly  from  the  five  watershed  charac- 
teristics, for  areas  up  to  500  square  miles.   As  shown  in  the 
previous  examples,  the  design  engineers  may  use  these  design  charts 
to  estimate  very  rapidly  the  25-year  peak  discharge  with  good 
accuracy.  The  peak  discharge  for  other  frequencies  may  be  obtained 
from  Fig.  Ik. 

k.  A  simple  formula  which  contains  only  three  geomorphological 

factors,  A,  H,  and  S  was  introduced  as  a  first  approximation.  H  is 
determined  by  the  a' -value  which  is  estimated  from  the  average 
hypsometric  curves  for  small  watersheds  in  Indiana. 
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PART  II 

DESIGN  HYDROGRAPHS 
FOR 
SMALL  WATERSHEDS  IN  INDIANA 
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INTRODUCTION 
In  the  design  of  many  hydraulic  structures,  the  engineer  is  con- 
cerned not  only  with  the  maximum  discharge  but  also  with  the  total 
volume  of  runoff  and  its  distribution  with  respect  to  time.  The  rout- 
ing of  floods  through  a  reservoir  to  determine  spillway  requirements 
may  be  cited  as  an  example  in  which  the  total  hydrograph  is  necessary. 
The  runoff  hydrograph  is  defined  as  a  graph  of  discharge  versus 
time.  Its  magnitude  and  shape  depend  upon  the  characteristics  of  both 
the  watershed  and  the  runoff- producing  storm.  The  ascending  limb  of 
the  hydrograph  is  influenced  primarily  by  the  storm  pattern  and  water- 
shed characteristics,  while  the  recession  limb,  after  the  point  of 
inflection,  represents  the  depletion  of  water  from  storage  and  its 
shape  is  influenced  by  the  characteristics  of  the  watershed  only. 

The  purpose  of  this  study  is  to  determine  the  general  shape  of 
hydrographs  for  small  watersheds  by  means  of  a  mathematical  expression 
containing  certain  parameters  which  can  be  correlated  with  identifiable 
and  readily  obtainable  physical  characteristics  of  the  watershed.  One 
may  then  develop  a  design  hydrograph  for  ungaged  areas  from  the  de- 
sired design  precipitation  and  certain  watershed  characteristics  which 
can  be  determined  from  a  topographic  map  of  the  basin. 

As  used  in  this  study,  the  term  "small  watershed"  means  those  of 
less  than  100  square  miles  in  area. 

The  earliest  concept  of  the  unit  hydrograph  was  presented  by 
L.  K.  Sherman  ^2^  in  1952.   It  is  defined  as  a  typical  hydrograph, 
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representing  one  inch  of  direct  runoff  generated  uniformly  over  the 
area  at  a  uniform  rate  during  a  given  period  of  time.  Since  the 
physical  characteristics  of  the  watershed  are  constant,  one  might 
expect  considerable  similarity  in  the  shape  of  hydrographs  from 

storms  of  similar  rainfall  characteristics. 

(2k) 
In  1939,  Snyder    '  presented  procedures  for  the  development 

of  synthetic  unit  hydrographs  for  ungaged  watersheds.  Formulas  were 
derived  for  three  elements  of  the  hydrograph,  time  to  peak,  peak 
discharge,  and  the  time  base  by  relating  them  to  watershed  charac- 
teristics. These  three  items,  plus  the  fact  that  the  volume  must 

equal  one  inch,  permit  the  sketching  of  the  complete  unit  hydrograph. 

(25) 
Edson,  v   '  in  1951,  derived  a  theoretical  expression  for  the 

unit  hydrograph.  He  considered  that  the  physical  characteristics  of 

a  watershed  exert  two  simultaneous  and  distinct  influences  upon  the 

resultant  unit  hydrograph: 

a.  whereby  the  runoff  is  brought  to  the  valley, 
so  that  Q  «  t  , 

b.  whereby  the  runoff  through  the  mouth  is  such  that 
Q  «  •-**. 

Combining  these  two  influences,  so  that  Q  «  t  e    and  the  total 
runoff  is  "1"  inch,  he  obtained  the  following  expression  for  the 
unit  hydrograph: 

c  A  y(yt)X  e_yt 


T(x+1) 


(3D 


In  Eq.  31,  Q  is  the  discharge  in  cfs  at  time  t,  t  is  the  time  in 
days  from  the  beginning  of  runoff,  A  is  the  drainage  area  in  square 
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miles  and  x,y  are  parameters  which  are  expressible  in  terms  of  peak 
discharge  and  the  time  to  peak,  c  is  a  coefficient  for  the  adjustment 
of  the  units  and  the  denominator  contains  the  Gamma  function  of  (x  +  1). 

Recently,  an  expression  for  the  unit  hydrograph  was  derived 
making  use  of  the  concept  of  unit  impulse  function  and  of  the  convolution 
integral.  The  unit  impulse  function  may  be  considered  as  an  instantaneous 
inflow  which  consists  of  a  burst  of  rainfall  of  constant  intensity  acting 
only  during  an  infinitesimally  small  period  of  time  so  that  the  product 
of  duration  and  intensity  is  one  inch  of  rainfall.  The  instantaneous 
unit  hydrograph  u(t)  may  be  regarded  as  the  impulsive  response  of  the 
watershed  to  an  instantaneous  inflow  as  shown  in  Figure  15. 

The  convolution  integral  or  Duhamel's  integral  is  used  to  de- 
termine the  response  Q(t)  of  the  system  to  any  arbitrary  excitation 

(26) 
or  rainfall  excess  i(t)  if  the  unit  impulsive  response  u(t)  is  known: 

t 
Q(t)  -  /   i(r)  u(t  -  T)dT  (52) 

o 

where  t  represents  the  time  and  is  the  dummy  variable  of  integration. 

For  a  one- inch  rainfall  excess  of  constant  intensity  and  of  duration  D, 

the  excitation  is  i(t)  =  =  and  Eq.  (52)  becomes: 

t 


U(D,t)  =  -J  /  u(t  -  r)dT  (55) 


o 

which  is  the  equation  of  the  unit  hydrograph  of  duration  D.  If  an 
instantaneous  inflow  i(o)  is  applied,  t  is  zero  and  the  time  t  ap- 
proaches At,  then  Eq.  52  yields  the  instantaneous  hydrograph  q(t): 


q(t)  =  i(o)  At  u(t)  -  Pe  u(t)  (5M 


%r 


t(hr.] 


INSTANTANEOUS    INFLOW 


INSTANTANEOUS   UNIT   HYDROGRAPH 


FIG.  15-  THE  INSTANTANEOUS  INFLOW  AND 
THE  INSTANTANEOUS  UNIT  HYD- 
ROGRAPH 
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where  P  is  called  the  instantaneous  "effective  rainfall"  to  cause 

e 

runoff.  The  expression  for  the  instantaneous  hydrograph  has  a  simple 
form  so  that  it  may  be  used  in  practical  engineering  design.  For  a 
given  effective  rainfall  P  the  instantaneous  hydrograph  gives  a  higher 
peak  than  any  hydrograph  of  a  finite  duration.  Thus  the  instantaneous 
hydrograph  yields  the  upper  limit  of  the  peak  discharge.  For  small 
watersheds  where  the  time  of  concentration  is  very  short  the  instanta- 
neous hydrograph  gives  a  good  estimate  of  the  peak  discharge. 

Expressions  for  the  instantaneous  hydrograph  were  derived  by 

(  ?7  }  (  2ft  ^ 

Nash     and  Dooge     with  the  assumptions  that  an  instantaneous 
inflow  is  applied  to  n  equal  linear  reservoirs  with  the  same  storage 
coefficient  K.  They  considered  that: 

a.  for  a  linear  reservoir,  the  storage  S  is  related  to  the 
outflow  q  by 

S  =  Kq  (35) 

b.  for  an  instantaneous  inflow,  the  outflow  from  a  linear 

reservoir  is  given  by 

t 

q  =  ■£  e  (36) 

where  V  «  AP  is  the  total  volume  of  runoff  from  the 

e  . 
basin  of  area  A. 

Then  for  a  number  of  n  equal  linear  reservoirs  in 
series  with  the  same  storage  coefficient  K,  the  outflow 
is  given  by: 

•     q  .  V   (t/K)"'1  e^  (57) 

K        (n  -  1)! 

Eq.  37  may  also  be  written,  by  expressing  the  factorial  in  terms  of 
the  Gamma  function,  as 
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V  (t/K)  D"1  e"1^ 
q=  K    '  '  r(n) (58) 

The  instantaneous  hydrograph  is  thus  given  by  an  expression  containing 

a  single  term  with  two  parameters  n  and  K  which  determine  the  shape  of 

the  hydrograph. 

Gray/  *'   in  i960,  has  shown  that  the  formula  derived  by  Edson 

can  be  transformed  to  the  form  of  Eq.  38  by  a  change  of  variables, 

although  the  assumptions  are  different  from  those  of  Nash  ana  Dooge. 

He  also  noted  that  Eq.  38  is  Pearson's  type  III  frequency  distribution 

which  is  known  as  one  of  the  most  common  frequency  curves  and  which  has 

been  used  numerous  times  in  the  analysis  of  hydrologic  data.  Gray,  *' 

in  the  same  year,  developed  the  synthetic  unit  hydrographs  for  small 

watersheds  by  using  the  formula  as  derived  by  Edson,  Nash  and  Dooge. 

The  measurable  characteristics  were  correlated  with  two  hydrograph 

parameters,  the  Gamma  function  argument  n  and  the  storage  coefficient  K. 
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HYDROGRAPH  STUDY 

Since  the  recession  curve  of  the  hydrograph  represents  the  with- 
drawal of  water  from  storage ,  there  must  be  some  definite  relationships 
between  the  recession  curve  and  the  watershed  characteristics.  Measurable 
watershed  characteristics  can  be  evaluated  by  a  geomorphological  study 
of  the  watershed,  and  the  relations  between  the  recession  curve  and  the 
watershed  characteristics  can  be  obtained  by  means  of  multiple  correlation 
techniques,  as  outlined  in  the  following. 

Derivation  of  Working  Formulas 
From  the  formula  for  the  instantaneous  hydrograph  Eq.  38,  if  a 
total  delay  time  K  (storage  coefficient)  is  assumed,  the  effect  on  the 
outflow  hydrograph  of  varying  the  value  of  n  gives  the  interesting  result 
that  is 

t  =  (n  -  1)K  (39) 

where  t   =  time  to  peak  discharge. 

This  can  be  proven  by  taking  the  first  derivative  of  Eq.  38  with 
respect  to  time  and  setting  ~  equal  to  zero. 

Then,  letting  X  =  -g     and  substituting  t  ■  (n-l)K  and  AR  -  V 
where  A  is  the  watershed  area  in  acres  and  R  is  the  depth  of  runoff,  in 
inches,  into  Eq.  (58),  we  have 

qt     /   ,  »  Yn-1  -X 

_LE  =  (n-1)  X"   e  (^} 

AR         r(n) 
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where  X  =  —  =  —  (n-1) 

P 

Eq.  ^0  is  the  general  formula  for  the  instantaneous  unit  hydro- 
graphs  used  in  this  study.   The  function  of  n  and  t  in  the  right 
member  of  Eq.  kQ   will  be  denoted  by  f(n,t)  and  Eq.4o  becomes 

qt 

j^   -  f(n,t)  (40-a) 

if  t  =  t  ,  then  q  »  q 

^  -  f(n,tp)  (*0-b) 

where  a  is  the  maximum,  or  peak  discharge  of  the  instantaneous  hydrograph. 

This  is  the  equation  to  determine  the  value  of  n  if  the  other 
factors  are  known,  where  n  is  the  argument  of  the  Gamma  function.   It  is 
the  major  factor  in  determining  the  shape  of  the  hydrograph.   It  is  also 
interesting  to  mention  here  that  the  gamma  distribution  itself  has  a 
shape  similar  to  that  of  the  hydrograph. 

A  study  of  Eq.  1*0  discloses  the  following  points  of  special 
interest: 

(1)  Instantaneous  Unit  Hydrograph. 
By  integrating  Eq.  kOt   the  total  volume  under  the  hydrograph  is 

r  a*       t     AR  (n-D  y?'1   e"X  dt  ,.* 

/q  dt  -  /  ^ ^ (U) 

P 
AR   ,   n-1   n-1   -X  ,.  /..,   v 

'run  J   t  x      e    dt  {kl-a) 

p 


AR     ,  .Jl-1   -X 

rjn)     !  ^      e   d 


(n-l)t 

t 
P 


(*l-b) 
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=  -^5,  r(n)  (Vi-d) 

=  AR  (kl-e) 

Since  AR  is  the  volume  of  runoff,  when  R  equals  1  inch  of  runoff 

over  the  total  area  then  Eqs.  ho   and  kl   refer  to  the  instantaneous  unit 

hydrograph  Eq.  (Uo-a)  becomes  u=  —  f(n,t)  .  (4l-f) 

P 

(2)  Peak  Discharge  Computation. 
Since  (Eq.  ko-b) 


if 


then 


Sn  =     -  f<n>   tp) 

p 

f(n,tp)  =  0.75 

64o  x  0.75  AR 

464  AR 

S»  "        t 

P 

t 
P 

e*2) 


(*2-a) 

Eq.  (42-a)  is  the  formula  for  computing  the  peak  discharge  used  by  the 
United  States  Soil  Conservation  Service,  A  being  the  watershed  area 
in  square  miles. 

(3)  Rational  formula. 

Since  from  Eq.  (40-b) 

'n.i 
P 

by  substituting  t  -  c.t   where  t  is  the  time  of  concentration 
and  c.  is  a  constant, 


P 
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then 


'1  c 


This  has  the  rational  form. 


\*    ^rf^V  (U5) 


\   =  CIA  (43-a) 


where 


f(n,  t  ) 

c- — — £-  (U>b) 


ci 


and 

I  =   5_  =  in/hr  (^5-c) 

c 

The  above  studies  identify   Eq.  (ko)      to  the  characteristics 
of  the  hydrograph  and  of  the  rational  formula.   The  empirical  formula 
used  by  the  United  States  Soil  Conservation  Service  is  only  a  special 
case  where  f(n,t  )  =  0.75. 

Development  of  Dimensionless  Instantaneous  Hydrograph 

The  dimensionless  instantaneous  hydrograph  is  defined  as  a  graph 

of   —  and  — —  as  ordinate  and  abscissa,  respectively.  The  mathe- 

td        p 
matical  expression  can  be  derived  from  Eq.  (4o)  as  follows: 
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qtI 
"IF 


12£ 


Yn-1  -X 
X    e 


(-D  I- 


when 


V  4 "  "» 


AR 


(n-1)  Yn-1 

m     n-1 


-X 


(kk) 


therefore, 


SL 
\ 


n-1 


(n-1)  f 
E. 


(n-1) 


t/t 

-I    '   T 


(n-1) 


n-1 


or 


TO      p     L 


(n-1) 


(f  -1) 


(*5) 


(^5-a) 


The  above  equation  gives  the  relation  between 


Tn      p 
any  given  value  of  n,  and  of  course,  it  can  be  plotted  as  a  dimension- 


less  hydrograph  for  any  given  value  of  n.  However,  the  computation  is 

rather  complicated,  and  hence  a  more  convenient  way  to  determine  the 

dimensionless  instantaneous  hydrograph  is  provided  by  the  data  in 

qt 
Appendix  B-l  which  gives  the  value  of   ttk*  for  different  values  of 

qt  qt 

n  and  t/t   .  The  value  of    7-5^       will  be  equal  to        -.J?     for  t  ■  t   , 

'    p  AR  *  .  AR  p' 

qtP         w 

The  ratio  of  -rs*"  and 


that  is  r-  -  1. 


AR 


AR 


gives  the  value  of 


Then  the  dimensionless  instantaneous  graph  can  be  plotted  as 


A. 

shown  in  Figure  16. 
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FIG.    16  - 


DIMENSIONLESS        INSTANTAN  EOU  S 


HYDROGRAPH 
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Recession  Curve  Study 

This  study  forms  the  basis  for  the  concept  that  n  can  be  corre- 
lated with  recession  curves,  that  the  value  of  n  can  therefore  be 
determined  from  the  recession  curve,  and  conversely,  that  the  reces- 
sion curve  can  be  determined  from  the  value  of  n. 

The  recession  curve  of  the  hydrograph  starts  from  the  point  at 
which  surface  inflow  to  the  channel  system  ceases  and  the  flow  is 
derived  wholly  from  the  withdrawal  of  water  from  storage.  If  a 
linear  storage  relationship  is  assumed,  the  recession  curve  will  plot 
as  a  straight  line  on  semilogarithmic  paper.  Linear  storage  means 
a  linear  relationship  between  the  storage  and  the  outflow.  The  storage 
coefficient  K.  can  be  determined  by  the  formula 

Kx  -  &—  (k6) 

2.5  log -^2 

where     q   ■  flow  at  any  time 

q.   =  flow  at  At  time  units  after  t 
1  o 

At  =  time  interval  form  t  to  t, 

o     1 

K1  =  storage  coefficient  for  a     to  q.,  obtained  from 
analysis  of  the  recession  curve. 

The  recession  curve  of  the  dimensionless  hydrograph  can  be  expressed  as, 

K..      At/t 

^  . L* (l*6-a) 

Kl 

Therefore,  —  can  be  determined  by  plotting  the  recession  part 

p 

of  the  dimensionless  hydrograph  on  semilogarithmic  paper. 
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By  introducing  the  above  recession  formula  into  the  mathematical 

formula  for  the  instantaneous  hydrograph  (Eq.  ho)   a  relation  between 

K  /t  and  n  can  be  derived  as  follows: 
1'  p 

from  Eq.  (U^-a) 

n-1  r 


—  -  (  r-  ) 

Til        p 


IV- 1)1 


(r-1' 


then 


t   n-1  r  .  i  (  -2  .  1) 

(o       re-(n.i)rt 

q     * 
^o     p 


t>> 


n-1 


(n-1) 


p 


<*7) 


taking  the  logarithm  of  both  sides, 


Ln 


Ln 


1       p         p 


i-  (n-1)  Ln^*  £     (n-1) 
ql  1     p 


-(n-1) 


(*7-a) 


(Vf-b) 


Introducing  Eq.  (k6-a) 


Ln 


-2  =*  (n-1)  Ln  t£  +  (n-1)  Ln 
q„  t. 


(i»7-c) 


By  simplifying, 
(n-1) 


1  _  (n_i)  _i   Ln  -2  *  (n-1)  Ln   ^ 
ql  1 


C+7-d) 
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1  -   (n-l)     <i 

L                                 p   . 

Kl 

1  -   (n-l)    xi  - 

P 

Kl 

t     * 
P 

IPg 
(n-l) 

1 
n-l 

-2  =  (n-l)  Log  -2 
ql                           Xl 

t 
log  +2 
*1 

t 

T                    ° 

Logr 

X 

T        % 

Log- 

<*7-e) 


(*7-f) 


(*7-g) 


The  final  expression  shows  the  relationship  between  n  and  K  /t  . 
Since  the  second  term  of  the  right  member  of  the  equation  is  not  a 
constant,  the  relationship  cannot  be  linear.  A  semi- logarithmic  plot 
of  Eq.  (U6-a)  is  shown  in  Fig.  17.   In  addition,  in  the  top  left  corner 
of  Fig.  17,  the  recession  curves  of  the  dimensionless  hydrographs  of 
Fig.l6  have  been  superimposed,  using  n  as  parameter.  Fig.  17  may  be 
used  as  a  tool  to  determine  the  Gamma  function  argument  n  which  gives 
the  shape  of  hydrograph  corresponding  to  the  two  hydrograph  parameters 
K.  and  t  for  any  given  watershed. 

Geomorphologlcal  Study  for  Small  Watersheds 
The  purpose  of  the  geomorphological  study  of  small  watersheds 
is  to  evaluate  various  factors  which  can  be  determined  directly  for 
the  topographic  map  in  order  to  make  correlations  with  the  hydrograph 
parameters  and  thence  to  the  design  flood  hydrograph  for  ungaged  water- 
sheds . 
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FIG.     17  -  RELATIONSHIP        BETWEEN        GAMMA 

FUNCTION        ARGUMENT       n     ft     K,/fD 
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Five  factors  were  evaluated  for  this  study,  as  follows: 

(a)  Drainage  area  (A): 

The  drainage  area  of  the  watershed  is  measured  directly  from  the 
topographic  maps  with  a  planimeter.   It  is  expressed  in  square 
miles. 

(b)  Length  of  main  stream  (L): 

The  length  of  main  stream  is  measured  in  miles  with  a  map  measurer 
along  the  main  stream  as  shown  on  the  topographic  maps. 

(c)  Mean  slope  of  main  stream  (S): 

The  slope  of  the  main  stream  can  be  obtained  from  a  study  of  the 
topographic  maps.   It  is  not  the  same  throughout  the  whole  stream. 
Usually,  the  slope  of  the  upper  reaches  of  the  stream  is  steeper, 
and  that  of  the  downstream  reaches  is  flatter.   It  would  not  be 
sufficient  Just  to  find  the  slope  of  the  straight  line  which 

connects  the  upper  and  lower  extremities  of  the  stream  profile. 

(l8) 
A  method  introduced  by  Taylor  and  Schwarz  v   ;  for  determining  the 

mean  slope  for  a  stream  is  used  in  this  study.   The  formula  for 

determining  the  mean  slope  is  as  shown  in  Eq.  16. 

(d)  Watershed  shape  factor  (f ): 

Differences  in  the  shape  of  the  watershed  will  result  in  different 
times  of  concentration  and  hence  will  affect  the  shape  of  the 
hydrograph.   The  watershed  shape  factor  is  designed  to  compare 
the  irregular  shape  of  the  watershed  to  an  idealized  shape,  as 
shown  in  Fig.  18.   It  is  defined  as  the  ratio  of  the  length  of 
the  perimeter  of  the  actual  watershed  to  the  length  of  perimeter 
of  a  circle  which  has  the  same  area  as  that  of  the  watershed  studied. 
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A  =  AREA 

P=  PERIMETER 


ACTUAL   WATERSHED    SHAPE 


A  =  AREA 


P-2Wr-  2Tr*/A. 
1T 


IDEALIZED    WATERSHED   SHAPE 


FIG.    18-  THE      ACTUAL     WATERSHED        SHAPE 


AND      THE       IDEALIZED       WATERSHED 


SHAPE 
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Let  A  be  watershed  area,  P  its  perimeter  and  P'  the  perimeter  of 
a  circle  of  equal  area,  then,  the  watershed  shape  factor  is 
defined  by  the  relation 

f  -  |r  <*8) 

where      f  ■  the  watershed  shape  factor 

P  «  the  length  of  perimeter  of  actual  watershed 

P1  «  the  perimeter  of  a  circle  having  the  same  area 
as  the  actual  watershed. 

The  watershed  shape  factor  is  always  larger  than  one,  since 
the  perimeter  of  a  circle  is  considered  the  smallest  as  compared 
to  the  perimeter  of  any  irregular  shape  with  the  same  area  as  the 
circle.  Therefore,  it  can  be  said  that  the  larger  of  the  water- 
shed shape  factor,  the  more  irregular  the  shape  of  the  watershed, 
(e)  Valley  shape  coefficient  (v): 

It  may  readily  be  conceived  that  the  shape  of  the  valley  has 
an  influence  on  the  hydrograph  as  a  result  of  the  effect  of  valley 
storage . 

The  quantitative  determination  of  a  coefficient  to  express 
the  influence  of  valley  shape  is  admittedly  difficult.  This  study 

is  based  on  the  so-called  hypsometric  analysis  which  was  developed 

(15) 
by  Langbein  and  others.  v   '   The  valley  shape  coefficient  is 

determined  as  follows: 

1.   Plotting  the  hypsometric  curve  for  the  given  watershed, 
shown  in  Fig.  5, 

where      A  «  total  area  of  the  watershed 

H'  «  total  height  of  the  watershed 

a  ■  area  enclosed  in  a  given  contour  line 

h  «  the  height  corresponding  to  "a" 
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The  values  of  A,  H',  a  and  h  can  be  measured  directly  from  the 
topographic  map. 

2.  Measuring  the  area  under  the  hypsometric  curve  for  different 

—  ratios. 

The  area  under  the  hypsometric  curve 

/    J     *  (  J7   )     "     «J  0*y*l  (1*9) 

where  y  «    — 


or 


h 


1         * 

jfT  J     a  dh-oj  (lio-a) 


so,  the  land  mass     is 


h 
/  a  dh     =     (V     AH»  (^9-b) 

o  y 

and  the  corresponding  storage  V  is 

V  =  (y  -  o£)AH'  »  0y  AH'  (1^9-c) 

h        V 
3.  Plotting  -=7  vs.  jgy  on  log- log  paper,  a  straight  line 

was  obtained  for  the  watersheds  under  study.   The  formula  of  the 

straight  lin«  is 

|r»c(I|7)V  (1,9-d) 

where  v  is  defined  as  the  valley  shape  coefficient. 

(22) 
Multiple  Correlation v   ' 

Multiple  correlation  is  a  statistical  method  to  find  the 

relationship  between  one  dependent  variable  and  a  number  of  independent 

variables.   If  a  linear  relationship  exists,  the  method  of  fitting  is 
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to  make  the  sum  of  the  squares  of  the  deviations  of  actual  observations 
from  the  theoretical  linear  relation  to  be  a  minimum.   This  is  called 
the  method  of  least  squares. 

The  general  formula  for  multiple  correlation  and  the  general 
formula  for  estimating  the  variance  are  as  expressed  in  Eqs.  (19), 
(20),  and  (21). 

HYDROGRAPH  STUDY  FOR  SMALL  WATERSHEDS 
IN  INDIANA 

Location  of  8mall  Watersheds  Studied 
Twenty-one  small  watersheds  distributed  throughout  the  state 
were  selected  for  the  hydrograph  study.  The  so-called  small  watershed 
is  defined  as  that  having  an  area  of  less  than  100  square  miles.  A 
map  of  Indiana  showing  the  location  of  the  twenty- one  watersheds 
studied  is  shown  in  Figure  19,and  Table  9  shows  the  names  of  the 
watersheds,  their  drainage  areas  and  their  assigned  study  numbers. 
For  convenience,  these  numbers  will  be  used  hereafter  instead  of  the 
name  of  the  watershed. 

Actual  Storm  Hydrographs 
Data  for  storm  hydrographs  were  obtained  from  the  Indianapolis 
District  Office  of  the  U.  S.  Geological  Survey.  These  data  consisted 
of  the  stream  stage  graphs  and  corresponding  rating  tables  from  which 
the  actual  storm  hydrographs  were  plotted.  For  each  watershed,  five 
or  six  good  hydrographs  were  selected,  based  on  the  hydrographs  with 
the  highest  peak,  caused  by  a  single  storm,  and  with  a  good,  smooth 
recession  curve. 
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tS^° 


SEE     TABLE    I    FOR   WATERSHED    DATA 


FIG.    19  -  LOCATION       OF      GAGING       STATIONS       FOR       STUDIED 


WATERSHEDS        IN        HYDROGRAPH        ANALYSIS 
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Table  9 
Study  Watersheds 


Watershed  Gaging  Station  Watershed 

number  area 


(sq.ml.) 


1  Bean  Blossom  Creek  at  Dolan,  Indiana  100.0 

2  Clifty  Creek  at  Rartsville,  Indiana  88.8 

3  North  Fork  Vernon  Fork  near  Butlerville,  Ind.  87.5 
k  Bean  Blossom  Creek  at  Bean  Blossom,  Indiana  Ik. 6 

5  Bear  Creek  near  Trevlac,  Indiana  7.0 

6  Graham  Creek  near  Vernon,  Indiana  77.6 

7  Brush  Creek  near  Nebraska,  Indiana  11.7 

8  Deer  Creek  near  Putnamville,  Indiana  59. 0 

9  Lawrence  Creek  at  Fort  Benjamin  Harrison,  Ind.  2.86 
10  Mud  Creek  at  Indianapolis,  Indiana  1*2. 5 

U  Hinkle  Creek  near  Cicero,  Indiana  16.3 

12  Little  Cicero  Creek  near  Arcadia,  Indiana  kk.f 

13  Buck  Creek  near  Muncie,  Indiana  36.7 
l£  Salt  Creek  near  McCool,  Indiana  78.7 

15  Little  Calumet  River  at  Porter,  Indiana  62.9 

16  Cedar  Creek  at  Auburn,  Indiana  93 

17  West  Creek  near  Schneider,  Indiana  5U.5 

18  Iroquois  River  at  Rosebud,  Indiana  30.3 

19  Bice  Ditch  near  South  Marion,  Indiana  22.6 

20  Big  Slough  Creek  near  Collegeville,  Indiana  84.1 

21  Carpenter  Creek  at  Egypt,  Indiana  1*8.1 
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In  studying  the  actual  storm  hydrographs,  it  was  noted  that  the 
time  to  peak  t  does  not  vary  radically  for  the  same  watershed  and 
hence  an  average  time  to  peak  can  be  used  as  a  parameter.  Actually, 
the  time  to  peak  is  a  function  both  of  storm  patterns  and  the  basin 
characteristics.   If  the  storm  pattern  changes  drastically  it  is 
obvious  that  the  time  to  peak  will  also  change.  But,  by  assuming 
that  the  storms  selected  possessed  essentially  the  same  kind  of  patterns, 
then  the  time  to  peak  will  be  essentially  the  same  for  the  same  water- 
shed. 

Dimensionless  Hydrograph 
Using  the  data  obtained  from  the  actual  atorm  hydrographs,  the 
dimensionless  hydrograph  was  obtained  by  plotting  q./<L  against  t/t  . 
It  is  quite  interesting  to  note  that  the  dimensionless  hydrographs 
retain  almost  the  same  shape  for  a  given  watershed.  Therefore,  a 
typical  dimensionless  hydrograph  can  be  used  to  represent  the  hydro- 
graph  shape  for  the  watershed.   By  comparing  this  typical  dimensionless 
hydrograph  with  the  theoretical  hydrograph  shown  in  Figure  16,  the 
Gamma  function  argument  can  be  determined  for  the  watershed  studied. 

Recession  Curve  Analysis  for  Dimensionless  Hydrograph 
The  general  theory  of  the  recession  curve  of  the  hydrograph  in- 
dicates that  it  will  plot  as  a  straight  line  on  semi logarithmic  paper. 
However,  such  a  plot  for  the  hydrographs  studied  reveals  that  there 
is  not  one  straight  line,  but  two  or  three  straight  lines.  This 
can  be  explained  by  the  fact  that  the  flow  comes  from  three  different 
types  of  storage,  which  are  stream  channel  and  valley,  surface  soil, 
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and  ground  water  storage.   Since  the  latter  two  types  of  storage 
constitute  the  lower  part  of  the  recession,  they  can  be  neglected 
without  causing  serious  error  in  the  total  storm  hydrograph.  There- 
fore, the  storgge  coefficient  K   has  been  determined  from  the  first 
part  of  the  recession  curve,  which  is  derived  from  stream  channel 
and  valley  storage.   Table  10  lists  the  hydrograph  parameters:  time 
to  peak  t  ,  storage  coefficient  K.,  and  the  corresponding  Gamma  func- 
tion argument  n  of  the  actual  hydrograph  studied. 

Geomorphological  Study 
Five  quantitatively  measureable  watershed  characteristics  were 
determined  for  seventeen  small  watersheds  where  topographic  maps  were 
available  and  are  listed  in  Table  11. 

Multiple  Correlation 

Different  sets  of  multiple  correlations  were  made  between  the 
dependent  variables  t  and  K..,  and  the  independent  variables  A,  L, 
S,  f  and  v.  Table  12  shows  the  results  of  these  multiple  correla- 
tions. 

From  Table  12,  it  will  be  seen  that  the  valley  shape  coef- 
ficient v  and  the  watershed  shape  factor  f  do  not  provide  any  markedly 
better  degree  of  correlation  than  that  obtained  by  the  use  of  the 
first  three  variables  alone.  That  is  to  say,  if  only  the  three 
variables  A,  L,  aad  S  are  used  for  correlation,  the  result  does 
not  lose  much  in  accuracy  as  compared  to  the  use  of  all  five  factors. 
Since  both  the  valley  shape  coefficient  v  and  the  watershed  shape 
factor  f  can  only  be  obtained  by  rather  tedious  and  time-consuming  work, 


73 


Table  10 
HYDROGRAPH  PARAMETERS 
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and  since  they  do  not  materially  improve  the  correlation,  they  were 
deleted  from  further  consideration.  Therefore,  the  three  independent 
variables  A,  L,  and  S  were  selected  for  use  in  estimating  the  hydro- 
graph  characters  t  and  K..    Figures  20  and  21  were  derived  by 
plotting  t  and  K  against  the  watershed  characteristics  A,  L,  and  S, 
based  on  the  regression  formulas: 

tp  -  5l.te  A1-08'  ifl.233  s-0.668  (5Q) 

and 

h    -780A°-937L-1.^   8-l.*73  (51) 

Straight  lines  were  fitted  by  the  method  of  least  squares. 

Working  Charts  for  Estimating  the 
Hydrograph  Parameters  t  and  K. 

The  multiple  correlation  analysis  shows  the  relationship  between 

t  and  K.  and  the  watershed  characteristics:  therefore,  t  and  K.  can 

pi  *     p      1 

be  calculated  by  the  regression  formulas,  Eqs.  50  and  51.   Since  the 
varying  exponents  render  the  formulas  somewhat  complicated  for  practical 
use,  working  charts  for  t  and  K,  are  shown  in  Figures  22  and  23. 

An  example  is  presented  to  illustrate  the  use  of  the  charts,  as 
shown  in  dotted  lines  on  Figures  22  and  23.   The  charts  are  difficult 
to  read  for  an  area  of  less  than  five  square  miles,  so  it  is  suggested 
that  the  regression  formulas  be  used  for  these  small  areas. 
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RELATIONSHIP         BETWEEN       tp       AND 
BASIN        CHARACTERISTICS       A  ,  L  ,  a  S 
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FIG.     21    -  RELATIONSHIP        BETWEEN        K,       AND        BASIN 

CHARACTERISTICS  A,    L,    ft     S 
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FIGURE      22    -  WORKING         CHART         FOR        REG- 

RESSION FORMULA,     EQ.    (50) 
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FIGURE      23  - 


WORKING        CHART       FOR         REG 
RESSION        FORMULA,      EQ.    (51) 
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Example :       Watershed  No.  1 

Drainage  area  (A)         100  square  miles 

Length  of  main  stream  (L)  28  miles 

-k  k 

Slope  of  main  stream  (S)    x  10     9.8^  x  10 

The  values  of  t  and  K.  are  found  in  Figure  22  and 
Figure  23  to  be  17.2  and  ik.'J   hours,  respectively. 


Thus  far,  this  study  has  presented  the  procedure  for  determining 
the  shape  of  the  hydrograph  by  measuring  certain  watershed  character- 
istics from  the  topographic  maps,  that  is,  to  determine  the  hydrograph 
parameters  t  and  K,  from  the  multiple  correlation  working  charts 

from  the  known  values  of  A,  L,  and  S:  and  from  the  ratio  of  Kn /t  , 

'      '  1'   p' 

to  determine  the  gamma  function  argument  n  which  gives  the  shape  of 
the  instantaneous  hydrograph. 

Storm  Rainfall  and  Runoff  Relationship 

It  has  been  demonstrated  that  the  shape  of  the  hydrograph  can  be 

determined  from  certain  measurable  watershed  characteristics  and  hence 

that,  where  topographic  maps  are  available,  the  shape  of  the  hydrograph 

can  be  found  for  ungaged  areas.  The  remaining  factor  to  be  determined 

in  order  to  apply  the  hydrograph  is  that  of  the  total  runoff  "R" .  From 

Eq.  ko,   it  is  noted  that  the  total  runoff  R  is  directly  proportional 

to  the  rate  of  discharge   An  accurate  estimate  of  the  total  runoff  R 

is  therefore  of  equal  importance  to  that  of  the  proper  shape  of  the 

/ 
hydrograph.   A  study  was  therefore  made  to  determine  the  relationship 

between  storm  precipitation  and  storm  runoff,  involving  the  analysis 

of  eighty  storm  hydrographs.  The  study  was  based  on  the  following: 
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(a)  The  storm  precipitation  diagram  and  the  resulting  runoff 
hydrograph  as  shown  on  Fig.  2k 

where 

P  ■  total  storm  precipitation  »  P  +  P  . 

la  X 

P  ■  the  amount  of  precipitation  in  inches  lost  initially  due  to 
evaporation,  infiltration,  and  depression  storage,  no  runoff 
occurs. 

P  »  the  amount  of  precipitation  in  inches  after  the  time  of 
beginning  of  runoff. 

P  «  the  effective  precipitation  which  produces  runoff  =  R. 

R  ■  total  runoff,  in  inches,  obtained  by  measuring  the  area 
under  the  actual  storm  hydrograph. 

f  «  rate  of  infiltration  in  inches  per  hour. 

(b)  The  runoff  coefficient  r. 

The  runoff  coefficient  r  is  defined  as  the  ratio  of  the 
precipitation  P  to  the  total  runoff  R.   Since  the  total  runoff  R  is 
equal  to  P  minus  the  losses  occurring  after  the  beginning  of  runoff, 
this  relation  can  be  expressed  as  follows: 

R  =  P  -  Losses  (52) 

If  the  infiltration  loss  be  considered  as  the  major  loss,  then 
t 
R  «  P  -  /  f .  dt  -  Other  Losses  (52-a) 

t 

o 

Also,  if  the  evaporation  and  interception  losses  are  neglected, 

retaining  only  depression  storage,  then 

t 

R  =  P  -  /  f  dt  -  repression  Storage  (52-b) 

X  t    1 
o 

Letting  d  =  the  depression  storage,  the  formula  becomes, 
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FIG.    24  -  STORM         PRECIPITATION         DIAGRAM         AND         THE 


RESULTING        RUNOFF         HYDROGRAPH 


8k 


t 
R     -     P     -     /f    dt.d  (52-c) 

t       x 
o 

then,  the  runoff  coefficient  r 

R                       '  f i  dt  d 

r  -    f    -  1  -       — | f  (53) 

X  XX 

This  relation  indicates  that  the  runoff  coefficient  is  a  func- 
tion of  the  infiltration  rate,  the  depression  storage,  the  size  of 
storm  and  the  base  length  of  the  hydrograph.   If  it  is  assumed  that 
the  last  two  terms  vary  only  over  a  very  small  range,  or  are  a  constant 
for  a  given  watershed,  then  the  runoff  coefficient  will  be  very  nearly 
a  constant. 

(c)  Determination  of  the  runoff  coefficient. 

Values  of  the  runoff  coefficient  r  were  determined  for  the  small 
watersheds  studied.  Analysis  of  eighty  storm  hydrographs  demonstrates 
that  the  runoff  coefficient  does  vary  among  the  several  watersheds. 
Since  the  studied  watersheds  are  small  and  the  rainfall  gages  are  not 
closely  and  evenly  distributed,  it  is  difficult  to  determine  the  true 
average  precipitation  for  a  given  storm  over  a  given  watershed,  and 
hence  difficult  to  determine  the  value  of  P  .  Therefore,  the  derived 
runoff  coefficients  r  are  not  taken  as  a  fixed  constant,  but  rather  as 
falling  within  a  certain  range,  say  0.1-0.3,  or  0.5-0.7,  whichever 
seems  to  be  the  more  reasonable. 

The  runoff  coefficient,  as  defined,  is  the  ratio  of  R  and  P  . 
P  is  the  amount  of  rainfall  after  the  beginning  of  runoff,  the  effect 
of  antecedent  moisture  conditions  having  been  eliminated. 
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The  runoff  coefficients  were  determined  in  this  manner  for  all  of 
the  21  small  watersheds.   If  precipitation  losses  subsequent  to  the 
beginning  of  runoff  are  taken  as  being  mainly  due  to  infiltration,  the 
runoff  coefficient  r  must  have  some  relation  to  the  soil  physics  of  the 
watersheds.  Studying  the  general  soil  regions  in  Indiana  and  their 
subsoil  permeability,  it  is  found  that  the  runoff  coefficient  is  highly 

correlated  with  the  permeability  of  the  soil.  A  soil  survey  map  made 

(31) 
at  Purdue  University  w   seems  adequate  for  practical  use  for  this 

purpose.  This  map  shows  the  distribution  of  soil  regions  and  their 
subsoil  permeability  in  Indiana,  and  is  shown  in  Figure  25.   It  indi- 
cates the  principal  soil  types  in  each  region  and  the  degree  of  subsoil 
permeability  ranging  from  very  highly  permeable  to  very  slowiy  permeable. 

The  relation  of  the  runoff  coefficient,  type  of  soil  and  permea- 
bility is  shown  in  Table  15.  Since  the  relation  is  logical  and  con- 
sistent, the  runoff  coefficient  can  be  estimated  from  a  knowledge  of 
the  soil  types  of  the  watershed.   Hence  by  locating  a  given  watershed 
on  the  soil  map,  the  runoff  coefficient  can  be  readily  determined. 
Table  Ik   lists  the  recommended  runoff  coefficients  for  various  types 
of  soils  for  the  runoff  design  of  small  watersheds  in  Indiana. 
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LEGEND 

A,(H)-VERY    PERMEABLE 

O.H.O- MOSTLY    PERMEABLE 


%    (BLC,E,G,M,P-  MODERATELY    PERMEABLE 
~]    K,L,N- SLOWLY    PERMEABLE 

t:::::::::|  b,f,i,j-very  slowly  permeable 

(    _    )  STUDY  watersheds 

ADAPTED     FROM     HIGHWAY     RESEARCH     BULLETIN    NO-  10,   PURDUE    UNIVEBSITY 


FIG.     25    -  SUBSOIL       PERMEABILITY        MAP        FOR         INDIANA 
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Table  13 
The  Runoff  Coefficient,  Type  of  Soil,  and 
Degree  of  Permeability  of  Subsoil 


Watershed       Runoff 
number        Coefficient 


0.5  - 

0.7 

0.5  - 

0.7 

0.8  - 

1.0 

0.8  - 

1.0 

0.6  - 

0.8 

0.8  - 

1.0 

0.8  - 

1.0 

0.4  - 

o.6 

0.6  - 

0.8 

0.5  - 

0.7 

0.7  - 

0.8 

0.6  - 

0.7 

0.5  - 

0.6 

0.2  - 

0.3 

0.4  - 

0.5 

0.5  - 

0.7 

0.3  - 

0.4 

o.l  - 

0.3 

0.1*  - 

0.6 

0.1  - 

0.3 

0.5  - 

0.7 

Type  Degree 

of  of 

soil**       Permeability** 


1  0.5-0.7  M,  L,  I  Moderately  and  slowly 

2  0.5-0.7  G  Moderately  and  slowly 

3  0.8  -  1.0  J  Slowly 

4  0.8  -  1.0  L  Slowly 

5  0.6  -  0.8  L  Slowly 

6  0.8  -  1.0  J  Slowly 

7  0.8  -  1.0  J  Slowly 

G,  I  Moderately  and  slowly 

9  0.6  -  0.8  E  Moderately 

10  0.5-0.7  E  Moderately 


11  0.7  -  0.8  E  Moderately 

12  0.6-0.7  E  Moderately 

13  0.5-0.6  E  Moderately 

14  0.2-0.3  A,  F         Mostly 

15  O.k   -  0.5  A,  F         Mostly 

16  0.5-0.7  F  Slowly 

17  0.3-0.4  A,  F         Mostly 

18  0.1-0.3  A  Very 

19  0.4-0.6  C  Moderately 

20  0.1-0.3  A,  C         Very 

21  0.5-0.7  B,  C         Moderately 


**  Refer  to  Figure  25. 
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Table  Ik 

Recommended  Runoff  Coefficients 

for 
Various  Types  of  Indiana  Soils 


Type  of  Soil 

A, 

H 

D, 

H, 

0 

c, 

E, 

G, 

M, 

K, 

L, 

N 

B, 

I, 

J 

F 

Runoff  Coefficient 


0.50 
0.50 
0.70 
0.80 
1.00 
0.50-0.80  * 


*  The  F  type  of  soil  is  as  slowly  permeable  as 
types  B,  I,  and  J,   but  the  losses  due  to  depression 
storage  are  quite  different.  There  are  svaraps  and 
lakes  in  this  region  where  the  amount  of  runoff 
passing  through  the  outlet  is  decreased. 
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DESIGN  STORM  RAINFALL 
Since  there  is  a  relationship  between  rainfall  and  runoff,  as 
demonstrated  in  the  preceding  section,  one  might  expect  to  estimate  the 
runoff  for  a  certain  area  through  the  design  rainfall.  Recent  (1961) 
data  as  to  precipitation  depth- duration- frequency  relations  can  be  found 
in  Technical  Paper  No.  Uo,  published  by  the  Weather  Bureau.   Figures  26 
and  27,  which  are  based  on  data  from  Weather  Bureau  Technical  Paper 
No.  UO,  show  the  six-hour  duration  rainfall  for  return  periods  of  25 
and  50  years.   Table  15  gives  the  applicable  ratios  to  convert  the 
six- hour  duration  rainfall  to  other  durations. 

For  the  design  of  many  small  hydraulic  structures,  a  design 
runoff  having  return  periods  of  25  to  50  years  would  seem  to  be  ade- 
quate. The  remaining  question  is  whether  the  return  period  of  runoff 
can  be  taken  as  equivalent  to  the  return  period  for  rainfall.  This 
problem  is  considered  in  a  subsequent  section  of  this  study. 
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ADAPTED    FROM     u  S  W  B. 
TECHNICAL     PAPER     NO     40 


FIG      26    -  25-YEAR,      SIX        HOUR       RAINFALL      IN       INCHES 

FOR'     INDIANA 


V 
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4  3 


ADAPTED    FROM     U.S.W.B 
TECHNICAL    PAPER    NO.    40 


FIG.     27-  50-  YEAR,     SIX       HOUR       RAINFALL       IN        INCHES 

FOR       INDIANA 
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PROCEDURES  OF  STORM  HYDROGRAPH  DESIGN 
From  the  above  studies,  the  general  procedure  of  storm  hydrograph 
design  can  be  outlined  as  follows: 

1.  Study  of  watershed  characteristics. 

The  delineation  of  the  watershed  on  the  topographic  map  and 
the  determination  of  the  watershed  area  in  square  miles  (A),  the  length 
of  main  stream  in  miles  (L),  and  the  mean  slope  of  the  main  stream 
(S  x  10  )  are  the  first  steps  in  the  hydrograph  design. 

2.  Determination  of  the  hydrograph  parameters  t  and  K. . 
The  time  to  peak  t  and  the  storage  coefficient  K  can  be 

determined  from  the  multiple  correlation  charts,  Figures  22  and  23, 
or  calculated  from  the  regression  formulas,  Eqs.  50  and  51,  if  the  area 
of  the  watershed  is  less  than  five  square  miles. 

3.  Determination  of  the  shape  of  the  instantaneous  hydrograph. 
Calculate  the  ratio  of  K./t  ,  find  the  Gamma  function  argument 

n  from  Figure  17,  and  determine  the  shape  of  hydrograph  by  the  "n"  from 
Figure  16.  A  dimensionless  instantaneous  hydrograph  can  be  plotted. 

k.     Determination  of  the  runoff  coefficient. 

Locate  the  given  watershed  on  the  soil  map,  Figure  25,  and 
select  the  runoff  coefficient  by  reference  to  Table  Ik. 

5.  Design  storm  rainfall. 

The  amount  of  storm  rainfall  in  inches  with  return  periods 
of  25  and  50  years  and  various  durations  can  be  found  in  Figures  25  and  27 
and  Table  15.  The  selection  of  the  duration  of  storm  rainfall  is  based 
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Table  15 
Factors  for  Conversion  of  Six-Hour  Rainfall 
Duration  to  Longer  Duration 


Duration  hours  Ratio* 

6  1.000 

7  1.035 

8  1.065 

9  1.090 

10  1.115 

11  l.l4o 

12  I.160 

13  1.185 

14  1.200 

15  1.220 

16  1.235 

17  1.255 

18  1.270 

19  1.280 

20  1.300 

21  1.315 

22  1.325 

23  1.340 

24  1.350 

25  1. 360 

26  1.375 

27  1.385 

28  1.395 

29  1.410 

30  1.420 

31  1.425 

32  1.435 

33  1.445 
31*  1.455 

35  1.465 

36  1.470 

*  From  the  Engineering  Handbook,  Hydrology, 
Soil  Conservation  Service,  U.S.D.A. 
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on  the  time  of  concentration  of  the  watershed,  considering  that  a 
storm  rainfall  duration  equal  to  the  time  of  concentration  will  result 
in  the  maximum  rate  of  discharge.   Since  the  time  of  concentration  of 
the  watershed  is  difficult  to  determine,  one  may  assume  the  design 
rainfall  is  uniformly  distributed  over  the  entire  watershed,  in  which 
case  the  time  of  concentration  is  equal  to  the  time  to  peak.   It  is 
therefore  suggested  that  the  design  storm  duration  be  made  equal  to  the 
value  of  t  ,  which  may  be  found  from  Figure  22.   For  conservative 
design,  it  is  suggested  that  a  six- hour  duration  be  used  as  a  minimum  for 
times  to  peak  less  than  six  hours,  and  that  the  duration  be  made  equal  to 
the  time  to  peak  when  it  is  longer  than  six  hours. 

£.  Determination  of  total  runoff. 

The  design  runoff  can  be  determined  from  the  formula 

R  =  Px-r  (5*0 

that  is,  the  design  storm  rainfall  times  the  runoff  coefficient. 

7.  Computation  of  maximum  discharge. 

The  maximum  discharge  can  be  computed  from  Eq.  ko   when  t  =  t  . 

AR 
a     =  -r-   f(n,t  ),  where  q   is  the  instantaneous  hydrograph 

P      P 
peak  discharge  and  is  a  good  estimate  of  the  maximum  discharge  for  small 

watersheds.  The  values  of  f(n,t  )  for  various  values  of  the  Gamma 

'   P 

function  argument  n  are  as  follows. 
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f(n,tp) 


l.k 

1.6 
1.8 
2.0 
2.2 
2.5 
3.0 
4.0 
5.0 

6.0 
7.0 
8.0 
9-0 
10.0 


0.210 
0.272 
0.323 
0.368 
0.^09 
0.U63 
0.5^0 

0.672 
0.782 
0.878 
0.972 
1.0U1 
1.150 
1.180 


8.  Plotting  the  storm  hydrograph. 

From  the  dimensionless  hydrograph,  the  time  to  peak  t  ,  and 
the  maximum  discharge  q  ,  the  instantaneous  hydrograph  can  be  easily 
plotted,  which  for  small  watersheds  may  be  taken  as  the  runoff  hydro- 
graph.  Figure  28  is  a  sketch  diagram  showing  the  sequence  of  the  steps 
employed  to  obtain  the  design  hydrograph.  An  example  of  the  computation 
of  a  design  hydrograph  is  shown  as  follows:  Watershed:  Pleasant  Run 
at  Arlington  Avenue,  Indianapolis,  Indiana. 
Watershed 


characteristics :   Drainage  area  (A) 

Length  of  main  stream  (L)       ^ 
Mean  slope  of  main  stream  (S)  10* 


Hydrograph 
parameters: 


From  Figure  22 
From  Figure  23 


7.67  square  miles 
3.6  miles. 
32. k   x  10" 

t  =5.8  hours 
K?  «  4.7  hours 


Gamma  function  argument  n: 

Since  K  /t  =  0.81,  from  Figure  17,  n-5;  from  the  table  in  step  7 


f(n,t  )  =  0.782. 


Design  rainfall,  Figure  26: 
25-year,  6-hour  rainfall 


P  «  3.5  inches 
x 
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Weather  Bureou 
Tech   Poper  No    40 


Design    Rolnfoll 


Subsoil     Permeobili 
Mop 


Runoff    Coefficient 


Topogrophic    Mop 


Design    RunoH 
\ 


/ 


Watershed 
Characteristic! 
4.   L.   S 
/ 


Hy drogf oph 
Pofometers 


Maximum    Discharge  "    t        Gommo    Function 

Argument     n 


Design    Hydrograph  Shape  of   the  Dtmensionless 

Hydrograph 


1m 


FIG.    28   -  SEQUENCE       OF        COMPUTATIONS        TO       DESIGN 


STORM        HYDROGRAPH        FOR       SMALL       WATERSHED 
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Design  runoff: 

From  Figure  25  and  Table  14  Runoff  coefficient  r  =  0.70 
Runoff  R  -  0.7  x  3-5  ■  2.1+5  inches 

Maximum  discharge,  q 

AR  Wn  t    \         7.67  x  640  x  2.45    w  n  7ft5 

%a  '     t"  f(n'V         5^ X  °*782 

P 

=  l64o  cfs. 


Design  Hydrograph: 
l64o  cfs, 


a 
-m 


5.8  hours 


n=5 


Dimensionless  hydrograph  * 


Storm  hydrograph 


t/t 


q/qm 


'(br.) 


L(cfs) 


0.2 

3.96 

0.4 

28.4o 

0.6 

64.40 

0.8 

91.10 

1.0 

100.00 

1.2 

93.^ 

1.4 

77.70 

1.6 

59.60 

1.8 

43.00 

2.0 

29.40 

3.0 

2.68 

4.0 

1.54 

1.16 

65 

2.32 

465 

3.48 

1,057 

4.64 

1,495 

5.80 

l,64o 

6.96 

1,530 

8.13 

1.274 

9-30 

976 

10.4o 

704 

11.60 

482 

17.40 

44 

23.20 

25 

■"The  dimensionless  hydrograph  can  be  obtained  from  Figure  16, 
or  from  Appendix  B-2.  The  derived  design  hydrograph  is  shovn  on 
Figure  29. 
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FIG.      29    -  DERIVED       DESIGN      HYDROGRAPH  ,     PLEASANT       RUN 

AT      ARLINGTON       AVENUE,     INDIANAPOLIS 
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Determination  of  Return  Period  at 
a  Given  Peak  Discharge 

This  study  also  provides  a  means  for  solving  another  kind  of 

practical  engineering  problem.  This  is  that  of  the  determination  of  the 

return  period  of  a  given  peak  discharge.  This  is  accomplished  by  a  more 

or  less  inverse  procedure  to  determine  the  frequency.  An  example  is  given 

as  follows: 

Watershed:   Pleasant  Run  above  State  Road  100,  Indianapolis 

Watershed  area  (A )  3. 31  square  miles 

Length  of  main  stream  (L)        .  2.k   miles. 

Mean  slope  of  main  stream  (S)x  10"  32.5  x  10" 

From  Eqs.  50  and  51 

Kl 
t  =  3.85      Kx  =  3-87        ^  -  1.0 


P 


From  Figure  17 


n=4  ^(n,t  )  -  0.672 


If  q^  =  2280  cfs, 

AD 

From  Eq.  i+0  q^  =  —     f(n,t  ) 

P 


2280  =     3<31X6^RXO'672 


R     6 
Runoff  coefficient     r  =  0.70         Px  *  r  =  0~7~  =  8'6' 

Design  Duration        ■  6  hr. 

Since  100  -  yr.6-hr.P   =  k.k" 

The  return  period  of  $  »  2280  cfs  is  more  than  100  years 
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CONCLUSIONS  AND  DISCUSSION 

1.  The  methodology  developed  in  this  study  for  the  design  of  a 
storm  hydrograph  is  semi- theoretical  and  semi- empirical.  The  derivation 
of  the  dimen6ionless  hydrograph,  the  Gamma  function  argument  n  relating 
to  the  shape  of  the  hydrograph,  and  the  recession  curve  study  developing 
the  relations  between  t  ,  K  and  n  are  all  theoretical.  The  multiple 
correlation  between  the  hydrograph  parameters  t  ,  K..  and  the  watershed 
characteristics,  and  the  study  of  rainfall  and  runoff  relationships 
relating  to  different  soil  properties  in  order  to  determine  the  runoff 
coefficient  are  based  on  empirical  data. 

2.  This  study  considered  rainfall  with  different  durations  and 
frequencies,  soils  of  different  permeabilities,  and  the  various  shapes 
of  the  hydrograph  resulting  from  the  different  geomorpho logical  features 
of  the  watersheds. 

3.  For  convenience  in  practical  engineering  design,  this  study 
has  been  directed  toward  making  the  design  procedure  as  simple  as  pos- 
sible.  Most  of  the  required  data  can  be  obtained  from  topographic 
maps,  and  from  the  charts  and  tables  presented  herein.  The  whole  design 
does  not  require  more  than  three  hours  to  complete. 

k.     Since  the  small  watersheds  used  in  this  study  range  in 
area  from  2.86  to  100  square  miles,  the  use  of  the  procedures  developed 
herein  for  watersheds  having  an  area  of  less  than  two  square  miles  or, 
say,  1000  acres,  is  questionable. 
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5.  As  a  check  on  the  reliability  of  the  results  obtained  from 
the  application  of  these  procedures,  values  of  a     were  computed  for  11 
small  watersheds  for  which  flood  magniture  and  frequency  studies  had 
previously  been  made  in  the  frequency  analysis.  Results  obtained  by 
the  two  methods  may  be  compared  from  the  data  in  Table  16. 

Table  16  shows  that  the  maximum  discharge  obtained  by  these  two 
methods  agrees  quite  well,  except  for  one  or  two  cases  in  which  the 
disagreement  may  be  due  to  poor  estimates  of  the  runoff.  In  general, 
however,  the  method  developed  in  this  report  seems  to  give  results 
which  are  adequate  for  the  design  of  hydrographs  for  small  ungaged 
watersheds. 

6.  Since  this  study  is  based  on  the  assumption  that  the  runoff 
frequency  may  be  taken  as  equivalent  to  rainfall  frequency,  the  question 
may  well  be  raised  as  to  the  validity  of  this  assumption.  It  is  under- 
stood that,  due  to  different  antecedent  moisture  conditions  and  other 
factors,  the  100-year  return  period  rainfall  might  very  well  not  produce 
a  100-years  return  period  runoff,  or  that  a  25-year  rainfall  might  cause 
a  100- year  flood  under  favorable  conditions.  However,  the  situation 

in  this  study  is  somewhat  different,  in  that  the  design  rainfall  con- 
sidered herein  is  taken  only  as  that  which  occurs  subsequent  to  the 
beginning  of  runoff,  thereby  eliminating  the  effect  of  antecedent 
moisture  conditions.  Therefore,  one  might  assume  that  the  frequency  of 
rainfall  can  be  used  as  the  frequency  of  runoff  without  causing  much 
error.   It  is  also  shown  in  Table  16  that  the  use  of  rainfall  fre- 
quency in  the  design  hydrograph  to  find  the  maximum  discharge  agrees 
quite  well  with  the  maximum  discharges  obtained  from  the  study  of  runoff 
frequency . 
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Table  16 

Comparison  of  results  of  the  maximum  discharge 
obtained  from  this  design  hydrograph  study  and 
the  frequency  study 


Watershed 

100- 

-years 

25-yeare 

number 

Design 
hydrograph 
(cfs) 

Frequency 
study 
(cfs) 

Design 
hydrograph 
(cfs) 

Frequency 
study 
(cfs) 

1 
2 

3 
Ik 

15,^00 

12,600 

28,500 

3,^50 

16,000 

17,000 

27,500 

5,840 

15,800 

10,500 

25,600 

2,790 

12,000 

15,000 

22,000 

2,950 

15 
16 
17 
18 

4,720 

2,570 

2,600 

575 

4,500 

2,000 

2,600 

570 

5,700 

2,160 

2,100 

1+65 

5,500 

1,620 

2,100 

465 

19 
20 

21 

1,750 
3,100 
6,380 

1,000 

3,100 
5,6oo 

1,1+00 
2,480 
5,000 

830 
2,1+60 
4,200 
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7.  As  mentioned  before,  the  design  runoff  R  is  based  on  the 
design  rainfall  and  the  runoff  coefficient  corresponding  to  the 
watershed  location.  Obviously,  the  worth  of  the  derived  design 
hydrograph  hinges  in  large  measure  upon  the  estimates  of  the  value  of 
R.  Since  the  runoff  coefficient  is  not  a  fixed  value,  the  estimate 
of  the  total  runoff  may  well  vary  with  the  Judgment  of  the  individual. 
The  suggested  runoff  coefficients  in  Table  Ik   are  considered  to  be 
conservative. 
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PART  III 


HYDRODYNAMICS  OF  OVERIAND  FLOW 
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INTRODUCTION 
The  theoretical  determination  of  the  peak  discharge  and  of  the 
synthetic  hydrograph  for  small  watersheds  must  be  based  on  the  hydro- 
dynamics of  overland  flow  or  sheet  flow.  However,  a  theoretical 
analysis  is  possible  only  for  small  areas  with  relatively  uniform 
slope.   Larger  watersheds  which  have  a  more  complicated  topography 
must  be  divided  in  small  zones  or  small  subwatersheds  in  which  each 
zone  is  subjected  to  an  overland  flow  process.  The  total  outflow 
may  be  obtained  by  considering  those  small  subwatersheds  in  series 
and  routing  their  outflow  through  the  channel  to  the  outlet.  This 

latter  operation  is  similar  to  the  process  of  flood  routing  which 

(32) 
makes  use  of  the  equations  of  unsteady  flow  in  open  channels. 

The  relation  between  rainfall  and  runoff  for  overland  flow  is 
derived  from  the  basic  equations  of  hydrodynamics:  the  momentum 
equation  and  the  equation  of  continuity.   They  contain  the  rainfall 
as  inflow,  the  depth  of  overland  flow  and  its  mean  velocity  are  the 
dependent  variables,  the  independent  variables  are  the  distance  of 
overland  flow  and  time.  A  set  of  two  qua si- linear  partial  differen- 
tial equations  of  hyperbolic  type  is  obtained.   They  may  be  solved 
by  the  method  of  characteristics. 

This  study  attempts  to  develop  the  outflow  hydrograph  from  a 
design  storm  precipitated  on  a  given  area  of  uniform  slope  and  rough- 
ness.  The  relation  between  design  storm  inflow  and  the  runoff  will 


io6 


yield  the  conditions  under  which  the  unit  hydrograph  assumption  of 
linear  superposition  made  by  Sherman  in  1952  may  be  used  as  a  good 
engineering  approximation. 

THEORETICAL  ANALYSIS 
The  theoretical  analysis  of  overland  flow  may  be  considered  in 
three  parts: 

a.  laminar  sheet  flow  approximate  solution 

b.  turbulent  sheet  flow  approximate  solution 

c.  general  hydrodynamic  equations  of  overland  flow. 

The  laminar  sheet  flow  occurs  at  the  early  stage  of  overland  flow 
where  the  depth  and  velocity  of  flow  are  small.  When  the  Reynolds 
number  of  overland  flow  reaches  the  magnitude  somewhere  larger  than 
1,000,  it  may  be  considered  as  turbulent  flow. 

Laminar  Sheet  Flow  Approximate  Solution 
There  is  no  definite  upper  critical  value  of  Reynolds  number 
of  the  laminar  flow.  Woo*5  '  in  1956  found  that  the  upper  critical 
value  of  Reynolds  number  for  sheet  flow  with  rainfall  is  not  a  con- 
stant  but  ranges  from  UOO  to  900.   Levicir  '   in  1962  explains  the 
phenomena  of  flow  of  thin  liquid  film  that  there  are  three  different 
laminar  flow  regimes: 

1.  At  Reynolds  numbers  R  ■  — y  that  do  not  exceed  20  to  30 
(where  u  is  the  average  velocity  over  the  cross  section  of  the  film 
and  y  is  its  thickness),  there  exists  the  usual  viscous  flow  and  the 
film  thickness  is  constant. 
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2.  at  Reynolds  numbers  >  50  to  50,  a  so-called  wave  regime 
appears  in  which  a  wave  motion  is  superposed  on  the  forward  motion 
of  the  film. 

3.  at  Reynolds  numbers  Re  £  1500,  the  laminar  regime  is  re- 
placed by  a  turbulent  motion. 

The  derivation  of  the  approximate  equations  of  laminar  sheet 
flow  is  based  on  the  work  done  by  Izzard.  y'         For  slow  steady 
laminar  sheet  flow,  by  equating  the  component  of  the  gravity  force 
in  the  direction  of  motion  to  the  viscous  drag,  the  average  velocity 
is  found  to  be 

2     S  ®o   2  5 

u  «  ky  =   y   and  q  «=  Ky  (55) 

3  v 
where  v  is  the  kinematic  viscosity  and  S  is  the  slope  of  the  plane. 
The  discharge  at  the  outlet  and  the  water  surface  profile  over  the 
plane  can  be  obtained  as  follows: 

1.  General  sketch  of  laminar  sheet  flow  produced  by  rainfall. 
Figure  (30)  shows  the  laminar  sheet  flow  produced  by  a  given  rainfall, 
y  is  the  depth  of  flow  at  the  outlet,  py  is  the  average  depth  over 
the  plane  where  p  is  a  constant  less  than  one  and  x  is  the  length 

to  the  outlet  from  the  starting  point  where  x  ■  0. 

2.  determination  of  the  value  of  p. 

From  Fig. (30),  the  total  volume  of  detention  over  the  plane  is, 

letting  q  be  the  discharge  at  the  dovnstream  end  of  the  plane: 
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Ve*pyox  «    /  y  dx 
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P(  j*  )  5*      -    /  (  i  )  *i 
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P(^f)3x       -    /     (£*)5dx 


(56) 


3.  The  continuity  equation  is  obtained  for  the  sheet  flow, 
taking  OAB  (Fig.  50)  as  the  control  volume.  The  basic  equation  of 
continuity, 

00A  pV^dA  «  -  |g  /  p  dv  (57) 

which  expresses  that  the  net  mass  influx  is  equal  to  the  rate  of 
change  of  mass  inside  the  control  volume  signifies  that  the  inflow 
minus  the  outflow  is  equal  to  the  change  of  storage.  The  inflow 
in  a  time  dt  is  ixdt.   Assuming  that  equation  (55)  for  the  velocity 
distribution  remains  valid  for  the  unsteady  flow,  as  suggested  by 
Levichv>%  the  outflow  in  a  time  dt  is 

qodt»   Tv  *odt  (58) 

and  the  change  in  storage,  according  to  equation  (56),  is 

dVe  "  I  X  ^o  (59) 

The  equation  of  continuity  thus  becomes 
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8° 
i  x  dt     -      —2    y5  dt  +    |  x  ^  (6o) 


or 


#°     3^+     _ 


(i  x  -  Tv  yo)dt  *  "5 x  dyo  ^6o-fi) 

for  a  finite  time  increment  At,  the  Ay  can  be  obtained  as 

*o-   3(i-  TO  yo)At    *  (6l) 

Equation  (6l)  may  be  used  for  a  numerical  calculation  of  the  depth 
of  overland  flow  at  the  outlet.  The  corresponding  average  velocity 
at  the  outlet  can  be  found  easily  from  Equation  (55).   By  using  the 
different  length  for  x  in  Equation  (6l),  the  overland  flow  profile 
can  thus  be  determined. 

Turbulent  sheet  flow  approximate  solution 
The  approximate  equations  for  turbulent  overland  flow  are 

derived  in  a  similar  fashion  as  those  for  the  laminar  flow.  The 

(36) 
1/7- th- power  velocity  distribution  lawv^' 

v     V   X/7 

2.  m  (  A-  )    where  v,  y'  are  local  velocity  and  depth  and 

U     y   '  ' 

U,  y  are  maximum  velocity  and  depth, 

is  used  to  express  the  relation  between  the  law  of  friction  and 
velocity  distribution.  The  mean  velocity  of  the  turbulent  sheet 
flow  can  be  derived  as  follows: 

1.  For  slow  steady  turbulent  sheet  flow,  by  equating  the 
component  of  the  gravity  force  in  the  direction  of  motion  to  the 
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drag,  ve  have: 

7  y  SQdx  «   Tq  dx  (62) 

2.  Expressing  the  shear  stress  t  by  the  Blasius  empirical 
equation 

to  -  0.0225  p  U7/U  (  ^  )XA  (65) 

where  U  is  the  maximum  velocity.  For  the  1/7- th- power  velocity 
distribution  law,  the  ratio  of  the  mean  velocity  u  to  the  maximum 
velocity  is  0.817,  then  Eq.  (63)  becomes 

tq  «  0.0521  p  u  7/4  (  2  )1A  (64) 

3.  Combining  the  equations  (62)  and  (64),  the  mean  velocity 
can  be  expressed  as 

u  =  k  y5/7  (65) 


and 


q  -  k  y1277  (66) 


where  4/7 

(83  ) 
k-7.15    A  (67) 

v  ' 

The  coefficient  p,  the  ratio  of  the  mean  depth  of  overland 

to  the  depth  of  flow  at  the  outlet,  for  the  turbulent  sheet  flow 

can  then  be  determined  by  referring  to  Fig.  (50) 
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p   s  X353   -  °-651  (68) 

The  equation  of  continuity  can  be  expressed  as 

1  71*5 
ixdt-ky     '   ^  dt  +  0.631  x  dy  (69) 

or 

[ix  -  ky1,Tl5]dt  =  0.631  x  dy  (69-a) 

for  a  finite  time  increment  At  ,   the  Ay   can  be  obtained  as 

Ay  «  1.583  [i  -  I  y1,715]  At  (70) 

Eq.  (70)  has  the  same  form  as  Eq.  (6l)  for  the  laminar  sheet  flow, 
and  it  may  be  used  to  obtain  the  approximate  solution  for  the  over- 
land flow  when  the  Reynolds  number  exceeds  1,000. 

General  hydrodynamical  equations  of  overland  flow 
The  differential  equations  of  the  overland  flow  during  the 
rainfall  is  derived  with  the  following  assumptions: 

1.  The  rainfall  intensity  is  uniform  and  constant. 

2.  The  surface  is  impervious,  or  the  infiltration  rate  is 
assumed  to  be  constant,  so  that  the  effective  rainfall 
rate  is  constant,  that  is,  the  rainfall  rate  minus  the 
infiltration  rate. 
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3.  The  surface  has  an  infinite  width,  so  that  the  unit  width 
is  assumed  in  the  computation. 

k.     The  slope  of  the  surface  is  uniform  and  small. 

5.  The  energy  coefficient  and  momentum  coefficient  are 
assumed  to  be  unity. 

The  physical  features  of  overland  flow  produced  from  a  uniform 
and  constant  rainfall  are  shown  in  Fig.  51,  where  i  is  the  rainfall 
intensity  in  inches  per  hour,  y  is  the  depth  of  overland  flow  in  feet, 
u  is  the  mean  velocity  of  overland  flow  in  the  direction  parallel  to 
the  ground  surface  in  feet  per  second,  v  is  the  component  of  rainfall 
velocity  in  the  direction  parallel  to  the  overland  flow  in  feet  per 
second  and  %is  the  slope  of  the  ground.   By  taking  a  small  section 
a,  b,  c,  d  as  a  control  volume,  the  equation  of  continuity  and  the 
equation  of  motion  can  be  derived  as  follows: 

i.  The  continuity  equation  is  based  on  the  conservation  of  mass. 
The  mass  inflow  minus  the  mass  outflow  is  equal  to  the  rate  change  of 
storage  of  mass  in  the  control  volume.  As  shown  in  Fig.  (52) 


dt 


inflow:     y  u  dt  +  i  dt  dx 

outflow:    (y  +   |^L  dx)(u  +  g  dx) 

rate  change  of     > 

storage:    ^  dt  dx 

the  conservation  of  mass  is  written  as 

y  u  dt  +  i  dt  dx  -  (y  +  |*  dx)(u  +  |jj  dx)dt  -  |*  dt  dx 

and  neglecting  second  order  terms,  it  becomes 
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FIG.    3!     —  DEFINITION  SKETCH      OF      OVERLAND      FLOW 
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du 


idtdx-y^dxdt-u^  cbcdt  a|Z  dt  dx 
thus,  the  continuity  equation  is 

ii.  The  equation  of  motion  is  derived  from  the  momentum  theorem 
as  follows: 

Considering  a  control  volume  a,  b,  c,  d  as  in  Fig.  (33) 
where  F,  and  F  are  hydrostatic  forces 


2 

F  «  y  — 
Pl    ~  2 

7(y  + 

F2- 

2 

the 

frictional  resistance  is 

t  dx  , 

o     ' 

the 

body  force  is 
7   ySQ  dx 

and  the  sum  of  the  forces  is 

2  r(y+  ^**)2 

7  Y    ♦  y    y  SQ  dX  -   g TQ  dX 

or  neglecting  second  order  terms 

-  y  y    |*dx     +7ySodx-Todx 

The  momentum  equation  is 

K   *  ///  B*P  dv  -  #   T(p  u*dA)  ♦  |^  ///"u*(p  dv)  (72) 

cv  cS  cv 

where  the  left  side  is  the  sum  of  the  external  surface  forces  and 


116 

of  the  body  forces,  the  first  term  on  the  right  hand  side  is  the  net 
momentum  flux  through  the  control  surface  and  the  last  term  is  the 
time  rate  of  increase  of  momentum  in  the  control  volume.  For  the 
control  surface  a  b  c  d  of  the  net  momentum  influx  can  be  evaluated  as 

00  ^u(pudA)»p(  u+  ^  dx)2  (y  +  |  dx)  -  u2y  p  -  p  i  v  cbt 

or  neglecting  second  order  terms 

»     p[(u     +  2u  ^  dx)(y  +    g  dx)  -  u2y  -   i  v  dx] 

2      ,    a    du,  2dy.  2 


p[u  y  +  2uyr^  dx  +  u      ^dx-  uy-   iv  dx] 
p[2u  y  |S  dx  ♦  u2  |g  dx     -ivdx] 


pudx[y~     +    |]-pivdx 


and  the  rate  of  increase  of  momentum  in  the  control  volume  is 


-§£  ///  u(p  dv)  -   ^pydx  +  ^[  updx 


The  momentum  equation  is  then 

ry  S  dx  -  y  y  —^-   dx  -  t  dx 

'       o      '  '  dx       o 


■|gpy  dx  ♦  p  yu    -^dx  ♦  pu  |  dx  ♦     pu  f*  dx  -   piv  dx 


By  simplifying  and  rearranging 

,  dy  o  du 

o  "  g     dx  py       "dt       u  dx  y     dx  "     y  dt 


gS     -  g  -4y °-     -    ^  +  u  4^    +    -.    £  +    -A*    -     *= 
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or 


dt     u  dx    B  3£     8  o     py    y  *  dx    dt  ;  *  y 
since 

H  ♦«£  ♦«  I  "  «V   ^°  •*('--)  (75) 

Equations  (71)  and  (7?)  are  the  equation  of  continuity  and  the 

equation  of  motion  for  the  overland  flow  respectively  where  u  and  y 

are  the  dependent  variables,  t  and  x  are  the  independent  variables,  g 

and  p  are  constant  and  S  ,  t  ,  i.  v  can  be  evaluated  for  given 

o   o 

specific  rainfall  and  watershed  condition.  The  Eqs.  (71)  and  (73) 
are  second  order  quasi- linear  partial  differential  equations.  They 
can  be  solved  by  the  method  of  characteristics  if  they  are  of  the 
hyperbolic  type. 

The  differential  equations  governing  the  motion  of  thin  liquid 
film  was  given  by  Levich    .  The  equation  of  motion  of  laminar  flow 
which  contains  the  surface  tension  and  viscosity  wa6  derived  from  the 
Havier-Stokes  equation  under  the  assumption  that  the  velocity  profile 
remains  similar  to  that  of  steady  flow: 

du      q      du     o  dfy     3vu  ,  .* 

Bt   *  %  u  3x"  "  p  -^5  "  —  ♦  ■  (7U> 

and  the  continuity  equation  is 

H  ♦  ^  -  •  (75) 

where  o  is  known  as  surface  tension  and  u  is  mean  velocity.   Since 
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the  change  of  depth  y  vlth  respect  to  length  x  is  small  the  third 
derivative  of  the  depth  with  respect  to  x  is  very  small,  and  the 
term  containing  the  surface  tension  may  be  neglected.  With  the  rain- 
fall inflow  and  rainfall  momentum  terms  added  to  the  above  equations, 
the  differential  equations  of  overland  laminar  sheet  flow  for  the 
earlier  stage  of  runoff  can  be  expressed  as  follows: 

The  equation  of  motion 


y 


y 


and  the  equation  of  continuity 

The  above  equations  (76),  (77)  are  of  the  parabolic  type,  and  are 
valid  only  for  the  initial  stage  of  overland  flow  when  the  Reynolds 
number  is  less  than  1000. 
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PHYSICAL  PHENOMENA  OF  OVERLAND  FLOW 

_   __  _       (57) 
Runoff  Phenomena 

In  general,  the  runoff  phenomena  produced  by  a  given  precipita- 
tion may  be  represented  as  shown  in  Fig.  34.    This  represents  a 
half  section  of  a  small  drainage  basin  with  precipitation  taking 
place  at  a  rate  exceeding  the  infiltration  capacity  of  the  soil. 
The  precipitation  excess,  at  first,  flows  into  depressions,  and  at 
the  time  when  these  depressions  are  filled  and  begin  to  overflow  there 
will  be  some  depth  of  initial  detention.  When  the  initial  detention 
is  filled,  the  precipitation  excess  or  supply  will  flow  toward  the 
outlet.  The  time  required  to  build  up  the  initial  detention  is  called 
the  initial  time  lag  which  can  be  defined  as  the  time  elapsed  from  the 
beginning  of  precipitation  excess  until  active  surface  runoff  begins. 
The  initial  detention  and  the  initial  time  lag  can  be  expressed 
graphically  as  shown  in  Figure  (55),  where  P  is  the  total  precipitation, 
F  is  the  accumulated  infiltration  loss,  Q  is  the  mass  runoff  and 

y  and  t  are  the  initial  detention  and  the  initial  time  lag  respec- 
o      o 

tively.  This  shows  that  the  initial  detention  and  the  initial  time 
lag  are  two  existing  items  in  the  physical  phenomena  of  surface  runoff 
and  can  be  evaluated  from  the  experiment  test.   These  will  depend  on 
the  intensity  of  the  precipitation,  the  infiltration  rate,  the  slope 
and  the  roughness  of  the  watershed.   These  values  may  be  very  small, 
but  these  serve  to  evaluate  the  initial  condition  which  is  a  necessary 
condition  to  solve  the  derived  differential  equations  for  overland 
flow. 
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FIG.    34-  HALF- SECTION       OF      A       SMALL       DRAINAGE        BASIN 

ILLUSTRATING       RUNOFF       PHENOMENA 
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Evaluation  of  tractive  force 
One  of  the  most  important  factors  in  the  overland  flow  study  is 
the  retarding  tractive  force  on  the  ground  surface.   It  can  be  ex- 
pressed by  the  term  t  dx,  or  as  a  function  of  head  loss  due  to  the 

friction  slope  S„.  The  relationship  between  r     and  S  is 
1  of 

or,  for  overland  flow,  as  R  —  y 


\>  -  7  y  sf  (79) 


Since  the 


hf 
Sf  =  ~\  (80) 

hf  -  f  T   k  (do  ■  *R)  (81) 

o 

the  tractive  force  t  can  be  shown  as  a  function  of  friction  coeffi- 
o 

cient  f  and  the  mean  velocity  of  the  overland  flow,  it  is 

v"  r  (82) 

The  friction  coefficient  f  in  the  laminar  region  can  be  defined  by  a 
general  equation 

r--i_  (83) 

(52) 
where  K  is  a  constant,  Re  is  the  Reynolds  number.  Chow'    indicates 

that  Re  <  500  is  the  condition  for  laminar  flow  to  exist  in  open 

channels  and  for  which  there  is  a  linear  relationship  between  the  f 

(5  a) 
and  Re.   Straub  and  others    °  found  that: 

a)  for  smooth  rectangular  channels  of  infinite 

width   f  »  2^ /Re  (8ya) 
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b)  for  rough  triangular  channels  vith  90  vertex  angle 

f  »  60/Re  (83-b) 

(33) 
Woow^'  in  I956  found  tbat  a  linear  relation  exists  between  the 

friction  coefficient  f  and  the  Reynolds  number  Re  for  the  laminar 

flow  in  rough  rectangular  channel.   The  upper  critical  value  of  Re 

is  not  a  constant  but  ranges  from  *K)0  -  900.  The  test  data  also 

indicate  that  the  Reynolds  number  varies  inversely  with  the  slope  of 

the  plane,  and  the  roughness  seemed  to  have  no  consistent  effect  on 

the  upper  limit  of  laminar  flow.  Again  Woo  and  Braterv"   in  1962 

found  that  the  effect  of  the  rainfall  impact  could  be  interpreted 

as  a  change  in  the  roughness  coefficient.  Experiment  shoved  that  the 

relationship  between  the  friction  coefficient  f '  for  disturbed  flow 

under  a  particular  rainfall  intensity  and  Reynolds  number  could  be 

best  fitted  by  straight  lines  in  log- log  plot  parallel  to  the 

corresponding  line  for  undisturbed  flow  without  rainfall. 

For  turbulent  flow,  there  is  no  linear  relationship  between 

the  friction  coefficient  and  the  Reynolds  number.  The  tractive  force 

term  is  then  evaluated  by  using  empirical  formulas,  such  as  the 

Blasiue  equation 

0.223 


Re 
for  smooth  channel  or  Chezy'   or  Manning's  formulas,   which  give 

S,  -    tA  (95) 

1        C     R 

or 

sf  -    5^  (86) 

2.2  R  i* 
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Overland  flow  hydrograph 
The  overland  flow  produced  from  a  certain  precipitation  inflow 
can  be  divided  into  three  stages.  Each  one  has  its  own  individual 
characteristics.   This  can  be  expressed  with  the  help  of  overland 
hydrograph,  as  shown  in  Figure  (36).  The  three  stages  are: 

Stage  I.  This  is  the  rising  part  of  the  hydrograph  which 
starts  from  the  beginning  of  the  precipitation.  The  outflow  increases 
with  time  which  makes  the  flow  unsteady.  The  mathematical  formulas 
derived  in  the  previous  paragraph  are  especially  for  this  stage: 
Eqs.  (71)  and  (73)  are  the  equation  of  continuity  and  of  motion 
respectively 

dy     Su 


s*  +y  5x  +uig    =  i 


I?  ♦  u  |H  +  g  !*.   .  g(S  .  B  )  +  |E  „  £! 

dt     ox    B  dx     ov  o    f '        y     y 

where  S-  is  used  to  evaluate  tractive  forces  instead  of  t  . 
f  o 

Stage  II.   If  the  duration  of  precipitation  is  long  enough, 
the  overland  flow  can  reach  its  equilibrium  condition,  that  is, 
outflow  is  equal  inflow.  Thus,  the  outflow  and  the  profile  of  water 
surface  are  constant  irrespective  of  time,  that  is,  the  flow  is 
steady.  The  mathematical  formulas  for  the  overland  flow  in  this  stage 
can  be  expressed  as: 

*£  ♦»£  -1       -  fe  -»  ») 

»i  ♦«£  -*<so- V  r-%  (88) 
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This  is  shown  as  a  horizontal  line  in  the  overland  flow  hydrograph. 
The  equilibrium  outflow  is  obtained  by  integrating  Eq.  (87). 

**  °  To725o"  <89) 

where  L  is  the  length  in  feet  of  the  strip  of  land  subjected  to 
overland  flow  assuming  a  unit  width,  i  is  the  intensity  of  rainfall 
in  inches  per  hour,  and  ^3,200  is  a  conversion  factor  to  determine 
the  discharge  in  cubic  feet  per  second.   At  this  stage,  the  main 
point  of  interest  is  the  determination  of  the  water  surface  profile 
which  yields  the  maximum  detention  storage  over  the  ground  and  is 
an  important  information  for  airfield  and  highway  design.  The  profile 
equation  of  water  surface  of  overland  flow  in  equilibrium  stage  can 
be  derived  from  Eqs.  (87)  and  (88). 
By  substituting  Eq.  (87) 

du     _i     u   dy 
dx  *  y     y   dx 


into  Eq.  (88)  we  have 


g(So .  B , ,  m^ak 


§  -  * — - —  (90) 


or 

s(so  -  sf )  ♦  ^^ 
*SL    ,  (90-a) 

dx  2  2 

i  x 

y 

The  profile  can  be  obtained  by  numerical  integration  started  from  a 
known  control  section. 
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Stage  III.   When  the  precipitation  ceased,  the  overland  flow  reaches 
the  third  stage.  There  is  no  more  inflow  and  the  outflow  is  solely 
the  depletion  of  water  from  storage.   It  shows  as  a  recession  part 
in  the  overland  hydrograph.   Its  theoretical  formulas  are  of  the 
same  form  as  for  the  stage  I,  but  with  i  ■  0: 

H  ♦*£  ♦»&  -  •  «« 

&  ♦«!!  ♦«&  -  «<vsf>  <*> 

This  is  the  same  system  of  equations  as  that  for  unsteady  flow  in 

(32) 
open  channels    .   However,  the  boundary  conditions  at  the  upstream 

ends  are  different.   In  this  case,  the  detention  storage  remaining  on 
the  ground  right  after  the  end  of  the  precipitation  is  the  initial 
condition  and  the  boundary  condition  states  that  there  is  no  inflow  at 
the  upstream  end,  whereas  for  unsteady  flow  in  open  channels  the 
boundary  condition  is  usually  given  in  the  form  of  an  imposed  hydro- 
graph  at  the  upstream  end.  At  the  downstream  end,  in  both  cases,  the 
boundary  condition  is  usually  given  in  the  form  of  a  stage-discharge 
relationship. 

The  duration  of  precipitation  is  the  main  factor  to  develop  the 
How  in  Stage  II.   If  the  duration  of  precipitation  is  less  than  the 
time  of  concentration  of  the  watershed,  Stage  II  of  overland  flow 
can  never  be  obtained;  that  means  the  flow  can  never  reach  the 
equilibrium  condition. 
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Location  of  control  section 

The  control  section  is  the  section  where  the  critical  depth 
occurs,  the  flow  changes  its  condition  from  subcritical  to  super- 
critical or  vice  versa.  It  is  because  the  critical  depth  of  flow 
can  be  obtained  from  the  discharge  directly,  that  it  is  used  as  a 
boundary  condition  for  various  problems.  The  purpose  of  the  deter- 
mination of  control  section  is  thus  to  obtain  a  boundary  condition 
in  the  form  of  a  stage- discharge  relationship  for  integrating  the 
differential  equation  of  overland  flow. 

For  steady  uniform  flow  the  relation  between  discharge  and 


depth  at  critical  flow  is  q  «  I  gyf  where  q  is  the  discharge  per 

unit  width,  and  y  is  the  critical  depth.  For  overland  flow,  the 

C  x    <2  2 

5    i  X 

critical  depth  is  given  by  the  relation  y  =  where  i  is 

c     8 

the  intensity  of  precipitation  and  x  is  the  length  of  the  plane. 

It  has  been  assumed  by  many  persons  that  the  critical  depth 
occurs  at  the  end  of  the  plane  where  there  is  a  free  fall.  Tnis  is 
true  only  for  specially  varied  flow  when  the  slope  of  the  plane  is 
small  enough  that  the  subcritical  flow  occurs  throughout  the  whole 
plane.  The  situation  is  more  complex  for  the  overland  flow  produced 
by  precipitation,  since  the  discharge  is  a  function  of  the  length 
of  plane.  Thus,  the  location  of  the  critical  depth  point  is  not 
constant  over  the  plane  but  it  depends  upon  the  intensity  of  pre- 
cipitation and  the  length  of  plane. 

Li^  °'     derived  the  flow  profile  equation  from  the  equation 
of  spatially  varied  steady  flow  for  channel  with  parallel  walls. 
He  neglected  the  friction  loss  and  did  not  consider  the  change  of 


128 


momentum  due  to  rain  drops.  He  introduced  a  boundary  condition  ex- 
pressed  by  the  depth  yQ  and  the  Froude  number  Fq   at  the  end  of  the 
channel  of  length  xq.   The  differential  equation  of  the  water 
surface  profile  is  thus: 

n»  ( h  >  -  <  f  >5 ) 4  (  V )  - 

°       L  yo  yo 

-U*!  (JL.  )(»).0(JL   )5]d|  (95) 


yo         L  yo 


where 


and 


F' 


2  2 

2  IT 

o  =     ~ 
^o 


G     = 


S  L 
o 

yo 


A  simple  expression  can  be  derived  for  horizontal  channel,  in  which 
case  the  above  differential  equation  may  be  integrated  to  obtain 

<I  >2'  <1+  -^X-v"  >  ~h   (  y"  )?  (?U) 

2F     yo    2F    yo 
o  o 

This  shows  that  if  the  end  condition  is  a  free  fall,  that  is  F   -  1 

'  o 

the  water  surface  profile  can  be  determined  from 

(  *  )2  „  2      I  *  (  -JL  )5  .  (95) 

*  L  ;      2   y     2  *  yQ  ;  v^' 

But,  for  the  channel  with  sloping  bed,  a  general  explicit  equation 
of  the  flow  profile  as  Eq.  (9*0  can  not  be  obtained.   However,  it 
may  be  solved  by  numerical  integration  starting  from  a  known  control 
depth  of  flow.  A  general  diagram,  shown  in  Fig.  (57),  indicates 
the  relationship  between  the  Froude  number  F  and  G.   It  consists  of 
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G  =  2 
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four  zones  representing  the  different  flow  conditions,  according  to 
Li: 

Zone  A  represents  the  condition  where  flow  is  sub- 
critical  throughout  the  channel.  The  Froude  number  F 

increases  as  the  flow  proceeds  downstream,  that  is.  -7-  >  0. 

'  '   dx 

Zone  B  represents  the  condition  where  the  flow  is 
subcritical  throughout  the  channel  but  where  the  value  of 
F  will  first  increase  as  the  flow  proceeds  downstream, 
reaching  a  maximum  value  less  than  unity  and  then  decrease. 

Zone  C  represents  the  condition  where  the  flow  is 
supercritical  in  the  downstream  portion  of  the  channel  and 
there  is  a  hydraulic  jump  in  the  channel. 

Zone  D  represents  the  condition  at  which  there  is 
supercritical  flow  throughout  the  downstream  portion  and 
no  hydraulic  Jump. 

In  the  physical  phenomena  of  overland  flow,  a  free  fall  end 
condition  where  the  water  flows  from  the  overland  flow  to  the  channel 
is  usually  assumed  and  the  water  surface  in  the  channel  is  supposed 
to  be  low  enough  so  as  not  to  affect  the  end  condition  of  tne  over- 
land flow.  Three  types  of  overland  flow  may  therefore  exist; 

1.  The  flow  is  subcritical  throughout  the  channel  and  a 
critical  depth  occurs  at  the  end  where  there  is  a  free 
fall. 

2.  The  flow  is  supercritical  throughout  the  channel  when  the 


131 

bottom  slope  is  big  enough. 
3.  When  the  slope  of  the  bottom  plane  is  not  large  enough, 
critical  flow  occurs  somewhere  on  the  plane  dividing  the 
flow  into  a  subcritical  region  and  a  supercritical  region. 
There  will  be  no  difficulty  to  find  the  control  section  for 
the  first  and  second  types  of  overland  flow  as  it  is  easily  seen 
that  the  control  section  for  the  first  type  is  at  the  end  of  the 
plane  and  for  the  second  type  is  at  the  starting  point  of  the 
plane  where  the  critical  depth  y  is  zero.  For  the  third  type  where 
both  subcritical  and  supercritical  flows  occur  on  the  plane,  it  is 
necessary  to  find  the  point  where  the  critical  depth  occurs. 

The  position  of  a  critical  section  can  be  determined  by  the 

(kr) 

method  of  singular  point  as  was  clone  by  Keulegan     for  the  spatially 

variable  flow.  Both  the  numerator  and  the  denominator  of  the  dif- 
ferential equation  of  the  free  surface,  Eq.  (90)  are  made  equal  to 
zero.   When  the  denominator  is  made  equal  to  zero,  the  critical  flow 
profile  is  obtained: 

,     i2x2 

When  the  numerator  is  made  equal  to  zero,  the  relationship 

2 
S  -  S_  -  Q*     ,  0    (v  is  neglected)  (97) 

OI       c 

gy 

yields  the  quasi-normal  flow  profile.  The  control  point  is  located 
at  the  intersection  of  the  critical  flow  profile  and  the  quasi- 
normal  flow  profile  where  x  ■  x  and  y  «  yc 

<so>c-   <V=-  ^r  -°  (98) 

gyc 
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By  using  Manning's  formula  to  evaluate  the  term  S   the  previous 

equation  becomes ,  remembering  that  u  =,  Jgy 

c      c  y 


2  2 

gn  2i  xc 


(so),  ■  rnr   -      — er-  -  °  (99) 


2.2  y/        g  yc 


Since  the  bottom  slope  at  the  control  section  is  equal  to  the 
critical  slope,  that  is  S  »  S  ,  and  substituting  the  relation 
of  Eq.  (96),  we  have 


gn  2y 

S  -   =—   -   — £-  »  0  .  (100) 

2.2  yV3 


The  Eqs.  (96)  and  (100)  are  two  algebraic  equations  with  two  un- 
knowns y  and  x  .  The  position  of  the  control  section  can  thus  be 
c      c 

determined.  The  slope  of  water  surface  at  the  singular  point  is 
obtained  by  use  of  the  L' Hospital's  rule  to  take  the  limit  of  the 

value  dy/dx  at  x  ■  x  . 

(1*2) 
There  is  another  method  suggested  by  Hinds      for  the  de- 
termination of  the  control  section.  It  is  to  find  the  total  head 
loss  due  to  friction  and  the  impact  from  rain  drops  at  every  sec- 
tion of  which  the  plane  is  assumed  to  have  critical  flow  everywhere. 
Therefore,  a  water  surface  profile  and  a  bottom  profile  of  criti- 
cal flow  can  be  determined. 
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The  control  section  is  located  at  the  point  where  the  slope  of  bottom 
profile  of  critical  flow  is  equal  to  the  actual  bottom  slope  of  the 
plane.  The  drop  of  head  due  to  friction  is  evaluated  by  Manning  formula. 

We  have 

2  2 
n  u 


b'"  *****  7T77T  «■  (10l) 


2.2  y 

Since  it  is  critical  flow  everywhere,  we  can  write 

2 

h  «  S  Ax  -  n  f .    Ax  .  (102) 

2.2  y'0 
c 

The  drop  of  head  due  to  impact  loss  caused  by  the  precipitation  in- 
flow can  be  evaluated  from  the  equation  of  motion.  If  we  write  the 
equation  of  motion  in  finite  increment  form,  equ.  (88)  becomes 

u  Ax    B  Ax    •  o   *  t         y 
the  total  loss  of  head  can  be  expressed  as 

Ay  -  -  (  —^     +  —  )Ax  ♦  Sn  Ax  -  S_  Ax  (103) 

^       g  Ax      gy  *^         0       f 

the  drop  of  head  due  to  impact  is  Ay* 

Ay.  .  *    Au  +  AH  ^  (10U) 

g     gy 

by  introducing  the  relation  ix  •  uy,  oridx-udy+ydu 

Ay'   -    —    Au  +    —    Ay  .  (105) 

g  gy 

But,  by  considering  the  physical  phenomena  of  the  overland  flow,  the 
loss  of  head  due  to  the  impact  is  relatively  much  smaller  than  that 
due  to  the  friction  since  the  depth  and  velocity  are  small  (for 
example,  for  the  precipitation  intensity  2in./hr.,  the  inflow 
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i  »  0.00005  cf s /square  feet,  if  the  length  L  -  1,000  feet,  the  critical 
depth  yc  is  shout  0.01— 0.0^  feet,  and  for  n  -  0.015  the  velocity  is 
ahout  0.5— 1.2  feet  per  second).  An  approximate  solution  can  he  ob- 
tained hy  using  only  the  friction  term  Sf  and  neglecting  the  loss 
due  to  impact.  This  gives  a  very  simple  method  to  find  the  location 
of  the  control  section.  At  the  control  section 

2 

Sc  *    *  1/3  (106) 

2.2  y  X'? 

The  same  relation  can  he  obtained  from  Eq.  (100)  vhich  is  derived 

by  the  method  of  singular  point.   If  the  control  section  is  not  close 

to  the  starting  point  of  the  plane,  the  value  of  y  is  much  smaller 

than  that  of  x  .  Therefore  the  last  term  in  Eq.  (100)  can  be  dropped 


and  ve  have  Eq. 

(106) 

S     - 
c 

2 

2.2  yV5 
c 

Since  S     is  equal  to  S    at  tl 
c                            o 

by  Eqs.    (96)  and 

(106) 

S     = 
0 

2 
gn 

2  2     V5 

d'd  yc 

v5. 

42  2 
ixc 

This  simple  method  for  the  determination  of  the  control  section  is 
valid  only  for  the  overland  flow  produced  by  precipitation.  This  is 
not  valid  for  other  cases  such  as  lateral  spillway  channel  where  the 
depth  and  the  velocity  of  flow  are  large.   In  that  case,  the  loss  of 
head  due  to  the  impact  is  relatively  much  larger  than  that  due  to 
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friction. 

The  approximate  determination  of  the  control  point  as  given 
above  is  only  for  the  steady  case.  It  is  not  valid  where  the  flow 
is  unsteady.  The  critical  depth  will  shift  upstream  when  the  flow 
is  unsteady.  Since  the  unsteady  flow  profile  changes  with  time,  the 
location  of  critical  depth  is  also  a  function  of  time. 

Initial  and  boundary  conditions 
As  indicated  in  the  previous  paragraph  of  runoff  phenomena 
on  a  watershed,  there  is  an  initial  detention  and  there  is  an  initial 
time  lag  before  the  runoff  starts.  Theoretically  and  physically  these 
are  the  initial  conditions  for  the  overland  runoff,  but  no  definite 
experimental  values  are  available.   Therefore,  for  the  initial  con- 
dition it  is  assumed  that  the  mean  velocity  u  and  the  depth  y  of 
overland  flow  are  zero  when  time  t  is  zero,  this  may  be  true  especially 
for  small  areas  where  the  actual  initial  detention  and  initial  time 
lag  are  very  small. 

For  the  boundary  conditions,  it  is  assumed  that  there  is  no 
inflow  at  the  upstream  end, that  is,u  and  y  are  both  zero  for  x«0. 
For  the  downstream  end  conditions,  a  free  fall  is  assumed.   It  forms 
a  control  point  if  the  slope  of  the  plane  is  small  enough,  and  gives 
a  boundary  condition  u«c  at  x*L.   If  supercritical  flow  occurs 
at  the  downstream  end,  no  boundary  condition  is  needed,  as  the  char- 
acteristic directions  are  in  the  same  direction  in  this  case. 
The  initial  and  boundary  conditions  are  as  follows: 
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1.  t  -  0,       u(x)  »  0     y(x)  «  0 

2.  x  «  0,        u  »  0       y  =  0 

5a)  x  -  L,                   u  =  >/gy  «  c    subcritical  flow  and  free  fall. 

5b)  *  "  L,  no  boundary  condition  needed  for  supercritical  flow. 
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METHOD  OF  CHARACTERISTICS 
The  method  of  characteristics  is  a  mathematical  tool  commonly  used 
to  find  the  numerical  solution  of  systems  of  partial  differential  equa- 
tions of  the  hyperbolic  type.  This  study  considers  primarily  the  case  of 
two  simultaneous  quasi- linear  partial  differential  equations  with  two 
dependent  and  two  independent  variables.  The  material  in  this  chapter 
is  based  on  the  notes  taken  by  Prof.  J.  W.  Delleur  in  a  course  on  "Methods 
of  Characteristics"  taught  by  Professor  Julien  Kravtchenko  in  1961  at  the 
Fluid  Mechanics  Laboratories,  University  of  Grenoble,  in  Grenoble,  France. 

Quasi- Linear  Differential  Equation 
The  general  form  of  a  quasi- linear  system  of  equations  for  the  case 
of  two  independent  variables  x,  t  and  two  dependent  variables  u,  y  can  be 
written  in  the  form: 

where  A,  A   B,  B1,  ,  E,  E1  are  known  functions  of  x,t.   Since 

these  coefficients  are  figuring  linearly,  they  are  called  "quasi- linear" . 

The  Cauchy  Problem 
The  Cauchy  problem  is  that  of  solving  equations  (  107)  and  (108) 
from  prescribed  initial  and  boundary  conditions  in  an  open  domain, 
as  illustrated  in  Fig.  (  58  ) ■  The  generalized  initial  and  boundary 
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conditions  can  be  expressed  by  the  "€auchy  data"  which  specify  along 
an  arbitrary  curve  in  the  plane  of  the  independent  variables 

x  "  x(t)  (109) 

the  value  of  the  dependent  variables 

u  «  u(t)  (no) 

y  -  y(t)  (111) 

The  Cauchy  data  may  be  represented  graphically  by  tracing  through  the 
curve  x  ■  x(t)  in  the  x-t  plane,  a  cylindrical  surface  in  the  z  direc- 
tion on  which  the  curves  C  and  C  represent  the  given  conditions  of 
Eqs.(109),(H0),(iii)as  shown  in  Fig. (59). 

The  solution  of  the  system  of  quasi- linear  partial  differential 
equations  (10?)  and  (108)  satisfying  the  conditions  of  Eq. (109)- (111)  is 
represented  graphically  by  two  surface  u(x,t)  and  y(x,t)  passing  through 

the  curves  C  and  C  . 
u      y 

Taylor  Series  Expansion  of  the  Solution 

Suppose  that  the  functions:  A,  B,  C,  ,   A   B  ,  E  as 

well  as  u  and  y  are  analytic.  The  functions  u  and  y  may  be  represented 

then  by  a  Taylor's  series  expended  in  the  neighborhood  of  the  known 

point  x  ,  t  ,  z 
r  o   o'  o 

u(x,t)  .  u  (xot0)  ♦  (x  -  x0)(  §2  )o  ♦  (t  -  to)  ($j  )o  ...    (lis) 

y(x,t)  -  y  (xoto)  ♦  (x  -  xo)(  |  )0  *  (t  -  tQ)  (g  )0  ...    (ii„ 

where  the  high  order  terms  are  neglected.  In  order  to  evaluate  u  and  y, 
it  is  necessary  to  find  the  four  unknown  partial  derivatives  (  -^  )q, 


^  at  V  *  dx   ^o'    (  dt   V    As  we  ****  only  tw0  equations  for  which  the 
coefficients  are  known  at  xq  and  tQ  (Eqs.(l07)  and  (108  ),  it  is  necessary 
to  obtain  two  supplementary  equations  in  order  to  get  a  system  of  four 
equations  for  the  four  unknown  derivatives.   If  the  u«u(x,t)  and  y«y(x,t) 
are  solutions  to  the  quasi-linear  system  Eqs.(l07)  and  (108  ),  then 
the  two  supplementary  equations  can  be  written  as: 

*  -  I to  *  $  «*  ui5) 

Therefore,  there  are  four  equations  (107),  (108),  (H1*),  (H5)  witn  f°ur 
unknowns  which  are  the  four  partial  derivatives  (  -<r-  )  ,  (  7S7  )  >  (  X"  )  » 
and  (  -r^  V,  They  may  be  obtained  by  solving  these  four  simultaneous 
algebraic  equations,  and  hence  the  value  of  u  and  y  are  obtained  in  the 
neighborhoood  of  x  ,  t  ,  and  z  . 

Theory  of  Linear  Algebraic  Equations 
The  necessity  of  solving  simultaneous  equations  (107),(108)(llU), (115) 
for  the  derivatives  suggests  a  review  of  the  theory  of  linear  algebraic 
equations.   First,  let's  consider  a  system  of  n  linear  equations  with  n 
unknowns.  They  are 
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Let  D  denote  the  determinant  of  the  coefficient  of  the  n  unknowns 
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then  the  solution  of  n  linear  equations  expressed  by  determinant  are; 
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A  general  expression  can  be  written  as 

^i  "  Ki  (119) 

By  studying  Eq.  (119), one  can  find  that  there  are  four  different  cases 
of  solutions: 

1.  D  #>  and  K  ^  0 

The  unique  values  of  x,,  x_,  ,  x  can  be  obtained. 

j.   c        n 

2.  D  }   0  and  K  «  0 

All  n  unknowns  are  zero,  x  =  0 

3.  D  *  0  and  K  +   0 

There  is  no  solution,  since  x  *  « 
k.     D  =  0  and  K  =  0 

The  solution  is  indeterminate,  there  are  sets  of  solutions. 

By  studying  the  above  four  different  cases  of  solutions,  it  is 
obvious  that  the  cases  2  and  3  are  impossible  and  there  will  be  either 
a  unique  solution  as  given  by  case  1  or  sets  of  solutions  as  given  by 
case  k.     As  for  the  unsteady  flow  problem,  the  coefficients  of  the  quasi- 
linear  system  of  partial  differential  equations  still  contain  unknowns 
and  this  makes  it  impossible  to  evaluate  the  determinant  as  a  unique 
value.  The  only  way  to  avoid  the  mathematical  difficulty  is  to  set  the 
determinant  equal  to  zero  and  obtain  the  solution  as  given  by  case  U. 
This  also  can  be  explained  by  considering  that  a  flood  wave  is  composed 
of  a  great  number  of  infinitesimal  surges,  and  each  surge  has  a  discon- 
tinuous surface  profile  which  mathematically  is  t^  *  q  • 


1^3 


Classification  of  Systems  as 
Hyperbolic .  Parabolic  or  Elliptic 

Since  we  have  the  quasi- linear  system  of  two  equations  and  two  sup- 
plementary equations  (Eqs. (107), (108), (llU)  and  (115),  it  can  be  written 
as  a  simple  matrix  equation 
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If  the  matrix  vanishes,  that  is  D  ■  0,  we  have 

(AC1  -  A  C)dx2  -  (ADX  -  AD  +  BC  -  B^Cjdt  dx  +  (BD-  -  B.D)dt2-o  (121) 

or 

(A^  -  AXC)  (  |2  f   -  (AD1  -  A2D  +  BCX  -  B^C)  ^  +  (BD1  -  B^)  *0  (122) 

This  may  be  considered  as  a  quadratic  equation  for  the  slope  dx/dt. 

The  quadratic  equation  (  122)  gives  two  real  slopes,  one  real  slope,  or  a 

pair  of  complex  slopes  depending  on  whether  the  discriminant 


6  -  (ADX  -  AD  ♦  BCX  -  B-jC)2  -  *  (AC^  -  A^Jf^  -  B^D) 


(123) 


dx 


is  positive,  zero,  or  negative.  The  directions  given  by  ^    which  satisfy 

Eq.  (  122)  are  called  characteristic  directions.   If  the  discriminant  is 

such  that 

5  >  0  the  system  is  called  Hyperbolic 
5  «  0  the  system  is  called  Parabolic 
B  <  0  the  system  is  called  Elliptic 


Ikk 


For  the  hyperbolic  type  there  are  two  real  characteristic  directions,  and 
the  solutions   may  be  obtained  by  the  method  of  characteristics. 

Characteristic  Equations 
The  quasi- linear  system  of  equations  of  the  hyperbolic  type  may  be 
solved  by  the  method  of  characteristics.  Eq.  (123)   has  two  real  roots 
for  dx/dt  which  give  two  characteristic  directions  in  the  x-t  plane.  They 
can  be  expressed  as 

xl  "  f  (t,*,u,y)  (12U) 

x2  =  f  (t,x,u,y)  (125) 

which  are  the  first  two  characteristic  equations,  or  the  slopes  of  the 
characteristic  curves  of  the  first  and  second  family  respectively. 

By  substituting  the  two  supplementary  equations  (ll1*),  (115)  into 
the  two  quasi- linear  equations  (107)  ,    (108)  ,  we  have 

(B-Ax-jg  ♦(D-CMfjj-l-All-cJI  (126) 

<Bi  -  V'»  5  +  (Dr  cix,)  I '  Ei  -  Alli  "  °1  i  W 

The  system  of  equations    (126)  ,    (127)     is  a  system  of  two  equations  with 

two  unknown  du  and  dy  ,        if  the  solution  is  indeterminate,   it  is  necessary 

cix  ^x 

to  have  the  following  relation 

E      a    iS  .  c  & 
B  -  Ax'  D  -   Cx'  dt  dt 

B.    -  A.x'      °     D.    -   C.x'      "  .  Hv  (12P,)> 

11  11  ra^.p^ 

Bl      Al  dt       Cl  dt 


or 


(B-Ax.XE,  -  *!  g  -  Cx  g)  -   (S,  -  A^'XE  -  *  $  -  C  $)  -  0  (1;v,; 
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which  can  be  expressed  as 

(AB1-BA1)du  +  [(CB^B)  ♦  (AC^Cjx'  ]dy  ♦  (BE^Ejdt  ♦  (EAj-AE^dx  -  0. 

(130) 
Therefore,  we  have  two  other  characteristic  equations 

(AB-pAjBjdu  +  [(CB-j-K^H  (AC^  -  A^Jx^dy  +  (BK1-EB1)dt  +  (A1E-AE1)dx=  0 

(131) 

(AB^A^Jdu  +  [(BXC  -  BCX)  +  (AC^A^Jx^ldy  +  (BE^B^dt  *   (A-jE-AE^dx  =  0 

(132) 
The  equations  (124)  ,   (125),  ( 131)  and  (132)  are  four  characteristic 

equations.  The  first  two  characteristic  equations  give  the  directions 

of  the  characteristic  curves;  the  last  two  characteristic  equations  must 

be  satisfied  along  the  characteristic  curves. 

Method  of  Calculation 
A  method  of  approximate  integration  is  used  to  find  the  solution 
of  the  quasi- linear  system  of  equations  (I07)and  (lo8)with  the  given 
Cauchy  data  as  shown  in  equations  (109),  (110)  and  (ill).  The  problem 
can  be  solved  semigraphically  (as  shown  in  Fig.  ^O  ),  if  the  two  charac- 
teristic directions  are  known.  For  a  known  point  M(x  ,t  )  on  the  given 
curve  x  «  x(t),  the  corresponding  values  of  u  ,y  are  known,  we  therefore 
know  the  direction  of  tangents  to  the  characteristic  curves  at  this 
point.  The  tangent  to  the  characteristic  curve  of  the  first  family  may 

be  traced  at  x  ,y  .   The  same  procedures  may  be  done  for  the  point 
o'  o 

M'(x.,t.  )  where  the  tangent  to  the  characteristic  curve  of  the  second 
family  is  traced.  Since  the  directions  of  the  characteristic  curves 
are  different  they  intersect  at  a  point  N  at  which  the  values  of  u  and  y 
can  be  determined  as  will  be  shown  later. 

Instead  of  the  three-dimensional  representation  of  Fig.  *+0 


lk6 

it  is  convenient  to  consider  two  auxiliary  planes:  the  x-t  plane  of  the 
independent  variables  and  the  u-y  plane  of  the  dependent  variables.  The 
points  in  the  x-t  plane  map  in  a  univocal  or  one-to-one  correspondence 
onto  the  u-y  plane.  The  following  mathematical  expressions  will  show 
how  to  compute  the  value  of  u  and  y  corresponding  to  the  intersection 
point  N  on  the  x-t  plane.   As  we  have  the  four  characteristic  equa- 
tions (124),  (125),  (131),  and  (132),  if  we  let 


AB1  "  »l 


b  -  CB1  -  C^B 


AC1  "  <*1 


BE  -  EB 


e  «  EA  -  AE1 
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and  express  the  differential  form  to  a  small  increment  form 

du  *  u  -  u      dy  =  y  -  y     dt  »  t  -  t  (13^) 

p  p  n 

then  the  equations  (  131),  ( 132)  can  be  written  as 
a(M,p)(u  -  u  )  +  [b(M,p)  +  x^(M,p)c(M,p)](y  -  yp) 
+  [d(M,p)  +  e(M,p)  x^(M,p)](t  -  t^)   -  0  (13J) 

a(M»,p')(u  -  u  ,)  +  [b(M',p«)  +  X2(M',p')C(M',p')](y  -  ypl ) 

+  [d(M\p')  +  e(M',p')x2(M',p')](t  -  t^, )  -  0  (156) 

therefore,  there  are  two  equations  for  determining  u,  y  and  t.   But, 
since  the  two  characteristics  intersect,  the  point  N  found  in  the  x-t 
plane  can  be  used  to  determine  t,  that  is  t  -  tN>   In  the  u-y  plane, 
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SEMI-GRAPHICAL       SOLUTION       OF 
THE      CAUCHY        PROBLEM 
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FIG.  41    -  ZONE       OF       INFLUENCE       FOR        GIVEN 


INITIAL        CONDITIONS        ALONG        AB 
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Q  is  the  image  of  the  point  H,  we  can  therefore  write  vu,y  instead 
of  u  and  y  in  Eqs.  (155)  and  (156).  There  are  two  equations  with  two 
unknowns  uQ  and  yQ  from  which  the  point  Q  is  completely  determined. 

The  above  shows  the  procedures  to  solve  the  characteristic 
equations.  This  method  can  be  carried  on,  point  after  point,  to  obtain 
the  total  solution  of  the  zone  confined  by  the  two  general  characteristic 
directions.  This  is  shown  in  Fig  s.  1*0, 'U,  the  zone  ABC  which  indicates 
all  the  solutions  for  the  given  boundary  condition  will  not  be  affected 
by  the  other  boundary  conditions.  This  can  be  called  as  a  zone  of  in- 
fluence. The  solutions  outside  this  zone  can  be  obtained  by  some  other 
boundary  condition. 

Special  Problems  of  Cauchy  type  I,   II,  III 
Type  I 

The  given  condition  x=x(t)  coincides  with  one  of  the  characteris- 
tics. As  shown  in  Fig.  (^2  ),  if  the  characteristic  of  the  first  family 
coincides  with  the  given  curve  x»x(t),  the  two  characteristics  intersect 
at  point  M'  which  is  already  known,  therefore  a  corresponding  point  Q 
in  the  u-y  plane  can  not  be  determined  uniquely.  A  problem  of  this  type 
has  no  solution. 

A  mathematical  expression  can  be  used  to  interpret  this  type  of 
problem.  The  problem  is  unsolvable  if  the  Jacobian  or  the  functional 
determinant  equals  zero. 


D(x,t) 
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Fig.  W2.  Special  problem  of  Cauchy  type  I. 


Fig.  *0.   Special  problem  of  Cauchy  type  II, 


Xi  r  X|(t) 


Fig.  Uk.     Special  problem  of  Cauchy  type  III. 
C  =  characteristics  of  the  first  family 

C  «  characteristics  of  the  second  family 
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Type  II 
There  are  two  sets  of  given  conditions 


X±   -  xx(t)         U]L  «  ux(t)        yx  -  yx(t) 
X2  -  x2(t)         u2  -  u2(t)        y2  -  y  (t) 


(158) 


where  the  first  coincides  with  the  first  characteristic  and  the  second 
coincides  with  the  second  characteristic.  Suppose  that  these  two  sets 
of  given  conditions  have  a  common  point  which  is  known.  The  graphic 
representation  can  be  shown  in  Fig.  (^5  ).   In  x-t  plane,  two  charac- 
teristics can  be  traced  from  the  known  points  M  and  M'  and  then  N  can  be 
determined.  The  point  Q  in  u-y  plane  is  determined  from  Eqs.  (155)  and  (156). 
The  total  solution  and  the  zone  of  influence  is  shown  in  x-t  plane. 

Type  III 
If  the  given  data 

x  »  x(t)     u  «  u(t)     y  ■  y(t) 
represents  a  characteristic,  these  data  are  not  sufficient  to  obtain  the 
solution.  One  more  condition  which  is  not  a  characteristic  is  given  on 
the  x-t  plane,  that  is 

xl  *  xi^^ 
or 

F(t,u,y)  -  0  (159) 

but  neither  u(t)  nor  y(t)  are  given,  therefore  their  image  in  u-y  plane 
is  not  known.  The  solution  of  this  kind  can  be  obtained  as  shown  in  Fig.  ^ 
The  point  Mf  which  is  the  intersection  of  the  given  curve  and  another 
characteristic  on  the  x-t  plane  is  given  by  Eq.  (155) 
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a<UQ  "  Up>  +  (t>  +  CX')(yQ  ■  V  +  (d  +  "'MV  "  tM)  "  ° 
and  the  given  condition 


or 


*<V-  v  V  -  (  f  )Mp  <v  -  v +  <  I V  <\  -  •>> 

These  are  two  equations  for  the  determination  of  two  unknowns  xu   and  y  . 
The  total  solution  can  he  obtained  by  extending  the  characteristics  on 
x-t  plane. 

Characteristic  Variable 
As  shown  in  preceding  sections,  we  have  a  hyperbolic  system  of 
two  equations  with  two  unknowns  and  two  independent  variables.  The 
solution  may  have  the  form 

u  =  u(x,t) 

y  -  y(x,t) 

The  method  of  characteristics  gives  two  families  of  characteristic 
direction  as  shown  by  Eqs.  (12h)   and  (125),  and  the  Eqs.  (131),  (152) 
must  be  satisfied  along  the  first  and  second  families  of  characteristic 
curves  respectively.   The  surfaces  u  and  y  are  thus  generated  by  a 
double  network  of  characteristic  curves  each  of  which  depends  on  an 
arbitrary  constant,  a  for  the  first  family  and  p  for  the  second  family. 

Thus  the  two  families  of  characteristic  curves  on  the  x-t  plane  can  be 

expressed  as: 
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xx  -  xx(t,  a)         x2  -  x2(t,  p)  (lU) 

and  the  corresponding  u,y  are 

uL  -  u^t,  a)  u2-  u2(t,  p)  (ite) 

yx  «  yx(t,  a)         y2«  y2(t,  p)  (1*3) 

For  any  point  N  which  is  the  intersection  of  two  characteristics 

x  -  xx(t,  a)  -  x2(t,  p)  (Ikk) 

we  can  therefore  take  a  and  p  as  independent  variables 

x  -  x(a,  p)  (1^5) 

t  -  t(a,  p)  U1*6) 

and  also 

u  «  u(x,t)  »  u(a,  p)  W) 

y  -  y(x,t)  »  y(a,  P)  U1*) 

By  taking  the  characteristic  variables  a,  p  as  two  new  Independent 
variables,  the  complete  solution,  i.e.,  u  and  y  can  be  evaluated.   In 
some  cases  a  physical  interpretation  of  the  solution  may  be  more  readily 
obtained  in  terms  of  the  characteristic  variables  a  and  p  than  in 
terms  of  the  independent  variables  x  and  t. 

Thus,  one  can  transform  the  original  hyperbolic  system  of  equations 
into  a  new  one  by  using  the  two  characteristic  variables  a,  p.   Let's 
consider  a  small  displacement  along  the  characteristic  a  -  constant, 
that  is  oaly  p  varies,  then 

dt  «   —  dp  du  -  jn   ^ 

dt     d  P  ^  3  P  (1U9) 

dx  -  M-   dp  dy  -  §L  dp 
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From  characteristic  equation  (124)  as  p  only  varies,  it  follows  that 

dx    w        m.  \     dt 

dp  •  fiu>  y'  x'  *>  dp  (150) 

Along  the  characteristic  p«constant,   only  a  varies,  and  from  the  charac- 
teristic equation  (125),   it  follows  that 

||  -  f2(u,  y,  x,  t)    g  (151) 

The  two  equations  ( 131)  and  ( 152)  which  hold  along  the  two  characteris- 
tics become: 


■I?  +(b  +  cv!  +<d  +  efi>!  -°  (152) 

'  35    +  (b  +  Cf2»    fl      +  (d  +  *f2>    I    "  0  (153) 

The  equations  (150),  (151),  (152)  and  (153)  form  a  system  of 
four  equations  with  four  unknowns,  5-   ,  St  ,  §•  and  ^  .  All 
solutions  for  this  system  are  the  solutions  for  the  original  hyperbolic 
system. 

Jacob lan,  and  the  Theory  of  Singularity 
Since  the  solution  of  a  hyperbolic  system  can  be  evaluated  in 
terms  of  the  two  characteristic  variables  a,  P   we  have  the  relations: 
x  «  x(a,  p) 
t  -  t(a,  p) 
The  above  transformations  are  possible,  and  there  is  a  one-to-one 
correspondence  between  a,   P,  and  x,t  if  the  Jacobian  or  the  functional 
determinant 
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The  necessary  and  sufficient  condition  that  x  and  t  can  be  expressed  as 
functions  of  a,  0  is  that  the  Jacobian  does  not  vanish. 
In  case  that  the  Jacobian  vanishes,  we  have 

J  =  0 
or 


dt   '   dt   ' 

5a   2   5p  1 

/  '    '  v  dt 

dt 

"  (x2  "  xi  >  55 

5e 

dt      dt 

5a     5p 

U5M 


This  is  to  say  that  J*0  when  x'  ■  x'  or  £_  ■  0,  °-r  «  0.  If  x'  »  x*, 
that  means  there  is  only  one  characteristic  direction  instead  of  two 
in  the  x-t  plane.  The  problem  is  no  longer  of  the  hyperbolic  type, 
but  of  the  parabolic  type,  the  method  of  characteristics  cannot  be 
applied . 

Consider  now  the  case  |-  »  0,  or  £g  ■  0  and  x^,  x£  are  not 


zero.  Since  we  have 


55  x2 


dx 

55 


at   ■       ax 

5p  xi  "  5p 


then 


f-  =  0  also  requires  that  ^-  •  0 
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and  <**  _  o  also  reauires  that  ^  »  0 


2|  *  0  also  requires  that  gj 
As  we  know  that  x  -  x(t)  is  a  curve  in  the  x-t  plane,  the  slope  of  the 
curve  can  he  expressed  as 


dx 
tan  9  =   -^ 


if  we  write 

dx/da 
tan  e  "  at/da 

then 

0  (155) 

tan  9  ■   -^  • 

This  forms  a  singular  point,and  the  solution  is  undetermined.    This 

shows  that  when  Jacohian  is  equal  to  zero,  the  two  characteristics 

directions  reduce  to  only  one  characteristic  direction,  or  as  55  -  0, 

|£  .  0,  singularities  exist  on  the  x-t  plane. 
^    Now,  let  us  consider  the  envelope  of  a  family  of  plane  curves 

C  on  the  x-y  plane.  The  equation  of  a  family  of  plane  curves  C  is 

given  hy 

(Kx,  y,  x)-o  or  sr-fCx,  O  (156) 

where  X  is  an  arbitrary  parameter  which  changes  the  for.  and  position 
of  the  curve  C  and  plays  the  sa*e  role  as  the  a    or  p  for  the  charec- 
teristlc  curves.   if  each  of  the  positions  of  the  fs*ily  of  curves  C 
ere  tangent  to  a  fixed  curve  E,  the  curve  E  is  ceiled  the  envelope 

n       c^^  the  envelope  c\irve  is  no  longer 
of  the  given  family  of  curves  C.   Since  the  env  j.  p» 

a  function  of  X,  we  have 

d  x 
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therefore,  the  envelope  curve  E  can  be  determined  by  eliminating  the 
parameter  \  between  equations  (156)  and  (157),  and  found  to  be  in 
function  of  x  and  y  only, 

R(x,y)  -  0  (158) 

If  we  consider  a  certain  point  MQ  which  is  x  ,  y  as  the  intersection 

point  of  the  two  curves;  the  envelope  curve  E,  and  the  curve  C  which 

is  a  particular  curve  corresponding  to  a  value  x  of  the  parameter. 
For  the  envelope  curve,  we  have 

S  <  s ».  +  %  <  IK  ♦  <  8  >8  -  •  •  <«»> 

Since 

M  (  p  )  ♦  M  A  )  .  0  (160) 

dx  v  dx  'o    dy  vd\  'o  i*w/ 

For  the  curve  C  ,  since  \  in  equation  (157)  has  a  constant  value  for 
the  particular  curve,  we  shall  have 

g*x*  |«*  -0  (161) 

By  comparing  the  equations  (160)  and  (l6l),  it  shows  that  the  tangent 
to  the  curve  C  coincides  with  the  tangent  to  the  curve  dexcrlbed  by 
the  point  (x,y),  at  least  unless  |£  and  ^2   are  both  zero,  that 
is  the  point  M  is  a  singular  point  for  the  curve  Cq.   It  follows 
that  the  equation  R(xy)  ■  0  represents  either  the  envelope  of  the 
curves  C  or  else  the  locus  of  singular  points  on  these  curves. 

As  mentioned  before,  there  will  be  no  one-to-one  correspondence 
between  a,  3,  and  x,y  when  the  Jacobian  vanishes.  The  Jacobian 


D(a,p) 


-  0 
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vhen  the  derivatives 


dor.  g  ta 


are  linearly  dependent.  In  this  case  the  function  of  x,y  can  be 
expressed  as 

0(x,y)  -  0  (162) 

Since  the  singularity  exists  vhen  the  Jacoblan  vanishes,  0(xy).  0 

could  be  the  locus  of  singular  points  and  also  the  envelope  of  the  families 

of  characteristics. 

When  the  Jacobian  is  equal  to  zero,  the  results  of  the  given 
hyperbolic  system  can  also  be  shown  graphically.  Let  us  consider 
the  representation  of  a  surface  in  space 

x  «  x  (a,  p)        y  »  y  (a,  p)        z  -  z  (a,p) 
The  tangent  plane  to  this  surface  has  the  parametric  derivatives 


3a 

^x 


OY 
5p 


dz 

Sa 

dz 

3p 


The  equation  of  the  tangent  plane  to  the  surface  is 


x-x. 


y-yi 


z-zn 


dx 

Sa 
oy 

dz 

5a 


dx 

dy 

dz 


-  0 


(163) 


or 


<"i>  $$7   *  <"i>    &     *  "-xJ  «)  "  ° 
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This  equation  permits  us  to  interpret  geometrically  the  conditions 
in  which  the  Jacobian  J«0  (that  is  the  coefficient  of  z-z).  As  J-0 
the  equation  of  the  tangent  plane  reduces  to 

[  <«i>  Sfeff    ♦  <«x>  i&   -  •  <*» 

vhich  means  no  matter  how  the  z.  changes  the  equation  of  the  tangent 
plane  remains  constant.  Therefore,  the  tangent  plane  is  vertical,  and 
also  the  projection  of  the  tangent  plane  on  the  x-y  plane  is  repre- 
sented by  ^(x,y)  «  0  which  is  a  locus  of  singular  points.    The  graphic 
interpretations  are  shown  as  follows:  (See  Figure  k^>). 
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FIG.  45  - 


GRAPHICAL        REPRESENTATION  OF 

CAUCHY       PROBLEM        WHEN       J    =     0 
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NUMERICAL  SOLUTION  AND  EXAMPLES 

The  main  purpose  of  this  section  is  the  numerical  solution 
of  the  equations  ,(71)   and  (73),  describing  the  runoff.  The  last 
paragraph  gives  the  complete  study  of  the  hyperbolic  equations  by  the 
method  of  characteristics.  Since  the  derived  overland  flow  equations 
(  71)  and  (73)   are  of  hyperbolic  type  they  can  in  principle  be 
solved  by  the  method  of  characteristics.  But  for  design  calculations 
it  is  found  there  are  some  mathematical  and  practical  difficulties 
which  make  this  method  difficult  to  use.  The  difficulties  are  as 
follows : 

1.  The  initial  condition  makes  the  Jacobian  equal  to  zero 
along  the  x-axis  in  the  x-t  plane. 

2.  The  upstream  condition  shows  that  the  Jacobian  is  equal 
to  zero  along  the  t-axis  in  the  x-t  plane. 

3.  Since  the  velocity  and  depth  of  overland  flow  is  small, 
the  Manning' 8  formula  can  not  be  used  to  evaluate  the 
friction  slope  S  .   It  is  necessary  to  evaluate  the 
friction  coefficient  in  each  computation. 

k.     The  friction  coefficient  in  the  transition  zone  is 
uncertain  and  its  evaluation  is  uncertain. 

The  first  and  second  are  the  mathematical  difficulties  which  can  be 
overcome  by  using  the  laminar  sheet  flow  approximate  solution  for 
the  early  stage  of  the  runoff  computation  and  then  use  the  method  of 
characteristics  to  solve  the  rest  of  the  problem.   But,  the  evaluation 
of  the  friction  coefficient  remains  uncertain. 
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The  approximate  solution  of  the  laminar  and  turbulent  sheet 
flow  given  in  the  previous  theoretical  analysis  seems  applicable 
for  practical  design.  The  computation  is  simple  and  the  results 
are  reasonable.  Two  complete  sets  of  design  calculations  are  given 
in  this  paragraph  and  also  a  computation  of  one  point  was  made  by 
using  the  method  of  characteristics. 

Numerical  methods 
1.  Laminar  sheet  flow  approximate  solution. 
The  solution  of  laminar  sheet  flow  is  obtained  by  using  Eq.  (6l) 

*  -  \  (1  "  15s  '> 

The  numerical  calculation  gives  the  depth  of  overland  flow 
at  the  outlet  in  any  time  interval  and  thus  the  mean  velocity  of 
flow  or  the  discharge  at  the  outlet  can  be  determined  from  Eq.(  55) 

"So  2 

u  «  y 

5v 


or 


5v 


2.  Turbulent  sheet  flow  approximate  solution 

The  solution  can  be  obtained  by  using  Eq.  (  70) 

Ay  -  1.585  [1  -  I  y  1.715]  At 


where  V7 

(gso) 

K  -   7.15   Yfi — 


v 
The  numerical  calculation  gives  the  depth  of  overland 

flow  for  the  outlet  in  any  time  interval  and  thus  tha  mean  velocity 
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and  discharge  of  the  outlet  can  be  determined  from  Eq.  (65)  and  (66) 

v    0.715 
u  -  Ky 

and 

5.  Method  of  characteristics 

As  derived  in  the  theoretical  analysis,  the  general 
hydrodynamical  equations  of  overland  flow,  Eqs.  (71)  and  (75)  are  the 
equation  of  continuity  ,  ^  +  u  g  +  y  g  -  i  and  the  equation 

of  motion,      |^  ♦  u  &      +  gP    -  g(S  -  S,)  -  ±-2   . 

•      dt      dx     ttdx    6Vo    f '    y 

Comparing  vith  the  tvo  quasi- linear  equations  (107)  and  (108) 

du  ^  a  du  .  „  by 


*$  ♦»£  -c§  ♦»£  •  * 


one  can  find  the  corresponding  value  of  the  coefficients: 

iu 
A  -  1     B»u    C  =  0   D«g    E»gSQ-gSf-  — 

Al  -  °    Bl"  y    °1  "  X   Dl  *  U  El  '  * 
The  characteristic  directions  may  be  found  by  Eq.  (122).  They  are 

£  «u+/gy  -u*  C  (166), (167) 

at     —         - 

For  y  >  0,  there  are  two  characteristic  directions  so  the  equations 
are  of  the  hyperbolic  type. 

The  two  characteristic  equations  may  be  determined  by  Eq.  (150) 
They  are  found  to  be: 
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£  (.  .  a.)  .  „.o  -  sf)  .  is.  is  {Jgj 

Using  the  finite  time  numerical  solution  used  by  Chow^    )for  un_ 

steady  open  channel  flow,  the  characteristic  equation  may  be  written  as 
(the  subscripts  u  and  d  refer  to  the  up-  and  down- stream  points  at 
time  t-dt) 
U  ♦  2C  »  Uu  ♦  2Cu  ♦  J   [g(So  -  Sf)  -  22  ♦  i£  ]dt  (170) 

U  -  2C  -  Ud  -  2Cd  ♦  /  [g(8o  -  Sf)  -  1H  .  l£  ]dt  (171) 

For  a  small  value  of  At,  Eq.  (170)  may  be  expressed  as 

U  ♦  2C  *  G  +  K,  (172) 


where 


G   »  U  +  2C  +  Y  (173) 

u      u     u    u  v  '  "" 


U      C 
Ku.i(gSo-gSfu-l7%l-»JAt  (17M 

Similarly,  Eq.  (171)  may  be  expressed  as 

u  -  2C  -  Gd  +  K2  (176) 


where 


Gd  '  Ud  -  2Cd  +  Kd  (177) 


^-^«B.-«Bfl-*  yV1  ^]Z"  U78) 

K2-|  [gSQ-  gSf  -  l|-  i  I  J  At  (179) 

Then  u  and  c  may  be  obtained  by  following: 

u  -  -|(Gu  ♦  Gd  ♦  Kx  ♦  K2)  (180) 


16  J* 

C  "  i   (Gu  "  Gd  *  Kl  "  K2>  (181) 

The  subscripts  u  and  d  refer  to  the  upstream  and  downstream  point 
passed  by  the  two  characteristic  directions. 

Numerical  examples 
1.  Laminar  sheet  flow  approximate  solution 

Assume  the  rainfall  intensity  i  «  3  inches  per  hour 
slope  of  the  smooth  plane   S  »  0.01 
length  of  the  plane        L  «  150  feet 
Find  the  outflow  at  L-25,  50,  100,  150  feet,  and  the 
overland  flow  profiles. 

The  solution  was  obtained  from  the  numerical  calculation  of 
Eq.  (61)  as  shown  in  Table  (  17  ).  The  outflow  at  L  -  25,  50,  100 
and  150  feet  was  shown  in  Fig.  (k6)   and  the  overland  flow  profiles 
was  shown  in  Fig.  (^7). 
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2.  Using  the  method  of  characteristics  to  solve  the  point 
(t-120,  x«150)  of  the  previous  example 

i.  Based  on  the  result  as  given  in  Example  1,  one  can  try 
two  characteristic  directions 

^  -  1.2,  0.23 
From  the  overland  flow  profile  given  in  Fig.  (1+7),  one  can  find  the 
u,  y  and  c  for  the  upstream  and  downstream  point 

u  -  0.510       y=  0.0072       c  «  0.U82 
u  u  u 

ud  «  0.5^        yd  «  0.00738       cd  -  0.1+87 

ii.  For  the  upstream  point,  the  Reynolds  number  may  be  cal- 
culated as  R  «  350.  Use  f  -  O.0V7  (referring  to  the  f-R  diagram 
given  by  Chow^52'  and  from  Eos.  (173)  and  (ilk),   it  is  found  that 

uu        ,  luu        .       icu 

s  f  75—  At  -  6.71     -7- At  =  0.11+7     — -  -0.139 
oyu  yu  u 

K  »  I.U7I         Gu  -  2.930 

iii.      Similarly,   for  the  downstream  point 

R  «  1+00,       Use  f  -  O.oM+5 

e     * 

From  Eqs.  (177)  and  (178),  it  is  found 

11  c 

f -A.  At  -6.57   i  — At  -0.151   i-^- At  -0.137 
8yd  yd  yd 

Kd  -  1.1+01  Gd  -  0.967 

iv.      From  Eqs.    (175), (179  ),(l8l)  one  can  calculate  c 
c  -  0.523 
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v.  From  Eq.  180,  use  trial-  and-error  method  for  Reynolds 
number  and  friction  coefficient  f ,  the  mean  velocity  may  be  deter- 
mined 

u  •  0.782 

vi.  Check  the  assumed  characteristic  directions,  it  is 
found  that 

|£  «  1.17  and  0.2, 

which  is  sufficiently  close  to  what  was  assumed.  The  result  of  this 
computation  agrees  well  with  the  result  given  by  example  1, 

3.  Turbulent  sheet  flow  approximate  solution 

Assume  the  rainfall  intensity  i  ■  3  inches  per  hour 
8 lope  of  the  smooth  plane     S  ■  0.01 
length  of  the  plana         L  ■  500  feet 
Find  the  total  overland  hydrograph  at  the  outlet,  L  «  500  feet. 

The  solution  was  obtained  by  numerical  solution  of  Eq.  (  6l) 
for  laminar  sheet  flow  and  Eq.  (  70  )  for  turbulent  sheet  flow.  After 
the  flow  reaches  the  equilibrium  state,  the  recession  curve  of  the 
overland  flow  hydrograph  is  calculated  by  setting  i»0  in  Eqs.  (  6l) 
and  (  70  ).  The  total  computation  is  given  in  Table  18  ,  and  the 
overland  hydrograph  is  shown  in  Fig.  (U8). 
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CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FUTURE  RESEARCH 

1.  The  use  of  the  method  of  characteristics  for  the  runoff  study  is 
theoretically  sound.  But,  for  practical  calculations  it  is 
necessary  to  evaluate  friction  coefficient  f .  Since  sufficient 
experimental  data  are  not  available  to  assure  a  definite  rela- 
tionship between  the  flow  and  the  friction  coefficient  f  for  the 
range  of  Reynolds  number  normally  encountered  the  computations 
are  uncertain. 

2.  The  approximate  solutions  for  laminar  and  turbulent  sheet  flows 
give  reasonable  results  for  the  overland  flow  profiles  as  well 
as  the  overland  hydrograph  at  every  point  of  the  plane.  These 
solutions  are  reliable  only  for  smooth  planes  and  may  be  used 

as  a  good  estimate  of  the  runoff  from  smooth  surface  such  as  paved 
highway,  street  surfaces,  and  airfield  runways. 

3.  The  rainfall  excess  does  not  need  to  be  constant.   Since  finite 
time  increments  are  used  in  the  numerical  solution,  different 
values  of  the  rainfall  excess  can  also  be  used  in  the  computa- 
tion by  using  average  values  for  each  time  interval.  Therefore, 
the  overland  flow  hydrograph  from  a  design  storm  with  variable 
intensity  and  infiltration  rate  may  be  determined  easily  without 
adding  much  difficulty  to  the  computation. 
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4.  The  approximate  solutions  for  the  laminar  and  turbulent  sheet 
flows  are  simple  and  good  for  practical  design.  However,  this 
solution  is  only  valid  for  smooth  surfaces.  The  same  approach 
could  be  extended  to  rough  surfaces  when  a  sound  knowledge  of 
roughness  effect  becomes  available.  Experimental  tests  and 
actual  field  investigations  are  necessary  before  the  theoretical 
calculation  could  be  made. 

5.  An  experimental  test  concerning  the  determination  of  runoff  from 
a  smooth  plane  produced  from  rainfall  is  the  only  way  to  prove 
the  accuracy  of  the  approximate  solution  of  the  laminar  and 

turbulent  sheet  flow  developed  in  this  study.  Figure  (**9) 

(35) 
shows  the  comparison  of  the  results  obtained  by  Izzard 

and  of  those  obtained  by  the  approximate  solutions. 

6.  The  computation  of  the  method  of  characteristics  will  be  simple 

and  the  results  will  be  reliable  for  high  Reynolds  numbers  for 

which  Manning's  formula  may  be  used  to  evaluate  the  friction 

(32) 
term  as  done  by  Chow  for  unsteady  flow  in  open  channel.  For 

high  Reynolds  number  the  friction  coefficient  may  be  assumed 

to  be  constant,  and  no  trial- and- error  calculations  are  needed 

to  estimate  the  friction  coefficient. 

7.  Since  the  computation  by  the  method  of  characteristics  is  compli- 
cated and  time-consuming,  the  use  of  the  digital  computer  should 

(  32) 
be  considered.  As  indicated  by  Chowv    ',  the  programming  can 

be  done  for  the  solution  of  the  flood  routing  problem  through 

open  channel  by  the  method  of  characteristics.  A  similar  program 
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may  be  developed  for  the  computation  by  the  method  of  charac- 
teristics for  this  study. 

8.  The  actual  physical  conditions  of  watersheds  are  substantially 
more  complex  than  the  ideal  plane  assumed  in  this  study  of  over- 
land flow.  The  natural  surfaces  are  very  irregular.  The  over- 
land flow  is  collected  in  small  rills  which  flow  to  larger 
streams. 

9.  As  mentioned  in  the  introduction  that  the  determination  of  the 
peak  discharge  from  a  watershed  requires  the  overland  flow  cal- 
culation for  the  small  subwatershed  and  the  flow  from  each  out- 
let must  then  be  routed  through  the  channels  to  the  outlet  of 
the  watershed. 

10.  Future  research  is  recommended  for  the  evaluation  of  the  fric- 
tion coefficient  especially  for  the  flow  with  Reynolds  number 
less  than  5,000.  This  will  require  sets  of  experimental  teats 
of  overland  flow  with  different  roughness  of  the  plane.   Infil- 
tration studies  are  necessary  for  actual  watersheds  in  order  to 
assure  a  good  accurate  estimation  of  the  effective  rainfall  from 
a  total  storm. 
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APPENDIX  A 
Results  of  frequency  analysis  of  the  32  watersheds  studied 
with  gaging  station  information 


(S-l)  Bean  BlosBom  Creek  at  Dolan,  Ind, 


(23) 


Location,—  Lat  39°14'30",  long  86°29'57",  in  SW^  sec.  2,  T.  9  N„,  R.  1  W   on 

dovmstreatn  aide  of  right  pier  of  highway  bridge  at  Dolan,  17,5  miles'upstream 
from  mouth,  v 

Drainage  area 100  sq  mi, 

Cage.-   Nonrecording  gage  Apr.  3,  1946,  to  Sept.  27,  1951;  recording  gage  thereafter 
D.itum  of  gage   is   576.41  ft  above   mean  sea  level,  unadjusted, 

age-discharge  relation.-Defined  by  current-meter  measurements.     Discharge  adjusted 
for  rate  of  change  of  stage  above   5  ft       Only  annual  maximums  adjusted  prior 
to  installation  of  recording  gage. 

Flood  stage, — 15  ft. 

Remarks.  .-Flow  regulated  since  April  1953  by  Bloomington  liesei-voir  (capacity, 

4.640,000,000  rallona)  7J  miles  upstream;  peak  discharges  probably  not  materiallv 
affected. 


Peak  Stages  and  Instants 


Annual  Pe3k   Discharge 


Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

V.'ater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1946 

Kay  16,  1946 

13,0 

1,830 

1953 

Mar    4,  1953 

LI  07 

1,320 

1947 

June  2,   1947 

17  8 

9,420 

1954 

Kay     2,  1954 

5  45 

361 

1948 

Mar. 27,   1948 

13- 5 

2,110 

1955 

Apr. 13,  1955 

LI  63 

1  390 

1949 

Jan,   5,  1949 

17,9 

9,060 

1956 

Kay  28,  1956 

L2.93 

1,740 

1950 

Jan.  4,  1950 

17.75 

8,740 

1957 

Kay  22,  1957 

15.78 

4,270 

1951 

Jan.21,  1951 

15.50 

3,700 

1958 

June  14,  1958 

3,040 

1952 

May  24,  1952 

16.12 

5,100 

1959 

Jan.  21,  1959 

5,480 

BEAN      BLOSSOM      CREEK      AT       DOLAN    ,     INDIANA 

Return  Period   (years) 
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(S-2)     CILfty  Creek  at  Hartsville,   Ind. 

Location.— Lat  39°16,25'',  long  85°42 '  10'' ,   in  NWj.  sec.  36,  T.  10  N. ,  R.   7  E. ,  at 

dovnai-ream  side  of  left  abutment  of  highway  bridge,  a  quarter  of  a  mile  north 
of  Hartsville,  and  5  miles  upstream  from  Duck  Creek. 

Drainage  area. — 88.3  sq  mi. 

Gage. — Nonrecording  gage  Feb.,   12,  1948,  to  Sept.  23,  1952;  recording  gape  thereafter. 
Datum  of  gage  is  677.34  ft  above  mean  sea  level,  datura  of  1929. 

Stage-discharge  relation. — Defined  by  current-meter  measurements  below  6, COO  cfs. 

Historical  data. -Flood  of  1913  on  Clifty  Creek  reached  a  stage  of  about  3  ft  higher 
than  the  KcKinley  (1897)  flood  according  to  a  report  in  the  Evening  Republican 
of  Columbus,  Ind.  dated  Bar.  25,  1913.  (The  preceding  statement  was  apparently 
for  the  Petersville  area,  about  6  miles  downstream  from  Hartsville). 


peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

V/ater 
Year 

Date 

Gage 
Height 

Discharge 

cfs 

1913 

Mar,  25,  1913 

25.1 

1954 

toy  27,  1954 

4.17 

635 

1948 

Mar,  27.  1948 

8-48 

3,710 

1955 

July  8,  1955 

6.24 

1,760 

1949 

Jan.  5,  1949 

13-4 

8,100 

1956 

June  22,  1956 

11,10 

5,890 

1950 

Jar..  4,  1950 

11 .8 

6,520 

1957 

July  4,  1957 

9,28 

4,270 

1951 

Nov.  20,  1950 

8  9 

3,910 

1958 

Kay  6,  1958 

2,700 

1952 

Jan  26,  1952 

11. 5 

6,250 

1959 

Jan. 21,  1959 

11,300 

1953 

Mar.  4,  1953 

5.57 

1,370 

CLIFTY       CREEK        AT        HARTSVILLE     ,    INDIANA 

p£turn  Period   (years) 
2  5  10  2S  50  I 


<  400C- 


■ft  ■  ■   ■  '  '■* 


Probability  -  percent 
Reduced  Variate       (Y) 


99"76         9979 


^ 


-i 


(S-3)     North  Fork  of  Vernon  Fork  near  Butlerville,  Ind. 

Location.—  Lat  39°02'55",  long  85°32<40",  In  SE£  sec.  17,  T.  7  H„,  H.  9  B.,  on 
left,  bank,  0  3  mile  downstream  from  Muacatatuck  State  School  dam,  lj  mlleo 
downstream  from  Brush  Creek,  and  2  miles  northwest  of  Butlerville. 

Drainage  area. — 37*3  sq  mi- 
Gage.— Nonre  cording  gage  Feb.    16,  1942,  to  Aug.  18,  1942;  recording  gare  thereafter. 
Datum  of  gage  is  669.40  ft  above  mean  sea  level,  datum  of  1929. 

Stage-discharge  relation. — Defined  by  current-meter  measurements. 
Flood  stags. --11  ft. 


Peak  Stages  and  Instaneous  Annual  Peak  Discharge 


Water 
Tear 

Date 

Gage 
Height 

Discharge 

Cf3 

Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

1942 

Apr     9. 

1942 

8  94 

2,560 

1951 

Nov. 

20,  1950 

15.98 

8,030 

1943 

Mar  16, 

1943 

17  79 

9,910 

1952 

Jan 

26,  1952 

13,18 

5,300 

1944 

Apr  11. 

1944 

12.63 

4,780 

1953 

Mar. 

4,  1953 

10,34 

3,260 

1945 

Mar.   6, 

1945 

18.72 

10,900 

1954 

Jan. 

1,  1954 

5  58 

840 

1946 

Feb  13, 

1946 

15-95 

8,030 

1955 

Feb. 

27,  1955 

12.05 

4.300 

1947 

June  2, 

1947 

14.30 

6,330 

1956 

Kay 

28,   1956 

16.23 

8,330 

1948 

Mar  27, 

1948 

15.12 

•7,130 

1957 

May 

22,  1957 

17  04 

9,080 

1949 

Jan  24, 

1949 

18,73 

10,900 

1958 

July 

22,  1958 

7,730 

1950 

Jan     4, 

195C 

17  90 

10,000 

195? 

Jan 

21,  1959 

26,200 

NORTH       FORK      VERNON        FORK       NEAR       BUTTERVILLE  ,     INDtANA 

Return  Period  (ytars) 
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Reduced  Variate 


(S-t)     Hart  ditch  at  Minister,  Ind. 

Location.—  Lat  41°33'40n,  long  87°28'50",   in  N  l/2  sec,  20,  T.  36  N,,  R.  9  W.,  on 
left  bank  at  city  limits  of  Munster,  a  quarter  of  a  mile  downstream  from  U,  3. 
Highway  41,   and  0.4  mile  upstream  from  mouth. 

Drainage  area  —69.2  eq  mi. 

Gage .  —Recording.     Datum  of  gage  is  591.21  ft  above  mean  eea  level,  datum  of  1929. 

Stage  discharge  relation. — Defined  by  current-meter  measurements.     Dredging  operations 
assumed  to  have  occurred  between  April  1944  and  April  1945,  and  subsequent  filling 
have  effected  high-water  rating.  Backwater  from  Little  Caiumat  River  and  possi- 
bly from  overbank  return  affects  stage  at  gage  at  times  during  periods  of  ex- 
tremely high  flow 

Flccd  Stage.-  -7  ft 

Remarks .  -  Hart  ditch  is  tributary  to  Little  Calumet  River   At  this  point  low  flow 

of  Little  Calumet  River  runs  west  into  Calumet  Sag  Channel  or  into  Lake  Michigan 
through  Grand  Calumet  River;  floodflow  at  times  runs  east  into  channel  storage 
or  through  Burns  ditch  to  Lake  Michigan, 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1943 

Mar.  16, 

1943 

6-95 

2,280 

1952 

June  14, 

1952 

4,39 

1,190 

1944 

Mar,  15, 

1944 

7  23 

2,420 

1953 

Mar.  15, 

1953 

3  84 

960 

1945 

May   8, 

1945 

3  73 

x,270 

1954 

Mar.  25, 

1954 

4  25 

1,110 

1946 

Jan,  6, 

1946 

2.38 

780 

1955 

Oct,  11, 

1954 

7,83 

2,600 

1947 

Apr   6, 

1947 

6,17 

2.490 

1956 

toy  11, 

1956 

5  27 

1,550 

1948 

May  11, 

1948 

5  60 

1,950 

1957 

July  14, 

1957 

7  60 

2.06C 

1949 

Feb.  13, 

1949 

3,00 

850 

1953 

June  10, 

1958 

960 

1950 

Dec  22, 

1949 

4.83 

1,570 

1959 

Apr-  28, 

1959 

2,670 

1951 

May  11, 

1951 

5  01 

1,430 

*3000 


HART        DITCH        AT        MUNSTER      ,  INDIANA 

Return  Feriod  (yvirs) 
2  ?  Ifl  25  5C 


500  ieuu 


5    16     26    30  tI»0  56   60     70       ETC 

51  }  Probability  -  percent 

■  Jm.,!,,,,^.,'?,..!....! — t 


9775         9979 


^- 


Reduced  Varlate        lYJ 


i""1""!, 7 


(S-5)  Salt  Creek  near  KcCool,  Ind. 


Location,  Lat  41J35-4S»;  long  87°03'40«,  in  SEi  sec  6t  T  36  N  ,  R  6  W   on  left 

bank  or.  downstream  side  of  highway  bridge,  50  ft  downstream  irom  New  York  Central 
Hailroad  bridge,  U  miles  north  of  IfcCool,  and  1  5  miles  upstream  from  Little 
Calumet  River 

Drainage  area  -  -78  7  sq  mi. 

Gage  -Nonrecordit*  gage  May  5,  1945,  to  July  24,  1955;  recording  gage  thereafter 
Da»um  of  gage  is  594-10  ft  above  mean  sea  level,  datum  of  1929  (levels  by 
Indiana  Flood  Control  and  Viator  rtesources  Commission)., 

Stage  discharge  relation. —Defined  by  current-meter  measurements  below  2,300  cfo 
Flood  Btage  -  10  ft 


Feak  Stages  and  instantaneous  A.-.nual  Peak  Discharge 


Water 
Year 


Date 


Gage   Discharge 
Height     efs 


V'ater 
Year 


1945 
1946 
1947 
1948 
1949 


June  29,  1945 
June  13:  1946 
Apr  5,  1947 
Kay  11,  1948 
Feb,  14,  1949 


10  SO 

11  ?.7 
11-83 

12  3 
9.28 


1950 

Deo  22,  1949 

12  02    1,' 

1951 

Kay  11.  1951 

10,78     ' 

1952 

Kov.  14.  1951 

10  63 

990 

1,280 

1,580 

1,910 

525 

,700 

970 

912 


Date 


Gage 
Height 


1953 
1954 
1955 
1956 
1957 
1958 
1959 


Kar  16,  1953 
Kar  26,  1954 
net.  11,  1954 
Apr,  29,  1956 
Apr,  27,  1957 
Nov.  15,  1957 
Apr  28,  1959 


8  16 
10  48 
14  12 
11,26 

9  81 


Discharge 
cfs 

454 

910 
3.180 
1,280 

725 

456 
1,200 


SALT      CREEK      NEAR      McCOOL        ,      INDIANA 

Return  Feriod  (years) 
2  5  10  25  50  100 


500  loop 


Reduced  Variate 


(S-6)     Little  Calumet  River  at  l-orter,   Ind. 

Location— Lat  41°37v18",   Ions  87°05'13",  in  XB  1/4  sec.  34,  I.  37  N,,  K    6  W\,  near 
center  of  span  of  downstream  aide  of  highway  bridge,  three-quarters  of  a  mile 
northwest  of  Porter,  and  4,5  miles  upstream  from  Salt  Greek, 

Drainage  area. — 62.9  sq  mi. 

Gage  — Nonrecording  gage  May  5,  1945,   to  June  25,  1952;  recording  gape  thereafter. 
Datum  of  gage  is  603.43  ft  above  mean  sea  level,  datum  of  1929* 

Stage-discharge  relation. — Defined  by  current-meter  measurements  below  2,500  cf s , 
Hating  subject  to  changes  throughout  ranpe  of  stage. 

Flood  etage,— 7  ft. 


feak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

V.'ater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1945 

June  2&t 

1945 

9.88 

2,440 

1953 

May     23, 

1953 

6. 64 

5a 

1946 

June  13, 

1946 

6  99 

715 

1954 

Apr     26, 

1954 

8.32 

1,170 

1947 

Apr       5, 

1947 

9. 42 

2,140 

1955 

Oct.  10, 

1954 

11  66 

3,110 

1948 

May     11, 

1948 

9. 10 

1,960 

1956 

Apr.  29, 

1956 

8  67 

1,370 

1949 

Kay     20, 

1949 

6,88 

690 

1957 

Ap:%   27, 

1957 

7.65 

848 

1950 

Dec     22, 

1949 

8.72 

1,720 

19  56 

Feb.   28, 

1958 

490 

1951 

May    11, 

1951 

8.11 

1,360 

1959 

Apr.   28, 

1959 

1,420 

1952 

Nov,  14, 

1951 

7.92 

1,060 

o 


LITTLE        CALUMET       RIVER      AT       PORTER    ,    INDIANA 

ii--rr.  Ftr- od  (;.-j=rsJ 
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Reduced  Varlate  IYJ 


(5-7)     Ceaar  Creek  at  Auburn,    Ind 

Location  --lat  41°2i',  lor.g  S5°03',   lr  SK  1/4  sec.   29,  T     34  K    .    A     13   S   .   ifar  center 
of  span  on  upstreara  sidt  of  Ninth  street  Bridge  in  Auburn  and  2  miles  upstream 
from  Peokbart  ditch 

Drainage  area  -~93  si  r.d. ,  appro::irnately, 

Cage       Nonrecordir.g  ga»e  July  30,   19i3,   to    ;ept,   3C,  1953;  recording  gage  there- 
after       Dfture  of  gage  iE   847.14  ft  above  aaan  sea  level  (city  of  Auburn  bench 
mark). 

Sta^e-discharge  relation.  — Defined  by  current-TT.etcr  measurements. 

?lood  stace„--6  ft.      *^v' 


Peak  Stages  and  Instantareous   Annual  Peak  Discharge 


Water 
Ye.ir 

Date 

Gage 
Height 

Discharge 
cfa 

'..rat9r 

Year 

Date 

Rage 
Height 

Discharge 
cfs 

1943 

May 

1943 

9  8 

1,470 

1952 

ICar 

11, 

1952 

8  45 

900 

1944 

Apr.  12 s 

1944 

9-3 

1,230 

1953 

Mar 

4, 

1953 

5  80 

471 

19s  5 

Kay  IP . 

3945 

9  13 

1,150 

1954 

tar 

25, 

1954 

7  57 

707 

1946 

June  13. 

1946 

8  53 

916 

1955 

Jan 

6.. 

1955 

7.61 

707 

1947 

Apr.  21 , 

1947 

9  02 

983 

1956 

A  or 

30, 

1956 

8,85 

1.050 

1948 

Feb,  28, 

^948 

8  53 

9C5 

1957 

ipr 

6, 

1957 

6  89 

651 

19*9 

Feb  16, 

1949 

9.21 

995 

1955 

Dec 

20, 

1957 

540 

1960 

Apr   5 , 

1950 

9  90 

1,5?0 

1959 

~eb 

14, 

1959 

390 

1951 

^eb  ??. 

1951 

s  -: 

9/0 

CEDAR         CREEK       AT       AUBURN       ,      INDIANA 

Return  Period  (years) 
2  5  10  Z'j  50  100  200 


500     1000 


Reduced  Varlate    IY) 


(S  3)     West    Cree* 


icrneider,    Inci 


Legation,-  -Lat  U°12*52;',  long  87029f36«.    in  NK  1/4  f;E  lA  sec.    L9,  T     32  11   ,   R-    9  W  , 
en  loft  bank  at  downstream  side   of  county  highway  bridge,   1.2  miles  upstream 
I'rcm  Singleton  ditch  and  2  3/U  rail3s  northwest  of  Schneider 

Drain' tit.   aroa,--5J*   5  sq  mi. 

Gage      -r  onrecordir.r  ga»e  July  29,   19i-»i    to  Dec     31,   1951,   and  Jan     1,   1954*  to  June  10 
1956;   recording  f,a,*e  since  Junr  11,   1956-      Datum  of  gage   is  627  36  ft  above 
rean   sea  level,   datum  of  1929   (levels  by  Soil  Conservation  Service). 

Sta»e  discharge  relation , --Defir.ee   by  current -meter  measurements 

Flo  >d  stare       7  ft. 


Peak  Stages  and   Instantaneous    \nnual  heak  Discharge 


Water 

Year 

Date 

Gape 
Feirht 

Discharge 
cfc 

'.'ater 
Year 

Date 

lage 
Height 

Discharge 
efa 

19c9 

Feb 

13,  1949 

4  53 

504 

1956 

Feb  25,  1956 

5  42 

710 

19S0 

jac 

22,  1949 

£  56 

1,050 

1957 

July  13,  1957 

7  02 

1,250 

1951 

Feb 

19,  lc51 

5  52 

733 

1953 

June  9,  1958 

794 

1954 

Mar 

25.  1954 

6.10 

i.COO 

1959 

Apr  2fl,  1959 

1,200 

1955 

Oct 

10,  1954 

e.06 

1  840 

CREEK       NEAR         SCHNEIDER 

Return  Period  (years) 
5  10  25 


"2 T 

Reduced  Variate         (Y) 


(S-9)     Iroquoi3   'liver  at  Rosebud,   Ind 

Location    ~Lat  41°02- ,  long  87°11- ,   in  S'..'  1/4  sec,   24,  T     30  N.,   R     7  W,,   100  ft 
dovnst.rea.-3  froa  bridge  on  county  road,    half  a  aile  north  of  Rosebud,  half  a 
mile  downstream  from  confluence  of  3wair.  and  Dexter  ditches.    1.5  riles  upstream 
from  Davidson  ditch,   and  2  miles  east  of  ;  arr 

Drainage  area     -30  3  sq  mi 

Oage,--Nonrecording  gap:  July  12,   194f,  to  Sept     )'.■,   1953;   recording  gage  thereafter 
Datum  of  gage   is  661  47  ft  above  mean  sea  level,   datum  of  1929, 

Sta^e-discharge  relation  —  Defined  by  current  zeter-  measurements  below  330  cfs 

Flood  stage  —10  ft 


Peak  Stages  ar.d   Instantaneous   Annual   Peak  Discharge 


V.'ater 
Year 

Date 

Rage 
Hel.Tht 

_ ischarge 
cfs 

Uater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1949 

Feb  If,, 

1949 

6  15 

254 

1955 

Jan   6,  1955 

4  84 

126 

1950 

Apr   ^, 

195C 

p  3 

422 

1956 

Apr  29,  1956 

6  65 

225 

1951 

July  9, 

1951 

7  2 

235 

1957 

Apr  28,  1957 

7  90 

29C 

1952 

Apr  23, 

1952 

7  3 

263 

1953 

June  1C,  1958 

308 

1953 

Mar  15, 

1953 

5.75 

:=>; 

1^59 

Feb,  10,  1959 

343 

1954 

Mar  25, 

1954 

4  59 

13C 

IROQUOIS       RIVER      AT       ROSEBUD    ,       INDIANA 

Return  Period  (years) 
2  5  10  25       50  130 


?0    66     76  '    'Bo 90 

}  Probability  -  percent 


Reduced  Varlate        (Y) 


.-•   10)     Bice  ditch  near  South  :'.ariun,   Ind 

Location-  Ut  40°52-,   long  S7°06',   on  line  between  sees     15  and  22 .  T     28  N 

K     OH  ,    on  left  bank  at  upstream  side  of  bridge  on  State  Highway  16     ;•  miles 
upstream  :rom  Big  Slough  Creek,    3  niles  southeast  of  South  ,',arion,   anc   5  miles 
southeast  of  Rensselaer 

Drainage  area. --22  6  sq  ml 

Sage  - -Nonrecording  gage  Dec     31,  1948,  to  Aug     4,   1955;   recording  gage  thereafter 
Datum  of  ^age  is  o>3-30  ft  above  moan  sea  level,  datum  of  1929 

Sta^e  discharge  relation. --Defined  by  current-meter  measurements 


peak  Stares  and    Instantaneous  Annual   Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

','ater 
Year 

Date 

Gage 

Height 

Discharge 
cfs 

1949 

Feb     15, 

1949 

9  09 

410 

1955 

June  11,  1955 

10  16 

353 

1950 

July  19, 

1950 

10,06 

490 

1956 

Apr     29,   1956 

10  75 

504 

1951 

July     9, 

1951 

11  43 

610 

1957 

July  13,  1957 

10  36 

458 

1952 

June  14 ? 

1952 

10  80 

556 

1958 

June  13,  1958 

780 

1953 

July     5, 

1953 

R  65 

'   374 

1959 

Feb.   10,  1959 

480 

1954 

Jane  22, 

1954 

8  12 

339 

DITCH         NEAR        SOUTH        MARION    ,   INDIANA 

Return  Feriod   (years) 
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Reduced  Varlate  (Y.) 


(S  11)     3ig  Slough  Creak  near  Collegevilla,   Ind. 

ville,   *  miles  upstream  from  ££  ^  3T X^o^  r^B^cT 
Drainage  area    --8/,  1  sq  ni 

Sap..  --Nonrecording  gaBe  July  28,   19tE„  to  Dec     31     1951tandOct     1     19-;?     t„  »„„     , 
1955;   recording  ga^e  since  Au,.    5,  1955       tit™  of  gag"  Is  637  75  ft  above  ^ 
mean  sea  level,  datum  of  1929,  aoove 

Stage-discharge  relation.    -Defined  by  current -mete.,  measurements 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Wa*ier 
Year 

Date 

Gage 
height 

Discharge 
cfs 

Vater 
Ysar 

Date 

Gage 
Height 

Discharge 
cfa 

1913 

March 

1913 

13   7 

195/* 

J-ine  22, 

195/. 

8  84 

390 

19-27 

1927 

12  5 

1955 

June  11, 

1955 

12  U 

1,100 

191.9 

Fer     15s 

1949 

11  26 

83  J 

1956 

anr     29, 

1956 

13  0 

1,1.70 

1950 

Apr       d 

lr50 

12.3 

1,180 

1957 

J  1-    It, 

1957 

12  96 

1,170 

1931 

July  10 

1951 

13   22 

l.i.50 

1953 

June  13, 

1958 

2,030 

1953 

\pr       1, 

1953 

10.50 

655 

1959 

A  r     23; 

1959 

1,030 

BIG        SLOUGH       CREEK       NEAR      COLLEGEVILLE    ,     INDIANA 

Return  Feriod  (years) 
2         5  10  Z'S       50  130  200 


'oji — i   '  J  ib  '  2b  ' 3'o " iL'o '  4o' £6"it '  'Bo""     90         95  9^  9 

,*  <  Probability  -  percent 


— r"-     '  l.1         '     't 


4 i 


Reduced  Variate  (Y) 


(3  12)     Carpenter    ;.-eek  at  Zcypi.,   Ind . 

Location  •  -Ut  40°52' ,   long  S7°12  ,   en  lir.e  between  S'^  sec.   15  and  JlWj  sec,   22. 
T.  28  l:   ,   h,   71   ,   on  left  bank   on  downstream  side  of  bridre  on  State  Highway 
16,    2  3/1   miles  upstream  from  ncuth,  and  4  miles  southwest  of  Collegevllle 

Drainage  Area    —48  1  sq  mi 

Gage      -Nonre cording  gape  July  26,  1943,  to  Dec.   31,  1951,  and  Oct.   1,  1952,  to  Sept     5, 
1955;   recording  gage  since  Sept.  6,  1955       Datum  01  fare  is  640-37  aoove  mean 
sea  level,   datura  of  1929. 

Stage  discharge  relation. —Defined  by  current -meter  measurements 


Peak  Stages  and   Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cf  s 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfe 

1949 

Feb  15,  1949 

10  14 

1  160 

1955 

June  8,  1955 

9  80 

984 

1950 

Apr   4,  1950 

10  3 

1,300 

1956 

Apr  29,  1956 

9  93 

1,040 

1951 

July  9.  x951 

10  92 

1,790 

1957- 

July  13,  1957 

942 

810 

1953 

July  6,  195.3 

9  21 

735 

1958 

June  10,  1958 

3,720 

1954 

June  22,  1954 

8,95 

655 

1959 

eb  10,  1959 

2,690 

CARPENTER       CREEK        AT         EGYPT     ,  INDIANA 

Return  Period  (years) 
5  10  25'  50 


130  200 sop  1000 


L     °- 


Reduced  Variat*        (Y) 


rpesanoe  P.ivr-r  nt  Oswego,   Ind 

Location       Lat  41°19'H",  long  85°W?J."1    in  KEr   KB?    sec     14,  "     33  11,  -  K  i  E       on 
left,  oank  1C  ft  downstream  .row  dan  at  Tippecanoe  Lake  outlet  in  6sw*ro,  3 
nr_lea  eas _  of  Le esburg. 

Draincpe  ar_a.--115  sq  mi 

Gage  - -Non-recording  Bare  Oct.    1,  1949,  to  Aug.   11,   1953;   recording  gage  thereafter, 
Datum  01   gape  ia  3,0  00  ft  a!:ove  scan  sea  level,  datum  of  1929 

Stage-discharge  relation    -nsrinsd  by  currsnt-™ter  measurements  below  680  cfs  and 
extended  to  1,050  cfs  b./  logarithmic  plotting 

Remarks. --Peak  discharges  affected  by  natural  storage  in  numerous  lakes  upstream 


peak  Stages  and   Instantaneous  Annual  .  eak  Discharge 


Water 
Tear 

';a^e 

Gage 
Height 

Dischsrpe 
cfs 

Water 

Year 

Date 

Gage 
-eight 

Discharge 
cfs 

1943 

Kay 

21, 

1943 

9  4 

1,050 

i  1955 

•ct       17,   1954 

9  65 

7~0 

1950 

Apr 

3-  10, 

195C 

3  62 

6  JO 

!  1956 

Hay  5     6,  1956 

3  08 

4  50 

1951 

Feb 

27  23, 

1951 

430 

!  1957 

pr  17  22,1957 

7  59 

315 

1952 

Jan 

30, 

195^ 

370 

|  1953 

Sept. 19-  13,1957 

383 

1953 

**ar 

22,23. 

1953 

6  ?_ 

179 

1  1959 

?B-      13,   1959 

548 

1954 

Apr 

29,30, 

1953 

7  60 

3=2 

1 

TIPPECANOE      RIVER       AT       OSWEGO     . 

Return  Period   (years) 
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Probability  -  percent 


Reduced  Variate         CY) 


J t,,,,|,..,t j 


(2)  Kississinewa 


'ear  Rldperille,  Ind 


Location  — Lat  ..0°17-,  long  8j°0C',   in  SH  1//.   »er     a     t     10  r       D     ,,   „ 

bank  10  ft  downstream  from  hirtway  bridre     0  «  in.  i      '  i  ,  E  '  on  rlght 

and  2  mile3  east  of  Ridg.nillf  '  ^  dovm3trea™  fron  Kud  Creek, 

Drainage  3rea    ■  130  sq  ni 

Gage    ■ -Konrecording  gage  Aug,   30.  1946.   to  Oct     3     TKn-   ,-„,. i< 

Datum  of  ea?e  js  96'  ??  -t   Ly^ll  ™  ,    '  ,        '   recorrting  gage  thereafter. 

"   *af!e   30   >°'  ZJ   ■»  above  mean  sen  level,  datum  of  1929. 

StaRe-diseharge  relation.- -Defined  by  current-meter  measurements  below  3,400  cfs- 
Flood   stare 10  ft 

Fi"°S  ?f  f;"1-0"1  r""ents  stated  that  »e  1913  flood  was  secondary  to  a 

flood   m  the  earljr  l930's  when  the   river  reached  an  estimated   stage  V°5  0  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

'..ater 
Year 

Date 

Gage 
height 

Discharge 
cfs 

1947 

Jan 

30, 

1947 

11  16 

2,490 

1954 

Mar    30,  1954 

7  80 

1,020 

1943 

Jan 

1, 

1948 

12  2 

3,480 

1955 

Jan.     6,  1955 

11  16 

2.490 

1949 

Jan 

5, 

1949 

131 

4,560 

1956 

Nov    16,  1955 

12  79 

4,200 

1950 

Feb 

14, 

2950 

13  4 

4,920 

1957 

June  28,  1957 

14,57 

8,830 

1951 

Feb. 

21, 

1951 

12  75 

4,200 

1958 

June  10,  1958 

13,900 

1952 

Jan. 

26, 

1952 

11.99 

3,250 

1959 

Jan.   21,  1959 

)  j.20 

1953 

Mar 

4, 

1953 

12  00 

3,250 

MISSISSINEWA      RIVER       NEAR        RIOGEVILLE      ,    INDIANA 
Return  Period   (years) 
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Reduced  Variate         (Y) 


>  .-.    car.    .13'-,    -  .    i  'eentc  rn,    -j.c 

Locatlor  -  -Lat  40°27  ,  long  35057',  en  lli;e  betwe.n  sacs  9  and  10,  T.  23  \'  .  K  5  E. 
on  left  bank  at  dowr.3tream  aid.;  o!  oridr.e  on  5ta„e  Highway  213,  1,5  nileo  south 
of  Oreent.own 

D^.-ilnage  area'- -•■  162  sq  ni;  172  sq  mi  prior  to  Junn   5,  195/. 

Gagea    -,','om-ecording  gsre  Feb.    20,  1945,   to  June  4,  1954;   recording  gage  thereafter 
Prior  to  Juno  5,  1954,  at  site  2.0  miles  downstream  at  datur  5.34  ft  lower  than 
present  datuw.      natuir.  of  present   gape  is  SC9.33   ft  above  mean  sea  level,  datum 
of  1929- 

Sta°e-discharge  relation. —Defined  by  current -reetar  measurements* 

Flood  stage.      11  it  at  both  sites. 

Historical  data.     The  following  statements  appear  in  old  newspapers  for  Kokoro, 
about,  9  miles  downstream.     September  1?66:      "Wildcat  Creek  raging;  no  train 
for  3  days."     1904:      "Greatest  flood  in  Kokono  history  " 

Flood  of  August  1371  reached  a  stage  3  inches  below  that  of  the  1913 
flood  at  a  bridge  X}_  miles  downstream  from  present  sif,e  according 'to  informa- 
tion by  local  resident  on  the  basis  of  remembered  high  water  marks  made  on  the 
same  tree 

Peak  Stages  and    Instantaneous  annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

'.Cater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1943 

May, 

1943 

15  C 

5,960 

1952 

Kar. 

11, 

1952 

11  52 

2,580 

1945 

Apr- 

1,  1945 

10-04 

1,680 

1953 

Mar 

U, 

1953 

10  34 

1  810 

1946 

Oct 

2.  1945 

9-94 

1,640 

1954 

Apr 

12, 

1954 

5  63 

450 

1947 

Apr 

30,  1947 

10  94 

2,140 

1955 

Jan 

7. 

1955 

10  00 

1  650 

1948 

Mar 

22,  1948 

11,63 

2,670 

1956 

Kay 

28, 

1956 

9-97 

1,650 

1949 

Jan 

19  1949 

13  19 

4,110 

1957 

June 

29, 

1957 

12  37 

2,260 

1950 

Jan 

4  1950 

15  3 

6,320 

1958 

June 

10, 

1958 

4,900 

1951 

Feb 

21,  1951 

10  7 

2,020 

1959 

Feb 

10, 

1959 

5,390 

WILDCAT      RIVER        NEAR         GREENTOWN    ,     INDIANA 

Return  Period  (years) 
2  5  10  25  50  100  200 


5  7C0C- 


500         inoo 


Reduced  Varlate         (Y) 


(/*;  Cicero  Creel'  near  Arcadia 

Location,- -La t  W>°11' ,  long  86°CC ,  on  line  between  sees.  IB  and  19,  T,  20  N 

K  5E-,  on  left  bank  on  downstrtam  side  of  county  bridge,  li  miles  east  of 
Arcadia,  Hamilton  County,  and  5  miles  upstream  from  little  Cicero  Creek 

Drainage  area.- -131  sq.  mi 

Gage  ---  Water  stage  recorder.     Datum  of  gage  is  815.12  ft  above  mean  sea  level 

Datui".  of  1929.      Prior  to   Dec.   ?,   1955,  wire  weight  gage  at  same  site  and  datum. 


Peak  Gtapes  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1955 

July  16,  1955 

1,280 

1958 

June  15,  1958 

2,740 

1956 

July  21,  1956 

1,540 

1959 

Feb,  11,  1959 

2,170 

1957 

June  29,  1957 

6,720 

CICERO       CREEK       NEAR       ARCADIA      ,        INDIANA 

Return  Period  (years) 

25      50 


500     1000 


Reduced  Varlate         (Y) 


(5)      Pall  Crte.     ie.u    .-'nrtville,  -  Ind.. 


Losation       La'.  39°.7  15\-   long  S5C5.= 


',   in  o;o 


17  .:       1    6  E  ,  on  right 


tank     t  downstream  side  of  brid?^  nr  3tate  Highway  233,.  1  mile  downstream  from 
lick   -Toe;  and  2  m-Ies  noiffchvest  of  "ortville 

Drainage  area       172  sh  mi 

Cage. --Nonreco. -ding  gage  Jul,"  ],   1941     to  June  16,   1942;    recording  page  thereafter 
Datum  of  jape  is  787  43  ft  above  scan  se?  level,  datum  of  1929  (levels  by 
Indianapolis  '  >ter  Co   ) 

Stage -diuchargG  relation  —  -Dei*:  r.ec  by  cuiTent  deter  measurements 

Flood  Dt^p,e        o  ft 

Historical  data, --Flood  of  1913  reached  a  stape  o-"  about  12  feet  according  to  informa- 
tion from  local  resident. 


Peak  Stares  and   Instantaneous   Annual  Peak  Discharge 


Water 
Year 

Date 

Cage 
Height 

Discharge 
cfs 

'.Vater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1942 

Mar 

17, 

1942 

7  39 

2,460 

1951 

Feb 

22, 

1951 

3.36 

4,250 

1943 

May 

18, 

1943 

9.77 

8,240 

1952 

Jan. 

27, 

1952 

6,85 

2.000 

1944 

Apr 

12; 

1944 

S  79 

5,000 

1953 

July 

6, 

1953 

S  14 

3,850 

1945 

June 

17, 

1945 

6  36 

1,840 

1954 

"far 

30, 

1954 

5  50 

1,140 

1946 

Oct 

2P 

1945 

6.87 

1,840 

1955 

Jan 

°, 

1955 

5  76 

1,260 

1947 

July 

15, 

1947 

7.25 

2,250 

1956 

Feb 

28, 

1956 

7  93 

3,130 

1948 

Mar 

24, 

1948 

7,97 

3,500 

1957 

June 

29, 

1957 

8. 12 

3.430 

1949 

Jan 

19j 

1949 

7-19 

2,350 

1958 

June 

14. 

1953 

5,040 

1950 

Jan 

4, 

1950 

• 

5,iro 

1959 

:Vo 

10 

1959 

3.010 

FALL        CREEK         NEAR        FORTVILLE 

P.eturn  Period  (years) 
2  5  10  25 


;oo 


5    lb  '  20    30  WWZ 


«  '* 


Probability  -  percent 

-J — t1-' V111 


99T9 


Reduced  Varlate 


Location   --Ut  39°46  40",   long  ii6°15  C2 ",    in  riW*   Bee     6,   T.   15  II       R     3  E.  ,   on  right 
bank  ao  downstream  side   of  bridge  on  Lynrhurst  Drive     3  0  miles  upstream  from 
Little  Eagle  Creek,    5  0  miles  west  o.'  Ilenunent  Circle   in  Lndianaoolis     and 
6  7  miles  upstream  from  mouth 

Drainage  are*--179  so.  mi. 

Gage   --Nonrecordijig  gape  Nov     18,   1938,    to   I'ar     19-   1939;   recording  r,age  thereafter 
Datum  of  gage   is  706  21    ft  above  m  :o  i  sea  lov«i ,   datum  of  1929 

Stage  discharge  relation  .—Defined  by  cur-eat  meter  measurements  below  9.000  cfs  and 
extended  above  on  basis  of  a  coctbi.net]   current Hnetor  measurement  ;md  slop-area 
measurement       Kigh-water  relation    iao   shown  a  tendency   to  shift        Discharge 
shown  for  the  1913   flood   is  an  app-v.  finite  value  based  on  logarithmic  extension 
of  an  early  rating  curve  above   9,000  cfs 

Historical  data     -The   following  in formation  was  obtained  from  a  report  on  'Tagle 

Creek  at  Indianapolis,  Channel  Improvements  for  Flood  ''ontrol"  by  Indiana  Flood 
Control  and  Water  Resources  Commission,  iated  February  1956:  "Investigations 
on  past  flooding  by  searching  newspaper  files,  interrogation  of  ?  ->cal  oersons, 
examination  of  old  surveys  and  reports,  indicate  that  flooding  occurred  along 
Eagle  Creek  in  1875,  1904,  1913,  1915,  1926.  1933,  1943,  and  1950  Mev.spaper 
accounts  of  flooding  on  other  strains  in  the  Indianapolis  area  Indicate  that 
flooding  probably  al30  occurred   in  l!  47,   1858,   1866,   1882  and  1883" 

"It  is  probable  that  the  floor's  of  1875  and  1904  were  among  the  greater 
floods  on  the  stream       f.'ewspaper  acccuntj  and  weather  records   indicate  that  the 
flood  of  July  191 5  was  nearly  aa   great  as  that  of  March  1913." 

Peak  Stages  and   Instantaneous   Annual   Peak  Discharge 


Water 
Year 

Date 

Cage 

Height 

Discharge 
cfs 

1,.  » 
i  V.ater 

1  Year 

Date 

Gage 
Height 

Discharge 
cfs 

1913 

March 

1913 

16  0 

19,000 

i  1947 

Jan. 

30, 

1947 

8  47 

3,370 

1938 

April 

1938 

14  5 

- 

|1948 

\pr 

6» 

1948 

12  26 

9,550 

1939 

Mar, 

12 

1939 

10  6 

6,610 

1949 

Jan. 

19, 

1949 

11  86 

7.250 

1940 

Mar, 

3, 

1940 

6  30 

1,860 

'1950 

Jan 

4, 

1950 

13  03 

8,670 

19a 

June 

12. 

1941 

5-77 

1.470 

|l951 

Feb 

21, 

1951 

8  47 

?  ,950 

1942 

Feb 

7, 

1942 

8.66 

4,120 

J1952 

Jf.n 

27- 

1952 

9  88 

5,520 

1943 

May 

11, 

1943 

12  17 

9,660 

!l953 

Mar 

5, 

lt53 

9  34 

4.920 

1944 

Apr 

11. 

1944 

10  61 

6,610 

|1954 

,\pr 

6, 

1954 

6  41 

2,250 

1945 

Mar 

31, 

1945 

S  80 

4,230 

1"55 

July 

16. 

1955 

6  91 

2,650 

1946 

May 

17, 

1946 

8  99 

3,860 

1956 
<1957 
;i958 
>59 

June 
Aug 
J  fin 

28 1 

28. 
8, 

22, 

1956 
1957 
1958 
1959 

13  62 

16  38 

9,920 
28.800 
8,560 
6,290 

CREEK       AT       INDIANAPOLIS    ,     INDIANA 
Return  Period   (years) 
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Reduced  Variate 


(7)     Youngs    Sreek  near  _:irbur£,    Ind 

Location. --Lao  39°25-08%  lonf  86°00  18"  in  S  .*  X/h  &ec  5,  7  11  N  .  R  5  E.  ,  on 
left  bank,  on  upstream  aide  of  highway  bridge  half  a  mile  southwest  cf  Amity, 
2  milea  upstream  from  movth_t  and   5  milee  northwest  of  Hdinburg- 

Drainage  area   -109  eq  rai 

Caere  --Nonrecordinf  pi«  Dec.  7?  19k?.,  to  June  29,  1955j  recording  gape  thereafter 
Datum  cf  gage   is  670  2'  ft  above  mean  sea  level,  datum  of  L.929. 

Stage -discharge   relation     -Defined   by  current -meter  measurements   below  7.000  cfs 
and  by  contracted -opening  meacurement  at  10;700  cfs 


Fear.  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Tear 

Date 

Ca?e 
Height 

Discharge 
cfs 

Y.'ater 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

1943 

Kar 

19, 

1943 

ic  to 

3,700 

1952 

Jan- 

27, 

1952 

13  4 

10.700 

1944 

Apr 

11. 

1944 

11  00 

4,290 

1953 

!lar. 

4, 

1953 

8  37 

2,080 

1945 

Mar 

6t 

1945 

11  00 

4,290 

1954 

Jan. 

27, 

1954 

3  27 

443 

1946 

Hay 

16, 

1946 

9  0 

2,510 

1955 

!!?.y 

28, 

1955 

6  2 

1,110 

1947 

June 

2, 

1947 

11-12 

4,3C0 

1956 

I.'ov  T 

16, 

1955 

12.20 

7,790 

1948 

Mar 

27, 

1948 

7-63 

1,650 

1957 

July- 

5, 

1957 

11.62 

6,510 

1949 

Jan 

5, 

1949 

11  9 

.5,190 

1958 

June 

11, 

1958 

4,350 

1950 

Jan 

4, 

1950 

10  8 

L,090 

1959 

Jan 

21. 

1959 

6.270 

1951 

Jan 

15, 

1-51 

9  15 

2,730 

YOUNG        CREEK        NEAR         EDINBURG     , 

Return  Feriod  (years) 
5  10  25 


0o.T 

i 


1  76        HO     "  •     9C 
>  J  Probability  -  percent 


-T-- 


J-L^. 


'1 7 


Reduced  Varlate       (Y) 


(0)     Blut.  River  at  Carthage,   Ind, 

Looat|™  ~Ut  -5°''6    '   long    i-°3/-   ■•   in  3ec     18>  T     "  !:  :  "     ?  S.a  on  right  bank 
500  ft  upstrc.m  lrom  highway  bridge,  hair  a  mile  went  or  Carthage,  and  2  i 
milas  downstream  from  Three  (tile  Creek 

Drainage  area . •   18/  sq  mi 

Gage  -Nora-ecording  gage  Oct  11,  1950.  to  July  13,  1951;  recording  gage  thereafter 
Frio;-  »o  July  19,  1951,  at  bridge  500  ft  downstream  Datum  of  gage  is  859  33  ft 
above  laean  sea  levels  datum  it  1929, 

Stage -discharge  relation. —Defined  by  current  meter  measurements 
Flood  stage     -7  ft 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfa 

191.9 

Jan       5,  1949 

10  6 

5,750 

1955 

Jan,     6,  1955 

5  97 

1,290 

1951 

Feb     21.    1951 

11,2 

6,650 

1956 

Nov     16,  1955 

11.52 

5,800 

1952 

Jan,   27,   1952 

11  02 

6,350 

1957 

June  18,  1957 

9  77 

3,900 

1953 

Mar       4,   1953 

9  17 

3,580 

1958 

Jure  14,  1958 

7,020 

1954 

Apr,     5,  1954 

10  02 

4,850 

1959 

Jan     21,  1959 

8,340 

BLUE       RIVER       AT       CARTHAGE      ,    INDIANA 

Return  Period  (years) 
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Reduced  Variate       CY) 


(9)     Sand  Creek  near  3rexer3ville,  Ind. 

Location  --Lat  39°05'05",  lone  B5°39(30»,   in  N>:{  3ec     5,  T.   7  N   ,  h     8  E  ,  on  left 

bank  at  downstream  side  of  count}    highway  bridge,   2^  miles  west  of  Brewer3ville B 
and  5  2  miles  upstream  from  Bear  Creek. 

Drainage  area     -156  aq  mi;  153  sq  ml   j>rior  to  Oct.   6,  1952,. 

Gage.-Nonrecording  gage  Feb,   11,  19i.fi,   to  Oct..    5S   1952,  at  bridge  1,7  miles  upstream 
at  datum  approximately  8  ft  higher.      Recording  gage  since  Oct,   6,  1952,  at 
present  site       Altitude  of  present  gage  is  630  ft  (by  altimeter)  =, 

Stage  --discharge  relation. — Defined  by  current-meter  measurements  at  former  site  and 
by  gage-height  relationship  with  former  site  at  present  location 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

Kater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1948 

Mar 

Z7, 

1943 

17  5 

9,980 

1954 

June  1, 

1954 

5  75 

1,240 

WW 

Jan 

5» 

1949 

19  0 

l:>,100 

1955 

Feb.  27, 

1955 

11  42 

4  300 

1950 

Jan 

4i 

1950 

19  2 

12,400 

1956 

May  28, 

1956 

15  45 

7,560 

1951 

Nov 

20, 

1950 

18,4 

11,100 

1957 

Apr   4, 

1957 

16.33 

8, ,480 

1952 

Jan 

26, 

1952 

13  4 

5,780 

1958 

July  22, 

1958 

7,150 

1953 

Mar 

4, 

1953 

10.19 

3,460 

1959 

Jan.  21, 

1959 

19  ,,900 

CREEK       NEAR        BREWERSVILLE 

Return  Feriod  (years) 
5  10  25 


Reduced  Variate 


(10)     North   for;:  S:.l:    "reek  near  9.-lno:.t,    1::6 

vocation      -Lilt  39"C9  CC" ,   long  -.6°/;  1.'  ,   in  IiV  ?./;    scr     5     T     S  :: . ,    \     2  3  ,   on 
right  ban!.  15  ft  covnstream  fro*.-,  oridge  on  jv-.te  Highway  46     100  ft  u-:>3tream 
fro-  Schooner  '.rork,  0-7  mile  northeast   jf  Reinrant,    6  1/2  -liles  upstream  from 
Brucmett  Creek .   and  20  miles  ujjstraa-  fron  mouth 

Drainage   irea      -120  so,  mi,   incljde:.  that  of  Schooner  Creek 

Sage    -    ,-wnrc-cor-iing  gage  Apr,   4,   194(-,   to  Oct     8,    1951;   recording  gage  thereafter 
Altitude  of  gape  is   546  ft  (fro:;  topographic  nap) 

.Stage -discharge  relation. --Defined  by  current  meter  measurerants  below  9,800  cfs 

Discharge  adjusted  for  rate  of  change  of  stage  above  7  ft       Only  annual  maximums 
adjusted  prior  to   installation  of  recording  gage 


Peak  Stares  and   Instantaneous  Annual  peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Uater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1913 

March 

1913 

25  7 

1953 

Mar   4, 

1953 

17- 76 

2,780 

1946 

May  16, 

1946 

20  3 

5,910 

1954 

Jan.  27, 

1954 

9  38 

825 

1947 

June  2, 

1947 

21  2 

10,100 

1955 

Mar  21, 

1955 

15  91 

2,220 

1948 

Mar  27, 

1948 

18.0 

3,010 

1956 

Kay  28, 

1956 

18  12 

3,030 

1949 

Jan   5, 

1949 

22  2 

13.300 

1957 

Apr,   4, 

1957 

19.92 

6,340 

1950 

Jan   4, 

1950 

21  7 

11,600 

1958 

June  14 1 

1958 

5,920 

1951 

Feb  21, 

1951 

19  53 

5,100 

1959 

Jan  21, 

1959 

12,000 

1952 

Hay  24, 

1952 

22  55 

15,200 

NORTH      FORK      SALT      CREEK      NEAR      BELMONT  ,    INDIANA 


5     10      20     30    43    W  60      76  '     'HO 

2  }  Probability  -  percent 

,      I     ,  ,  ,  |  ,  it  ,      ,      ?         ,  , 

Hedjced  Variate  C Y3 


~"*7 


(11)  Patoka  riiver  at  /asper.  Ind 

Location  --La t  38°24'49",  long  86°52  ;,6",  in  SEj  sec  20,  T  1  S  ,  R  4  w  ,  on  left 
bank,  0  3  mile  upstream  from  unneraed  outlet  of  Jasper  Lake,  1  0  mile  downstream 
from  Coon  Seltz  bridge,  1.2  miles  downstream  from  Beaver  Creek,  and  3  3  miles 
northeast  of  Jasper- 

Drainage  area  --257  sq  mi;  270  sq  mi  at  former  site 

Gage  --Nonrecording  gage  Nov  20,  1947,  to  Sept.  17,  1956;  recording  gage  thereafter 
Prior  to  Sept  18,  1956,  at  site  5  6  miles  downstream  at  datum  0  34  ft  lower; 
datum  of  present  gape  is  446  19  ft  above  mean  sea  level,  datum  of  1929 

Stage-discharge  relation. — Defined  by  current  meter  measureaients  below  5,000  cfs  at 
former  site  and  below  1,100  c's  for  present  site. 

Flood  stage  —14  ft;  9  ft  at  former  site. 

Historical  data. --Flood  of  March  1913  is  maximum  stage  known.  Maximum  stage  at 

present  site  for  period  1925-57,  20  ft  in  1925  (information  from  local  resident). 

Remarks.  — Flow  slightly  regulated  by  Beaver  Creek  Tieeervoir,  whose  outlet  enters  the 
Patoka  River  1.2  miles  upstream  from  the  gage;  peak  discharges  not  materially 
affected 

Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Kater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1913 

March 

1913 

15  9 

16,000 

1953 

Mar 

8,  1953 

7  90 

1,640 

1937 

January 

1937 

14  8 

12,100 

1954 

Mar 

2,  1954 

- 

950 

1948 

Apr  15, 

1948 

11  57 

4,920 

1955 

Mar 

3,  1955 

9  80 

2,940 

1949 

Jan.  28, 

1949 

11.13 

4,220 

1956 

Feb 

29,  1956 

9.98 

3,100 

1950 

Jan.  7, 

1950 

12.37 

6,300 

1957 

i  ly 

2'.s   1957 

17  87 

6,900 

1951 

Mar  21, 

1951 

11.46 

4,760 

1958 

Dec 

22,  1957 

4,250 

1952 

Mar.  14, 

1952 

10  78 

3,880 

1959 

Jan 

24,  1959 

9,150 

POTOKA      RIVER       AT      JASPER     ,         INDIANA 

Return  Feriod  (yyars) 
2  5  10  25  50 


Reduced  Varlate        <Y) 


(12)     Busseron  Creek  near  Carlisle,   Ind. 

Location — Lat  38°58  30",  long  87°25>35"     in  wj   aee     17    t    f.  »       u     a  ,, 

ll^and  ffS  *V**     "  ^  ^  *  ^-  nVrtn^stTL- 
iisxe,  and  o  3/i,  miles  upstream  from  mouth. 

Drainage  area, — 228  sq  ml. 

0&gB^lZTtT^egdT2,'f6^fl5^^l,3'   ^  KOV\7'   19,5°:   "cortln«  ™°  thereafter 
lXu!  bJnch  martf  ^  "  ""'  ""  "a  l6Vel   (Stat"  ^^  ".prt-ent  of 

St*B"^2!!J"'hr1*Ki0?,~D?f1™d  by  =urrent-«ter  measurements  below  4,500  cfe  and 
extended  above  by  logarithmic  plotting,  »«.»~  >-io  «™ 

Flood  stage. --12  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Water 

Tear 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Tear 

Date 

Cage 
Peight 

Discharge 
cfs 

1944 

Apr 

12, 

1944 

16  96 

4,700 

1952 

Mar 

11. 

1952 

16.17 

4,070 

1945 

Apr 

2» 

1945 

17,60 

5,500 

1953 

Kar 

4, 

1953 

16.0/, 

3,890 

.1946 

Hay 

20, 

1946 

14  90 

2,900 

1954 

Aug, 

4, 

1954 

6.31 

430 

1947 

June 

2, 

1947 

14-60 

2,720 

1955 

Apr 

13. 

1955 

13-13 

2,040 

1948 

Jan. 

3, 

1948 

15.15 

3,100 

1956 

June 

22, 

1956 

16.12 

3,980 

1949 

Jan. 

20, 

1949 

16.3 

4,200 

1957 

Hay 

23, 

1957 

17  61 

5,200 

1950 

Jan. 

5. 

1950 

20.05 

8,800 

1958 

Dec 

21, 

1957 

5,400 

1951 

Feb. 

21, 

1951 

14.75 

2,900 

1959 

Jan. 

22, 

1959 

3,100 

BUSSERON       CREEK       NEAR       CARLISLE     ,   INDIANA 

Return  Period  (years) 
2  5  10  25  50 


°o.T 


5    IB    20   Jo'iW  50  66    76      85'  90~ 

J  $  Probability  -  percent 

1 1 . ' . . ! .  1 . 1 1 1  u  ijs  1 1 1 !  a  1 1 1 1 1 1 . 


9"9"7o         9979 


~=r 


V 


-Y 


>""'""o J 


Reduced  VarUte  (Y) 


(131     tiast  Fork  Whitewater  River  at  Richmond,  Ind. 

Location. --Lat  39°48'24",  long  84°54<25",   in  SE  1/4  sec.   7,  T.  13  K.,  H„  !».,  on 
left  bank  50  ft  downstream  from  highway  bridge,  three-quarters  of  a  nolo  soutn 
of  Richmond,  and  2  miles  upstream  from  Short  Creek. 

Drainage  area. — 123  sq  mi. 

Gage.-Nonrecordlng  gage  Apr.  27,  1949,  to  July  26,  1949;  recording  gage  thereafter. 

Datum  of  gage   is  854  01  ft  above  -near,  sea  level,  datum  of  1929  (levels  by  Indiana 
Flood  Control  and  Water  Resources  Ccnmission) , 

Stage-discharge  relatlon.—Definad  by  jurrent-aieter  measurements  below  5,100  cfs  and 
by  slope-area  measurement  at  13,500  cfs. 

Flood  stage. — 10  ft. 

Hiatoriral  data. --Flood  of  September  1366  was  reported  by  the  Indianapjlis  Journal  to 

be  higher  than  ever  before  known.      Flood  of  torch  1913  is  the  maximum  stage  known 

according  to  information  by  local  residents. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 

Water 
Tear 

Date 

Gage 
Height 

Dischargs 
cfs 

Y.'ater 

Year        Date 

Gage 
Height 

Discharge 
cfs 

1913 

March    1913 

15  0 

- 

1955     Feb  21,  1955 

6  07 

2,540 

1950 

Jan  15,  1950 

12  49 

13,500 

1956    Nov  16,  1955 

10,70 

8,200 

1951 

Nov.  20,  1950 

10.82 

9,250 

1957     June  28,  1957 

10-54 

7,800 

1952 

Jan.  26,  1952 

10.66 

9,250 

1958     Aug,  2,  1958 

10,400 

1953 

Hay  22,  1953 

6.53 

3,160 

1959     Jan.  21,  1959 

14,100 

1954 

Mar.  30,  1954 

3.86 

1,160 

EAST      FORK      WHITE      RIVER        AT      RICHMOND  ,   INDIANA 
Return  Period  (years) 
2  5  10  25  50  100 


L    < 


Seduced  Varlate  (Y) 


i'lh)     Silver  >eek  r.ezr  Liel^ersbi-rR 

-        ♦•  ,.,    «°22il5«     long  05°«;J5'',   inSWf   Let  68,  Clark  Military  r.rant,  on 

Location  - -U-t  38  22  x_    ,   lor,        3  __,  ^a       Q       ^  dovmstrea:B  fr0It 

mouth 
Dra-n3 '3  area.— 188  sq     mi. 
Gaff  .-- Wire-weight  [■«?«  read  twice  daily. 


Peak 

Stares  a 

nd   Li  itar.iar 

eiu3    \nn 

uil  reak  Ulscrarp 

Water 
Year 

Hate 

Gage 
Height 

Disc  iarf  s 
c.'s 

Water 
Year 

Date 

Gape 
Feip.ht 

Discharge 
c:*9 

1955 
IS  56 
1957 

Feb     22, 
Feb.     2, 
May     23, 

1955 
1956 
1957 

L,221 
5,250 
6,250 

1959 
1 

Rev    19,   1957 
Jan     22,   1959 

5,080 
19,600 

SILVER      CREEK      NEAR      SELLERSBURG    ,    INDIANA 
Return  Period  (years) 
2  5  10  25 


oo        looo 


9979 


'     !■  -'""I" I 


Reduced  Varlate        CYJ 


0  5)     Big  Ir.dia-.i  Creek    r>.i:   Onrytion,  Ind 

Location. —La  t  38°L6'35",  long  86°C6':)5",  in  S3  1/4  sec.  6S  T  3  S.,  R  4  2.,  on 
upstream  sine  of  bridge  on  State  Highway  335,  0.6  mile  upstream  from  iaoooon 
Branch  ard  4  1/2  nilee  north  o"  Corydon. 

Drainage  area. — 12J    in  mi. 

Gaje.     Ilonrecordlnj  gage  Oct.    i.2,  194',  to  Doc     8,  1948;  racer  ling  Rage  thereafter, 
Datum  of  gage  is  577.12  ft  above  i«an  sea  level,  datum  of  1929. 

Stage  -d-scharge  relation. —  Defined  :>y  current-ceter  measurements  below  6,600  cfs  and 
emended  abcv;  by  logarithmic  plotting. 

Historical  data. --"lend  of  Mar.,  lc),  1°43,   is  the  na.-d.nuni  knoi.n  at  Corydon  aince 
beginning  of  knowledge  in  18.".5„ 


Peik  Soages  and   Insts  itaneiua  Annual  Peal  Discnarge 


Water 
Year 

Jats 

Gage 
Height 

D.'.9uh*.r  ;s 

CfL. 

V'ate: 
'■."ear 

Date 

Gage 
Height 

Discharge 
cfs 

19',? 

Bar  .   19, 

TQI   ■ 

22,4 

17,  000 

10;2 

!-«ar„  11, 

1952 

18  ..76 

10.,  400 

1944 

.'.pr     11, 

19.4.* 

14  0 

4,8(0 

1953 

Mar.     3, 

1953 

14  57 

5  400 

1945 

!fer.     6, 

194  5 

19,2 

11.0C0 

1954 

Sept ,20, 

1954 

14  04 

4,800 

1946 

Feb     14, 

194': 

15.0 

5,9J0 

1955 

liar,,   15, 

1955 

12,98 

3,900 

1947 

,j...     2j 

19V, 

14.3 

5.160 

1956 

Fob.    25, 

1956 

16  25 

7,180 

194? 

/or..  12, 

194.; 

19.3 

11.100 

1957 

Feb.    10, 

1957 

16.32 

7,300 

1949 

June  15, 

1949 

16.20 

7,180 

1958 

Kov.  19, 

1957 

6,590 

195C 

Feb .     9 , 

1950 

16  77 

7.-900 

1959 

Jan.   21, 

1959 

23,800 

;  1951 

•■  =  '■      '4, 

' 

If  -1 

6  a;o 

— 

.,     

i 

BIG      INDIANA 


1.01 


CREEK      NEAR       LORYDON 

Return  Period  (years) 
5  10  25 


Reduced  Varlate         CYJ 


i!6)     Kankakee  hivei-  r.ea:    North   Liberty,    Iivi 

Location— Lai.  / 1°  ';).',  long  :;i:2;'  JO",  rn  line  between  seci  11  and  23,  T.  36  I.  , 
K.  1  V  ,  m  loft  bank  at  cowr  stream  siae  rf  brirtre  or.  3t .  Joseph  County  highway 
nare.1  "New  Road."  i  miles  northwest  of  ::orth  Liberty. 

r-.'jtie  area,— 152  sq  mi. 

Gage,— Konrecordisr,  gaffe  Jan.   12,  1951,  tc    Jur.e  25,   1' 56;   reco:-dir..-i  gage  thereafter. 
Datun  of    »a?e  is  680. 0^  ft  above   nean  sea  level,  datum  of  19L9  (levels  by  Indiana 
Flood  Coir.rol  and  ','ater  llecource.-;  CoEi.ission). 

Stare  oischarrii  relation. •  -fl3lat  icn  affected  by  varying  arrun  ...  of  backwater  caused 

by  return  flow  from  overbank   storage.      Frequent  current  laeter  iraasureaients  neces- 
sary to  rtafine  relationship  durin;  this  period. 


Peak  ota?es  anc    li  star  tar.ecu?   Annual  i  eak  Dlscharj 


Water 
Year 

Date 

Gage 
Height 

Dischs  r,  e 
cfa 

1  'ater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1951 

Kay    11 , 

1951 

6.15 

520 

lV.-.b 

a-r     30,  1956 

6,92 

660 

195? 

Nov.  lit, 

195- 

f.97 

680 

1957 

\cr.   .57,  1957 

0.90 

660 

1953 

Mar.  16, 

1953 

4.i" 

260 

1958 

0;c.   20,  i>57 

560 

195a 

Apr     26, 

195's 

800 

1959 

"ar     27,  1559 

560 

1955 

net     10, 

195fc 

R.6i 

6R6 

KANKAKEE        RIVER       NEAR        NORTH       LIBERTY 
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Reduced  Vari.it.- 


(17)     Singleton  titch 


Schneider,   Ind, 


Location, --Lat  41°12'44",  long  87<:26'.J.",   on  line   oetveen  1"S  1/4  sec.  21  and  NVf 

l/lt  sec.  22,  T,  32  K.,  U.  9  I.'.,  on  left  bank  15  ft  upstream  from  bridge  on  U.  S 
Highway  41,  half  a  nile  upstream  fr:>m  3ruce  ditch,  1  1/2  miles  dovnstrean  from 
Cedar  Creak,  and  1  2/3  miles  nor.h  of  Schneider, 

Drainage  area. — 122  sq  mi. 

Cage  —  I onrecording  gage  July  23,  194!,  to  Aug.  13,  1951;  recording  ra-e  thereafter, 
Prior  to  Oct.   1,  1949,   at  datum  ?.00  ft  higher.      Datum  of  present  gage  is  623  67 
ft  above  niean  sea  level,   datum  oS  1929. 

Stage-discharge  relation. — Defined  by  current-ireter  measurements,      Dredging  in  1950 
and  subsequent  floods  and  channel-  deterioration  have  materially  affected  the 
stage-discharge  relation, 


Peak  Stages  and  Instantaneous  Ann-al  Peak  Discharge 


Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
"fear 

Date 

Gage 
Height 

Discharge 
cfe 

1949 

Feb  15,  1949 

55C 

1955 

Oct, 

11,  1954 

10,10 

953 

1950 

Apr  10,  1950 

- 

1,100 

1956 

Feb, 

25,  1956 

9  62 

338 

1951 

Feb  19,  1951 

8, 50 

841 

1957 

Apr 

28,  1957 

10  27 

979 

1952 

June  15,  1952 

9.82 

1,010 

1953 

714 

1953 

3,39 

312 

1959 

?eo0 

14,  1959 

992 

1954 

Xar-  25,  1954 

9.04 

SIC 

SINGLETON      DITCH       AT      SCHNEIDER    ,    INDIANA 

Return  Feriod  (years) 
2  5  10  25  50  100 


Reduced  VaxUte  (Y) 


(18)     Deep  Hiv-r  _t  Lai      George  Outlet  at  Kobart,    Ind. 

Location-—  Lat  41°32'i0".   ion?  87°li'  ;=."     i„  »n,  i /, 

left  ban*  at  upstream ^d.  If     i 'L^, l^oT 'fti2'  ?■   36  ■■•   »•    7  *.  - 
and  400  ft  downstream  from  Lake   -eorse  Dam!  upstream  fron  Duck  Creek, 

Drainage  area.  — 12=  eri  mi. 

present  datum.      DatuTof  present 'gte   is'  88^  £  SET  ^  "  "^  "»" 
of  1929  (levels  by  Indian^  Flood  S&S  S3£  H^urcTc^sfe  ^ 

Stage-discharge  relation. --Defined  bv  current  »»(■=. 

uoiajiea  oy  current-meter  measurements  below  3,300  cfs. 


Water 
Year 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


Gage       Discharge 
Height  cfs 


1947 

Apr, 

6j 

1947 

5  41 

2,410 

1948 

Hay 

11< 

1948 

5,86 

2,740 

1949 

Feb 

14, 

1949 

3,50 

620 

1950 

Dec. 

22, 

1949 

5  35 

2,390 

1951 

May 

Hj 

1951 

4,52 

1,440 

1952 

Nov. 

u, 

1951 

4.41 

1,340 

1953 

iter. 

16, 

1953 

3.86 

912 

Water 
Year 


Date 


Gage         Discharge 
Height  cfs 


1954 
1955 
1956 
1957 
1958 
1959 


Kar-   26,  1954  4.55 

Oct.  11,  1954  7.68 

Kay     LL,  1956  11,15 

July  14,  1957  12.3J 
June  10,  1958 
July  24,  1959 


1,440 
3,880 
1,320 
1,650 
720 
1,970 


DEEP      RIVER       AT       LAKE      GEORGE      OUTLET     AT     H08ART ,     IND 

Return  Period  (years) 
2  5  10  25 50  100  200 


5X> 


Reduced  Varlate        CY) 


(19)     Pigeon  Creek  at  Hogback  Lake  Outlet,  near  Angola,   Ind. 

Location..-- Lat  41°37'24",  long  S:'°05 ' 4.V" ,   in  t:E  1/4  NW  1/4  sec.  36,  T.   37  N.,   H.  12  E., 
on  right  bank  200  ft  north  of  lake  outlet,   2  miles  southeast  of  Flint,  and  5.1 
miles  west  of  Angola. 

Drainage  area. — 102  sq  mij  105  sq  mi  prior  to  October  1947. 

Cage. — Conrscording  gape  Oct.   16,  1915,  to  Aug.  3,  1953;  recording  page  thereafter. 

prior  to  Oct.   I,  1947,  at  site  1  1/2  miles  downstream  at  different  datum.     Oct.  1, 
1947,  to  Aug.  3,  1953,  at  site  600  ft  downstream  at  present  datura.     Datum  of 
present  gage  is  940,00  ft  above  mean  sea  level,  datura  of  1929. 

Stage-discharge  relation. — Defined  by  current-meter  measurements  below  240  cfs  at 
former  site  and  by  current-meter  measurements  at  present  site. 


Peak  Stages  and   Instantaneous  Annual  peak  Discharge 


Water 
Year 

Date 

Gage 
Height 

Discharge 

CfB 

Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

1946 

Feb.  19, 

1946 

- 

220 

1953 

tor.  19,  1953 

9,30 

122 

1947 

Apr.  24, 

1947 

10.71 

458 

1954 

Kar-  30,  1954 

11,31 

317 

1943 

Mar .  2, 

1948 

11.79 

355 

1955 

Oct.  17,  1954 

11,54 

339 

1949 

Feb.  19, 

1949 

11.93 

366 

1956 

Kay   4,  1956 

13.39 

548 

1950 

Apr.  8, 

1950 

14.95 

744 

1957 

Apr.  14,  1957 

11.29 

317 

1951 

Feb.  24, 

1951 

12.50 

448 

1958 

Sept. 21,  1957 

274 

1952 

Jan-  21, 

1952 

11.85 

370 

1959 

Feb ,17-19, 1959 

442 

PIGEON        CREEK       AT       HAGBACK      LAKE      OUTLET       NEAR       ANGOLA      ,        INDIANA 

Return  Period  (years ) 
1.01  2  5  10  25  50  100  200 


5    """"t f 


Raduced  VirUt«        IYJ 


(2C)     Laughery  Creux  r.e:-r  Fsritrs  .-tetreat,   Ind. 

Location.— Lat  38c57'05",   lcng  85°Oi-22",   in  eec.   2,  T.  1  I!.,  K.  3  H.,  on  right  bank 
2  miles  southeast  of  Farmers  Retreat  and  3  5 A  raJ.es  doxnstreara  from  Bear  Creek. 

Drainage  area, — 218  sq  ni. 

Ga;;e.— tlonrecording  ga;-e  Oct.  3,  1910,  to  Apr,   15.  19U;  recording  gage  thereafter. 
Altitude  of  gape  is  526  ft  (by  baroreter) , 

Stage-discharge  relation. — Defined  by  current -meter  nH5sureEent3  below  11,000  efa. 

Flood  stage.-  13  ft. 

Historical  data.— Flood  of  1397  reached  a  stage  of  about  18  feet  and  is  the  highest 
known  flood,   from  information  by  local  residents. 


Peak  Stares  and  Instantaneous  Annual  Peak  Discharge 


Water 
Tear 

Date 

Gage 
Height 

Discharge 

Cf3 

V.'ater 
tear 

Dlte 

Gage 
Height 

Discharge 
cfs 

1941 

June  9 

or 

10,1911 

9-62 

3,960 

1951 

Jan. 

3, 

1951 

13.50 

9,660 

1912 

Apr, 

9, 

1912 

13.91 

10,800 

1952 

Mar 

10, 

1952 

13  16 

9,660 

1913 

liar, 

19, 

1913 

11  50 

12,900 

1953 

Kay 

17, 

1953 

9  U 

3,960 

1911 

Apr 

U, 

1911 

13.16 

8,880 

1951 

Hay 

3, 

1951 

3  99 

610 

1915 

Mar, 

6. 

1915 

15:51 

17,000 

1955 

Mar, 

21, 

1955 

13-88 

10,800 

1916 

Feb 

13, 

1916  * 

12:  TP 

7,980 

1°56 

Hay 

28, 

1956 

14-45 

12,500 

1947 

May 

25, 

1917 

11.62 

13,300 

1957 

July 

5, 

195' 

16.15 

20,200 

1948 

Apr, 

12, 

1918 

13  01 

8,110 

1958 

July 

22, 

1958 

17.000 

1919 

Jan. 

21. 

1919 

15,23 

15,700 

1959 

Jan, 

21, 

1959 

17,800 

1950 

Feb 

3, 

1950 

If  03 

l  ,eo~ 

LAUGHERY       CREEK       NEAR        FARMARS       RETREAT  ,    INDIANA 

Return  Feriod   (years) 
25 


5oo         iooo 


.4O0C3- 


°0.T 


Hedjced  Varlate  lY) 


APPENDIX      B 
I.     Instantaneous  hydrograph 


II 


Dimensionless  instantaneous  hydrograph 
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APPENDIX  B-l 
Instantaneous  Hydrograph 


qtp           (n-l)       „n 
AR       "     T(n)         r 

-1     -A 

e 

X  = 

(n-l)  4 

P 

t/t 

'   P 

1.4 

1.6 

n 
1.8 

2.0 

2.2 

2.5 

0.1 

0.119 

0.099 

0.104 

0.090 

0.076 

0.056 

0.2 

0.152 

0.166 

O.I69 

0.164 

0.150 

0.131 

0.5 

0.170 

0.200 

0.215 

0.222 

0.221 

0.216 

o.4 

0.186 

0.225 

0.251 

0.268 

0.280 

0.298 

0.5 

0.195 

0.241 

0.277 

0.504 

0.526 

0.346 

0.6 

0.201 

0.254 

O.296 

0.529 

0.556 

0.392 

0.7 

0.205 

0.262 

0.509 

0.548 

0.584 

0.426 

0.8 

0.208 

0.267 

0.517 

0.559 

0.398 

0.447 

0.9 

0.209 

0.270 

0.522 

0.566 

0.4o6 

0.46o 

1.0 

0.210 

0.272 

0.525 

0.568 

0.409 

0.463 

1.1 

0.209 

0.270 

0.522 

0.366 

0.407 

0.460 

1.2 

0.208 

0.268 

0.519 

0.564 

0.400 

0.450 

1.5 

0.207 

0.265 

0.515 

0.555 

0.591 

O.457 

1.* 

0.205 

0.261 

0.507 

0.546 

0.378 

0.418 

1.5 

0.202 

0.256 

O.299 

0.554 

0.366 

0.400 

1.6 

0.199 

0.251 

0.292 

0.525 

0.350 

0.580 

1.7 

0.196 

0.245 

0.285 

0.511 

0.533 

0.358 

1.8 

0.195 

0.259 

0.275 

O.297 

0.317 

0.336 

1.9 

O.I89 

0.252 

0.261 

0.285 

O.299 

0.316 

2.0 

0.186 

0.225 

0.255 

0.270 

0.273 

O.292 

2.2 

0.178 

0.211 

0.252 

0.244 

0.248 

O.250 

2.4 

0.170 

0.198 

0.215 

0.218 

0.218 

O.209 

2.6 

0.162 

0.184 

0.195 

0.192 

0.188 

0.174 

2.8 

0.154 

0.170 

0.174 

0.172 

0.176 

0.146 

5.0 

0.146 

0.158 

0.157 

0.150 

0.139 

0.119 

5.5 

0.128 

0.128 

0.120 

0.105 

0.092 

0.071 

4.0 

0.110 

0.105 

O.O89 

0.072 

0.059 

0.042 

4.5 

0.094 

0.081 

O.065 

0.046 

0.057 

0.024 

5.0 

0.080 

O.065 

0.047 

0.054 

O.025 

0.016 

.216 


APPENDIX  B-l 


qt   /AR 
p' 


t/t 

'  p 

3.0 

4.0 

n5.0 

6.0 

7.0 

10.0 

0.1 

0.032 

0.010 

0.003 

0.001 

0.002 

0.000 

0.2 

0.107 

0.060 

0.031 

0.015 

0.0075 

0.008 

0.5 

0.198 

0.148 

0.104 

0.070 

0.0469 

0.013 

o.i* 

0.287 

0.260 

0.222 

0.180 

0.1445 

0.036 

0.5 

0.368 

0.376 

0.560 

0.535 

0.505 

0.208 

0.6 

0.455 

0.482 

0.505 

0.504 

O.494 

0.454 

0.7 

0.484 

0.565 

0.625 

0.661 

0.686 

0.712 

0.8 

0.517 

0.626 

0.715 

0.782 

0.857 

0.962 

0.9 

0.535 

0.662 

0.765 

0.853 

0.930 

1.130 

1.0 

0.540 

0.672 

0.782 

0.878 

0.972 

1.180 

1.1 

0.537 

0.662 

0.770 

0.857 

0.963 

1.140 

1.2 

0.522 

O.658 

0.730 

O.803 

0.866 

1.010 

1.3 

0.502 

0.599 

0.675 

0.726 

O.767 

0.841 

1.4 

0.477 

0.556 

0.608 

O.639 

0.658 

0.670 

1.5 

0.448 

0.506 

0.536 

0.547 

0.547 

0.511 

1.6 

0.418 

0.455 

0.466 

0.458 

0.441 

1.7 

0.586 

0.4o4 

0.397 

0.576 

0.549 

1.8 

0.554 

0.556 

0.336 

0.505 

0.269 

1.9 

0.525 

0.510 

O.279 

0.250 

0.205 

2.0 

O.295 

0.268 

0.230 

O.I89 

0.155 

0.075 

2.2 

0.238 

0.195 

0.151 

0.112 

0.0815 

2.4 

0.189 

0.159 

0.096 

0.064 

0.0414 

2.6 

0.149 

0.097 

0.034 

0.0201 

2.8 

0.116 

O.067 

0.019 

0.0095 

5.0 

0.090 

0.045 

0.021 

0.010 

0.0045 

0.017 

5.5 

0.046 

0.030 

4.0 

0.022 

0.012 

4.5 

0.010 

5.0 

0.004 
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Appendix  B-2 
Dimensionless  Instantaneous  Hydrograph 


t/t 

l.k 

1.6 

n 

1.8 

2.0 

2.2 

2.5 

0.1 

56.6 

36.  4 

32.2 

24.4 

18.6 

12.1 

0.2 

72.  U 

61.4 

52.4 

44.8 

36.7  • 

28.3 

0.3 

81.0 

73-5 

66.6 

60.4 

54.1 

46.7 

0.4 

88.6 

82.7 

77.7 

72.8 

68.5 

64.4 

0.5 

92.9 

89.6 

85.8 

82.5 

79.7 

74.8 

0.6 

95.7 

93.  4 

91.7 

89.4 

87.O 

84.6 

0.7 

97-6 

96.  4 

95.7 

9^.5 

94.0 

92.0 

0.8 

99.0 

98.I 

98.2 

97.5 

97.4 

96.5 

0.9 

99.5 

99-3 

99-6 

99-5 

99.3 

99.3 

1.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

1.1 

99.5 

99.3 

99.6 

99.5 

99.5 

99.3 

1.2 

99.0 

98.5 

98.8 

99.0 

97.9 

97.1 

1.3 

98.6 

97.5 

97.0 

96.5 

95.6 

94.4 

l.k 

97.6 

96.O 

95.1 

94.0 

92.5 

90.3 

1-5 

96.2 

94.1 

92.6 

90.7 

89.5 

86.4 

1.6 

9^-7 

92.3 

90.5 

87.8 

85.5 

82.1 

1.7 

93.4 

90.0 

87.6 

84.  5 

81.5 

77.4 

1.8 

91.9 

89.9 

84.5 

80.7 

77.6 

72.6 

1.9 

90.0 

85.4 

80.9 

77.5 

73.1 

68.3 

2.0 

88.6 

82.8 

78.4 

73.4 

66.8 

63.1 

2.2 

84.8 

77-5 

71.9 

66.4 

60.6 

54.0 

2.4 

81.0 

72.9 

66.0 

59-3 

53.4 

45.2 

2.6 

77.2 

67.6 

59.8 

52.2 

46.0 

37.6 

2.8 

73.4 

62.5 

53.9 

46.8 

43.0 

31.6 

3.0 

69.6 

58.1 

48.6 

4o.7 

34.0 

25.7 

3.5 

61.0 

U7.0 

37-2 

28.5 

22.4 

15.3 

4.o 

52.4 

37.9 

27.6 

19.6 

14.4 

9.1 

4.5 

44.8 

29.8 

20.2 

12.5 

9.1 

5.2 

5.0 

38.1 

23.9 

14.6 

9.2 

5.6 

3.5 
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APPENDIX  B-2 
Dimensionleas  Instantaneous  Hydrograph 


t/t 

1  p 

3.0 

4.o 

n 
5.0 

6.0 

7-0 

10.0 

0.1 

5.9 

1.5 

0.4 

0.1 

0.0 

0.0 

0.2 

19.8 

8.9 

4.0 

1.7 

0.8 

0.7 

0.3 

36.6 

22.0 

13.3 

8.0 

4.8 

1.1 

o.4 

53.2 

38.7 

28.4 

20.5 

14. 9 

3.1 

0.5 

68.1 

56.0 

46.0 

38.9 

31.2 

17.7 

0.6 

80.2 

71.8 

64.4 

57.4 

50.9 

36.8 

0.7 

89.6 

84.1 

80.0 

75.3 

70.6 

60.4 

0.8 

95.6 

93-1 

91.1 

89.O 

86.3 

81.5 

0.9 

99.1 

98.6 

98.0 

97.1 

95.7 

95.8 

1.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

1.1 

99.  ^ 

98.6 

98.5 

97.6 

99.1 

96.6 

1.2 

96.6 

95.0 

93.4 

91.5 

89.1 

85.5 

1.3 

93.0 

89.1 

86.4 

82.7 

79.0 

71.3 

1.4 

88.4 

82.9 

77.7 

72.8 

67.7 

56.8 

1.5 

83.O 

75.4 

68.5 

62.2 

56.3 

43.3 

1.6 

77.5 

67.8 

59.6 

52.2 

45.4 

1.7 

71.5 

60.2 

50.8 

42.7 

36.0 

1.8 

65.6 

53.0 

43.0 

34.6 

27.9 

1.9 

61.7 

46.2 

37.1 

26.2 

21.1 

2.0 

54.3 

39.9 

29.4 

22.7 

15.8 

6.4 

2.2 

44.  1 

29.0 

19.1 

12.8 

8.4 

2.4 

35.1 

20.7 

12.3 

7.3 

4.3 

2.6 

27.6 

14.4 

3.9 

2.1 

2.8 

21.5 

10.0 

2.2 

1.0 

3.0 

16.6 

6.7 

2.7 

1.4 

0.4 

1.4 

3.5 

8.5 

4.5 

4.0 

4.1 

1.8 

4.5 

1.9 

5.0 

0.8 

VITA 


.219 


VITA 

I-pai  Wu  was  born  in  Ching-kiang,  Kiang-su,  China,  on  June  23, 
1933.  He  finished  his  high  school  education  at  Tainan,  Taiwan,  China 
in  June  1951,  and  received  the  Bachelor  of  Science  degree  in  Agricul- 
tural Engineering  from  the  National  Taiwan  University,  Taipei,  Taiwan, 
China,  in  1955. 

After  his  graduation,  he  Joined  the  Shih-men  Development  Commis- 
sion as  an  assistant  engineer,  and  was  employed  in  the  dam  construction 
department . 

In  the  Spring  of  1958  he  enrolled  in  the  Department  of  Irrigation 
Engineering,  Utah  State  University,  Logan,  Utah.  In  the  Summer  of  the 
same  year,  he  transferred  to  the  Agricultural  Engineering  Department, 
Purdue  University,  West  Lafayette,  Indiana.  He  received  his  Master  of 
Science  degree  in  Agricultural  Engineering  in  February,  i960. 

He  started  his  graduate  study  in  Hydraulic  Engineering  in  Septem- 
ber, i960  in  the  School  of  Civil  Engineering,  Purdue  University,  since 
that  time,  he  has  been  working  on  the  small  watershed  research  sponsored 
by  the  Indiana  State  Highway  Department  and  administered  by  the  Joint 
Highway  Research  Project,  except  from  July  I96I  to  August  1962  when  he 
was  employed  as  a  hydraulic  engineer  at  the  Indiana  Flood  Control  and 
Water  Resources  Commission  while  his  Major  Professor,  Dr.  J.  W.  Delleur, 
was  on  sabbatical  leave  in  France. 


